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Resume 



L'ere de I'astronomie gravitationnelle a deja commence. La premiere generation 
d'interferometres terrestres a ete construite et a atteint la sensibilite experimentale 
prevue au niveau du projet. Pour I'instant, les ondes gravitationnelles n'ont pas ete 
dctcctccs, mais a un tel niveau dc sensibilite, mcme dcs rcsultats ncgatifs prcsentent 
un interet physique. Dans ce contexte, deux exemplcs sent trcs rcmarquablcs. 

Le resultat negatif d'une recherche visant a la detection d'ondes gravitationnelles 
d'origine stochastique a permis de deriver une contrainte comparable aux limites 
cosmologiques obtenues a partir de la nucleosynthese et du fonds diffus a micro- 
ondes. 

Le resultat negatif d'une recherche visant a la detection d'ondes gravitationnelles 
emises par I'etoile a neutrons situee dans le remanent de supernova Cassiopee A 
a permis, pour la premiere fois, de placer des limites superieures sur I'amplitude 
du mode d'oscillation de type «r» des etoiles a neutrons. 

II est tres important de souligner que ces resultats ont ete obtenus par le fait que 
les interferometres de premiere generation ont atteint leurs sensibilites prevues lors de 
la phase de projet ; en d'autres termes : il etait attendu que nous puissions explorer 
des regimes physiques bien precis d'une fagon propre aux techniques experimentales 
de I'interferometrie gravitationnelle. De meme, nous nous attendons aujourd'hui a ce 
que la deuxieme generation d'interferometres, en voie de realisation prochaine, puisse 
inaugurer I'epoquc de detection d'ondes gravitationnelles de fagon routiniere. Parmi 
les sources les plus probables qui seront explorees d'abord, il y a les systemes binaires 
d'objets compacts, comme les etoiles a neutrons et les trous noirs : en effet, remission 
d'ondes gravitationnelles est telle que Ic mouvement orbital descend en spiralc et Ics 
premiers signaux attendus a ctre dctcctcs sont emis par les deux objets lors de cctte 
phase. Ce mouvement en spiralc est caracterise par une vitesse relative v et une sepa- 
ration relative r, telles que la dynamique orbitale est non-relativiste : par consequent, 
ce regime peut etre decrit par une serie perturbative dans le petit parametre e d'ordre 
{v / cY ovl {mG N / r) , c etant la vitesse de la lumiere, la constante de Newton et m la 
masse typique des objets ; les corrections a I'ordre principale, representees par la dyna- 
mique newtonienne, sont classees selon les puissances du parametre e qu'ils contiennent 
et forment I'expansion dite post-newtonienne. En raison de sa nature perturbative, I'ex- 
pansion post-newtonienne peut aussi etre interpretee d'un point de vue de la theorie des 
champs et reformulcc en termes dc diagrammes de Feynman. Au cours des ces dernieres 
annees, I'application de ces instruments a la gravite classique a ete relancee grace a I'in- 
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troduction d'un approchc bascc sur Ic thcoric effective des champs. Dans cc contexte, 
I'etude d'un systeme binaire est traitee d'une fagon similaire an probleme des etats lies 
non-relativistes en electrodynamique et cromodynamique quantique : en raison de cette 
analogie, le traitement d'un systeme binaire d'objets compacts selon la theorie effective 
des champs a ete appele « theorie non-relativiste de la Relativite Generale». Dans ma 
these, j 'explore differentes utilisations des methodes typiques de la physique de parti- 
cules empruntees par la theorie non-relativiste de la Relativite Generale. Notamment 
j 'applique de facon etendue sa serie perturbative des diagrammes de Feynman. Outre 
leur convenance au niveau des calculs, la technique des diagrammes pent ctrc adop- 
tee pour faire des investigations autour des theories alternatives de la gravite et pour 
considerer les informations provenant des experiences d'un point de vue de la theorie 
des champs. Dans le premier travail original que je presente dans cette these, j'ai etu- 
die une extension de la theorie non-relativiste de la Relativite Generale au cas d'une 
theorie tenseur-scalaire de la gravitation. J'y ai traite la dynamique des sources de type 
ponctuel et a une dimension ; je me suis interesse en particulier a I'etude de la renor- 
malisation du tenseur energie- impulsion pour les deux types de sources. Le deuxieme 
et troisieme travail que je decris au cours de cette these se situent dans le contexte des 
tests de la gravitation. La Relativite Generale est une theorie fortement non-lineaire : 
jusqu'a present, elle a fait I'objet de plusieurs tests experimentaux et elle a ete confir- 
mee avec succes. Toutefois, ces tests ont explore des regimes dynamiques correspondant 
aux premiers ordres de I'cxpansion post-newtonienne seulement : en particulier, I'ordre 
le plus eleve ayant pu etre investigue est (v/c)^, qui rcpresente le tout premier terme du 
secteur radiatif cn Relativite Generale : a travers la detection des ondes radio emises 
par le pulsar de Hulse et Taylor, on a mesure ce terme et done prouve I'existence de la 
radiation gravitationnelle. Par ailleurs, au moyen de I'astronomie gravitationnelle, nous 
nous attendons a pouvoir tester la Relativite Generale jusqu'a {v/cY"^ ! Les effets qui 
correspondent a ces termes sont si faiblcs que seulement une technique dcdice comme 
I'interferometric gravitationnelle pcrmet de les dctecter. Concernant I'expansion post- 
newtonienne jusqu'a I'ordre {v/c)^, une grande quantite de travaux dans la litterature 
a ete consacree a comprendre quelles sont les caracteristiques de la Relativite Generale 
testees par une experience et comment ces caracteristiques different des predictions 
faites par des theories alternatives de la gravitation. En vue des premieres detections 
d'ondes gravitationnelles, il est done important d'cnvisagcr des cadres theoriques et 
phenomenologiques propres aux nouveaux regimes dynamiques. Au cours des deux 
derniers chapitres de cette these, je discute un exemple original d'un de ces cadres. A 
travers la theorie non-relativiste de la Relativite Generale, les non-linearites de la gra- 
vitation sont representees par une serie perturbative de diagrammes de Feynman dans 
lesquels les gravitons classiques interagissent avec les sources materielles ainsi qu'avec 
eux-memes. Nous avons choisi d'etudier les non-linearites de la Relativite Generale en 
designant les vertex d' auto-interaction gravitationnelle par des parametres : cela nous 
a permis, par exemple, de traduire la mesure de la diminution de la periode orbitale du 
pulsar de Hulse et Taylor en une contrainte sur le vertex a trois gravitons, contrainte 
precise au niveau de 0.1%. Cette attitude originale constitue un nouveau cadre pheno- 
menologique dans le contexte des parametrisations du regime radiatif et de champ fort 



de la Relativite Gencralc : un tel cadre est tres utile pour interpreter les experiences 
gravitationnelles presents et futures. 
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Introduction 



The era of gravitational-wave (GW) astronomy has aheady begun. First-generation 
ground-based interferometers have been buih and operated at their design sensitivities. 
No GWs have been detected yet but, at this level of sensitivity, even negative results 
have physical relevance. Most notable examples are the following ones. 

The negative result of a search for periodic GWs from the Crab pulsar [1] showed 
that no more than 4% of the energy loss of the pulsar is caused by the emission 
of GWs: this bound lies below the upper limit inferred from the spin-down rate 
of the Crab pulsar. 

The negative result of a search for stochastic backgrounds of GWs [2] allowed one 
to derive upper bounds comparable to limits of cosmological origin inferred from 
the nucleosynthesis [3] and from the cosmic microwave background 

The negative result of a search for periodic GWs from the neutron star in the super- 
nova remnant Cassiopeia A enabled one to put upper limits on the amplitude of 
r-mode oscillations of a neutron star for the first time [5j. 

It is important to stress that these scientific results have been obtained because first 
generation interferometers met their design sensitivities; in other words: it was expected 
that one would explore well determined physical regimes in a way which is unique to 
the experimental techniques of GW interferometry. Analogously, today it is expected 
that the soon-to-come second-generation ground-based interferometers will open the 
phase in which GWs are routinely detected. Among the most likely sources to be 
observed first are systems of compact binaries like neutron stars and/or black holes: 
in fact, the emission of GWs cause the orbital motion to inspiral and the first signals 
awaited are the waves emitted by binary systems in this phase. The inspiral motion 
is characterized by a relative velocity v and a relative separation r that make the 
orbital dynamics non-relativistic: this regime can then be described by a perturbative 
series in the small parameter e ~ {v/cY ~ "^GnI"^-, where c is the speed of light, 
Gtv is Newton's constant and m the mass of one of the objects; corrections to the 
leading order Newtonian dynamics are labeled by powers of e and form the so-called 
post-Newtonian (PN) expansion. 

Being perturbative in nature, the PN expansion can also be interpreted from a 
field-theoretical standpoint and reformulated in terms of Feynman diagrams. In the 
last few years the use of these instruments in classical gravity has revived thanks to 
the introduction of an effective field theory (EFT) approach [6]. In this context the 
dynamics of a binary system is studied in a way analogous to non-relativistic bound 



1 



2 



states in QED and QCD [7] : for such a reason the EFT framework of ref . [6] has been 
dubbed non-relativistic General Relativity (NRGR). The use of an EFT framework 
in the two-body problem in GR is motivated by the fact that the binary dynamics 
is a typical situation in which the disparity of physical scales allows one to separate 
the relevant degrees of freedom in subsets with decoupled dynamics. For example, 
because the extension of the radiating source is much smaller than the wavelength of 
gravitational perturbations, the modes corresponding to the internal dynamics of the 
compact object can be integrated out and the source is effectively described as a point 
particle. By means of NRGR, it has been possible to obtain previously un-computed 
spin contributions to the source multipoles at next-to- leading order [8]: this provides 
the last missing ingredient required to determine the phase evolution to 3PN accuracy 
including all spin effects. 

In this thesis I explore different uses of the EFT methods of NRGR, in particular 
of its perturbative series of Feynman diagrams. Beside being a convenient tool for 
calculations in the PN regime of GR, diagrammatic techniques can be employed to 
investigate alternative theories of gravity and to have a different point of view on the 
physical information that can be extracted from experiments. In the first original 
study I present in this thesis I have extended NRGR to the case of a scalar-tensor 
theory and treated the dynamics of point-like and string-like objects; notably, I have 
addressed the issue of the renormalization of the energy-momentum tensor for both 
types of sources. The second and third study I report in this thesis are set in the 
context of testing gravity. GR is a highly non-linear theory of gravity: so far it has 
been the subject of many experimental tests and it has passed them all with flying 
colors. However, these tests have probed dynamical regimes that correspond only 
to the first few orders of the PN approximation: notably, the highest such order is 
{v/cY and represents the very first term of the radiative sector in GR. The existence 
of this term has been assessed by timing Hulse- Taylor pulsar and has provided the 
fundamental result of the reality of gravitational radiation. In contrast, by means of 
GW astronomy, one aims at testing GR up to (f/c)^^! This corresponds to the very 
high 3.5 next-to-leading order in the PN expansion, which comprises a set of non- 
linear effects that are very peculiar to GR and cannot be probed by other means than 
gravity interferometers. For what concerns the orders up to (f/c)^, considerable work 
has been done in the past to understand the peculiar features of GR that are probed 
by an experiment and to confront them with the predictions of alternative theories of 
gravity. In view of the first GW detections, it is then relevant to envisage other testing 
frameworks that are appropriate to the new dynamical regime. The one explored in 
the last two chapters of this thesis constitutes a relevant example: here, by means of 
the NRGR description, gravity non-linearities are represented by a perturbative series 
of Feynman diagrams in which classical gravitons interact with matter sources and 
among themselves. Tagging the self-interaction vertices of gravitons with parameters, 
it is possible, for example, to translate the measure of the period decay of Hulse- Taylor 
pulsar in a constraint on the three-graviton vertex at the 0.1% level. What emerges 
from this attitude is a phenomenological framework that is suitable for parametrizing 
the strong- field/radiative regime of GR. 
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This thesis begins with an account of the theoretical and phenomenological frame- 
works that can be employed to test relativistic gravity according to the dynamical 
regime. Chap. [T] starts with a concrete example of alternative theories of gravity, 
scalar-tensor theories, which will be a reference point for many of the following dis- 
cussions. In the course of discussing the consistency of scalar-tensor theories with 
experiments, I introduce the various testing regimes; moreover, the derivation of the 
testable predictions of scalar-tensor theories is a first example of the use of diagram- 
matic techniques in classical gravity. These particle-physics tools are the leitmotiv of 
the present thesis. Afterwards, I move to a more detailed treatment of the regimes, 
each of which is characterized by a parametrized phenomenological framework. In the 
order, I will discuss: 

the IPN order accessible by Solar-System experiments and the parametrized post- 
Newtonian approach, 

the onset of the radiative dynamics of GR at 0{v/c)^ probed by binary pulsars and 
the parametrized post-Keplerian framework, 

the higher PN orders characterizing the GW phasing formula and one of the possible 
parametrizations for the latter. 

Chap. [1] constitutes the background material for what concerns the NRGR-based in- 
vestigations presented in this thesis. 

In Chap. O then, I discuss the EFT approach to inspiraling compact binaries: here 
I expose the relevant ingredients needed in order to present my original studies. The 
latter are the subject of the remaining chapters. 

In Chap. |3]I report the aforementioned investigation of mine in which, together with 
Riccardo Sturani, I extended NRGR to the case of a scalar-tensor theory: this 
allowed us to treat the dynamics of point-like and string-like objects and to 
address the issue of the renormalization of the energy-momentum tensor for both 
types of sources. 

Chap, m contains a study I have performed with the gravity group at the University of 
Geneva. In this work we expressed GR non-linearities in terms of the non-Abelian 
vertices of the theory and we derived the bounds that can be put on them by 
means of experiments of relativistic gravity: in such a way, we put forward a 
proposal for a parametrized field-theoretical framework for testing gravity. With 
binary pulsars timing data we found interesting constraints on the vertices at 
the first two perturbative levels; then, we assessed the possibility of deriving 
meaningful constraints with GW observations at a first qualitative level. 

In Chap. [5] I address the issue of bounding the vertices with GW tests in a more 
quantitative fashion: this is done from the point of view of parameter estima- 
tion and relies on a follow-up study that I have taken up in collaboration with 
Emanuele Berti at the University of Mississippi. In this chapter I treat some 
technical details of the calculations that are common with the investigation of 
Chap, m moreover, I compare our original parametrized field-theoretical frame- 
work for testing gravity with other approaches in the literature. In the course of 



discussion, I present different and complementary strategies of investigation tliat 
take advantage of the peculiarities of GW astronomy ... an era that has already 
begun. 



Chapter 1 



Theoretical and phenomenological 

frameworks for 

testing relativistic gravity 

What the numerous tests of relativistic gravity have probed so far is actually the 
post-Newtonian (PN) regime, i.e. the first orders of a perturbative series that can be 
used to approximate General Relativity (GR) as well as other theories of gravity; non- 
fuUy-relativistic systems are indeed describable by means of an expansion in a small 
parameter represented by two inter-related quantities: the velocities with which the 
objects move and the gravitational fields they generate. 

Until present, experiments conducted in the PN regime have confirmed GR and 
tightly constrained alternative theories. In the confrontation between theoretical pre- 
dictions and experimental tests, it is convenient to pursue two avenues: one which is 
theoretically minded, the other which is based on a phenomenological attitude. In the 
first case, one compares the observable predictions of GR against those of consistent 
alternative theories of gravity. In a phenomenological context, the attitude is more 
oriented towards understanding which features of GR as the theory of gravity are actu- 
ally tested by means of a specific experiment. "This is a bit similar to the well-known 
psycho-physiological fact that the best way to appreciate a nuance of color is to sur- 
round a given patch of color by other patches with slightly different colors" |9]. In the 
case of testing gravity, this is a very convenient way of proceeding. In fact, GR is a 
very peculiar theory where many symmetries imply that otherwise possible effects are 
strictly zero. For instance, GR requires that the inertial and gravitational mass coin- 
cide. Another example is the effacement of the internal structure of compact objects, 
i.e. the fact that in GR one can describe the external field of a self-gravitating body by 
means of a point particle approximation to a high degree of accuracy. These are only 
two examples of the many symmetries of GR that have been confirmed by approach- 
ing experiments with questions like "what would be the observable manifestation of 
a violation of Lorentz invariance?". In other words, one can entertain the possibility 
of violating GR symmetries even outside the domain of consistent alternative theories 
of gravity: for the purpose of testing a symmetry, one considers models, instead of 
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theories, that are directly built to exhibit a possible violation. 

Both the theoretical and the phenomenological attitudes make use of parameters. 
In the case of alternative theories, the parameters enter the couplings of the action 
functional; on the other hand, in the set of phenomenological models, the parameters 
stem from tagging different physical effects. While GR predicts a number for some ob- 
servable, both models and theories have this number multiplied by a parameter, whose 
range of values usually contains the one predicted by GR. Therefore, in the compar- 
ison with experiments, at times one may consider parameters coming from the PN 
expansion of a theory like GR, at times one may deal with parameters that belong to a 
model and have been introduced only to weigh a specific contribution to an observable. 
This last case applies to studies I have conducted personally, where one re-interprets 
gravity non-linearities in a particle-physics context and tries to constrain the values of 
the non-Abelian vertices of the theory. In this context, one makes use of instruments 
like Feynman diagrams, which were first applied to classical gravity in 1960 [10]. More 
recently a particle-physics attitude has been employed with the aim of testing gravity, 
notably through the confrontation of GR with a specific class of alternative theories. 
In refs. [HI [12] Damour and Esposito-Farese (DEF) have started a thorough investiga- 
tion of scalar-tensor theories and their observable predictions employing Feynman-like 
diagrams. I will begin this chapter by giving an account of the studies by DEF as an 
example of the theoretical attitude in testing gravity; moreover, discussing the adoption 
of diagrammatic techniques in gravity will serve as a basis for introducing my work in 
the subject. 

In describing the observable predictions of scalar-tensor theories, I will introduce 
the different testing regimes of gravity, defined according to the typical values of the 
velocities and the gravitational fields at play in a given system. A more detailed account 
of these dynamical regimes will be presented when I discuss the phenomenological 
frameworks employed in testing gravity. This discussion begins in Sec. 11.21 with Solar- 
System tests of the weak-field regime and the parametrized post-Newtonian (PPN) 
approach. In the Solar System gravity is so weak and velocities are so small that only 
the first PN order has been completely mapped out: the conclusion is that the slow- 
motion/weak- field limit of GR is fairly consistent with all the experiments performed 
so far. 

With the discovery and timing of a binary pulsar [13] , it has been possible to probe 
GR through 2.5 PN order, i.e. up to corrections of order {v^ /c') included. At this order, 
the dynamics of the system is made of both conservative and dissipative contributions: 
therefore, binary-pulsar tests are sensitive to both the strong-field and the radiative 
regime of the theory of gravity. I will present these tests in Sec. 11.31 where I discuss 
the parametrized post-Keplerian approach. Most notably, the detection of radio pulses 
from the Hulse- Taylor pulsar fT3] has allowed one to measure the decay of the system's 
orbital period, which is a loss of mechanical energy due to the emission of gravitational 
radiation. This constituted the first proof of the existence of gravitational waves (GWs) 
and confirmed the GR expectation for this effect at the 0.02% level: for these reasons, 
R. A. Hulse and J. H. Taylor have been awarded the Nobel Prize in 1993. An even 
more spectacular laboratory for testing GR is represented by the recently-discovered 
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double pulsar J0737— 3039 [H]. The uniqueness of this system is due to the fact that 
both constituent objects are neutron stars pulsating in the direction of the Earth. 
Besides being oriented favorably, this binary is so relativistic that it has been possible 
to measure more parameters than with Hulse- Taylor pulsar and in much less time. 

For the sake of testing the PN theory of gravity to orders higher than (f/c)^, one 
needs to detect objects that are more relativistic than binary pulsars, like inspiraling 
compact binaries: these are systems of neutron stars and/or black holes which are 
in the very last stages of their orbital motion. The high PN orders that govern the 
dynamics of the inspiral stem from the interplay of very peculiar non-linearities that 
are intrinsic to the specific theory of gravity. In this context, it is convenient to adopt 
the same point of view that originates from phenomenological frameworks like the 
PPN: here GR is thought of as one candidate in the space of parametrized theories 
of gravity and the confrontation among theories is done at the next-to- leading order. 
It is relevant to extend this phenomenological approach to the high PN orders of the 
GW regime, where strong-field and radiative effects come into play. I will conclude 
this chapter by discussing one such extension of the PPN: this approach stems from a 
parametrization of the GW phasing formula and has been studied in refs. fl5l [161 flT] . 
Another extension of the PPN framework is discussed in Chapters H] and [H this study 
is based on the aforementioned particle-physics attitude and constitutes part of the 
original work presented in this thesis. 

1.1 Tests of gravity through diagrammatic techni- 
ques - The case of scalar-tensor theories 

The first use of field theoretical diagrammatic techniques in classical gravity dates 
back to 1960 with the work of Bertotti and Plebanski [10] but a more systematic use 
had to wait until the 90's. It is at this moment that DEF have put forward a program 
to construct the PN limit of scalar-tensor theories in order to compare it with that of 
GR. This program has been pioneered in ref. [Ll\ and is well summarized in the recent 
review from Esposito-Farese [18] that we will follow in the presentation. 

To better address the issue of modified theories of gravity, it is useful to quickly 
remind which are the basis hypothesis of GR: 

• Matter fields are universally and minimally coupled to one single metric tensor g^j^^,. 
This enforces the Einstein Equivalence Principle (EEP), which embraces the so- 
called "universality of free-fall" together with local Lorentz invariance and local 
position invariance (see Will's review JT9J for a comprehensive discussion). 




• The metric gfj,,y propagates as a pure helicity-2 field, i.e. its kinetic term is given by 



1. In this chapter we adopt the "mostly plus" convention for the metric 77py=:diag(-, +,+,+). 




the Einstein-Hilbert action^ 




(1.1) 
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o = O material sources 

□ = • left hand side of the field equations 

□ ^ = propagator of the fields (Green function) 

O = • fields 

Figure 1.1: Diagrammatic representations of matter sources, fields and their propagat- 
ors. Figure taken from ref. [T8] . 



where is Newton's gravitational coupling constant, g is the determinant of 
the metric and R its Ricci scalar. 

The theory of GR can then be encoded in the following action: 

S = SEnlg^w] + Smatterli^i, Q^iu] (1-2) 

where the matter fields ipi should be taken to mean every non-metric field, therefore 
including bosonic fields. To envisage a theory which departs from GR, one can either 
change the matter sector or modify the gravitational one. In the first case, then, one 
would violate the weak equivalence principle, i.e. the universality of free fall: this 
is one of the best tested of GR, at the point that it can be taken to represent the 
foundation of every (classical) theory of gravity [19]. To consider deviations from GR 
that are not already tightly constrained from experiments, one can choose to focus on 
modifying the Einstein-Hilbert action (11. ip . This does not mean having more metric 
fields, but more gravitational interactions, i.e. matter can couple to other fields, for 
example scalars, and these helicity-0 fields participate in setting up gravity together 
with the Einstein helicity-2 field. Nevertheless, there is still one single metric field to 
which matter is minimally and universally coupled: of course, this metric field will not 
be given by g^y anymore but rather by g^^ = gf_iu[giiu,4>a], with (pa indicating one or 
more scalar fields. In full generality, one can represent the total set of gravitational 
fields (helicity-2 + helicity-0) with $ and the matter variables with a; the result is a 
compact reformulation of eq. (II. 2p : 

S = S4<l>]+S„,aUer[(r,<^] (1.3) 

The dynamics of the fields, as dictated by their equations of motion, can be conveni- 
ently represented through symbols like those of Fig. 11.11 by means of this "dictionary", 
one can expand the full gravity+matter action of eq. (II. 3p in the diagrammatic fashion 
of Fig. 11.21 This is similar to what one does in particle physics, where the perturbative 
expansion of interaction amplitudes is expressed in the form of Feynman diagrams [20j : 
we will extensively employ these diagrams starting from next chapter, where we present 
the EFT approach to binary systems. 

In their studies of scalar-tensor theories, DEF aim at deriving the action that de- 
scribes the motion of bodies subject to their mutual gravitational interaction. To 
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s= 




kinetic term 
of the fields 



+ '^o[a] + O • +y 

free bodies linear interaction ^ 

of matter and fields higher vertices 

Figure 1.2: Diagrammatic representation of the full action of the theory, eq. ( II .Sp . 
expanded in powers of the gravitational fields $, which comprise helicity-2 and helicity- 
modes and are described by black blobs. The first line corresponds to the field 
action and the second one to the matter action (describing notably the matter-field 
interaction). Figure taken from ref. [18]. 




Figure 1.3: Diagrammatic representation of the Fokker action (4), which depends only 
on matter sources (white blobs). The "dumbbell" diagram labelled "Newton" repres- 
ents the Newtonian interaction oc G^, together with all velocity-dependent relativistic 
corrections. The 3-blob diagrams labelled "1PM" represent first post-Minkowskian cor- 
rections, i.e. the 0{Gj^) PN terms as well as their full velocity dependence. The 4-blob 
diagrams labelled "2PM" represent second post-Minkowskian corrections 0{G%). Fig- 
ure taken from ref. [T8] . 

this purpose, DEF integrate out the full set of gravitational fields (helicity-2 + helici- 
ty-0) to obtain what they call a "Fokker" action, i.e. a functional which depends only 
on the matter variables: 

^FokkerM = 5$ [I>[(t]] + cS^atter [u, <I[ct]] , (1-4) 

where $[cr] denotes a solution of the field equation 6S/6^ = for given sources a, with 
S given by eq. (II. 3p . In diagrammatic terms, this action is represented by Fig. II. 3| 
where one has expanded the gravitational perturbations $ up to the second post- 
Minkowskian (PM) order, i.e. including terms of 0{Gj^) with respect to Minkowski 
space-time. 

In the EFT terminology of Chap. |2l the Fokker action of DEF is said to be obtained 
by integrating out the modes that are responsible for binding the system, the so- 
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called potential modes, as opposed to the radiation modes, which rather describe GWs. 
Another characteristic of the EFT approach is power counting, a generalized version of 
dimensional analysis: such a tool enables one to assess, ahead of any calculation, the 
order at which a given diagram will contribute to the perturbative series. For example, 
we will explicitly see in Chap. |2] how power counting explains why in calculating the 
A^-body action it is justified to neglect diagrams containing loops of gravitons. 

Coming back to Fig. 11.31 the post-Minkowskian[§ series needs to be further expanded 
in the small velocity parameter: this is because of the virial theorem, which makes terms 
of 0{Gi\i) appear at the same time as terms of ^) in the post-Newtonian expansion. 
Notably, to get Newton potential and its corrections in v"^, one has to take the diagrams 
labeled by "Newton" in the post-Minkowskian series of Fig. 11.31 and consider that the 
matter-gravity vertices depend on the velocity through the proper time of the particle 
(cfr. eq. fl232D l The diagrams of Fig. O labeled by "1PM" are 0{G%) and contribute 
at IPN through the zero-th order term of their velocity expansion. Summing up the 
relevant amplitudes gives a modified version of the Einstein-Infeld-Hoffmann (EIH) 
action [21], where the usual contributions from the helicity-2 sector appear together 
with those of the helicity-0 terms. At IPN, the deviations from the GR behavior 
produced by the scalars can be parametrized by two coefficients /3ppN and 7ppn, where 
the index "PPN" is due to the fact that they belong to the general framework of the 
parametrized post-Newtonian (PPN) approach. Historically, these parameters were 
introduced by Eddington [22], who was not considering a specific alternative theory 
with respect to GR but rather had a phenomenological attitude in mind. This idea is 
at the heart of the PPN approach that we will describe in more detail in Sec. 11.21 

Labeling the massive bodies with capital letters A, B, . . . , the modified EIH action 
reads: 



^1 



Fokkcr 



J2 I dtmA^^Jl--^\/c^ (1.5) 



1 \ ^ [ , GN'mA'tTlB 



I/O On 3 



1 / N/ N , 7PPN/ n2 

- l^y^AB ■ ^a)[t^AB ■ Vb) + -^(Va - Vij) 

where: 

- Tab denotes the instantaneous distance between bodies A and B; 

- hab is the unit 3- vector pointing from B io A; 

- is the 3-velocity of body A; 

- the sum over B ^ A ^ C allows B and C to be the same body. 



2. We revert for a while to the non-abbreviated versions of the PN and PM acronyms to better 
stress the difference between the post-Minkowskian and post-Newtonian expansions. 
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- the first hne describes free bodies in special relativity and corresponds to what 
has been noted by 5*0 [a] in Figs. 11.21 and [T73| 

- the second line describes the 2-body interaction, i.e. the dumbbell diagram of 
Fig. [L3l labelled "Newton": the lowest-order term of this diagram is the Newtonian 
potential, while the remaining terms are its first post-Newtonianian corrections 
of order 0(f^/c^); 

- the last line comes from the diagrams labelled "1PM" in Fig. II. 3[ at order zero 
in V because of their gravity non-linearities. 

The theoretical content that the parameters /3ppN and 7ppn have in the context of 
scalar-tensor theories can be conveniently addressed by means of a specific model. It 
turns out that the simplest possibility of a single scalar field is general enough to 
exhibit all the peculiarities of this type of theories. Let us then take for definiteness 
the action [23l[2il[25l 



S = — ^ / d^X^ {R - 2g''''d^^d,^) + Scatter Qf^u = e^^'^^^Q^,] , (1.6) 

where a potential V{(p) has been discarded on the assumption that the field is massless 
and that it is not quickly varying on the scale of the A^-body system. The field g^y is 
the physical metric, the one to which matter is universally coupled: it is the metric that 
defines the lengths and time intervals as measured by rods and clocks. As we anticip- 
ated, this metric depends on both the Einstein metric and the scalar field through the 
product of g^y and a function of the scalar e'^"'^'^\ This function is what characterizes 
the coupling between matter and the scalar field. To enable a perturbative study, it is 
convenient to expand the coupling function around the background value that the 
scalar field assumes far from any massive body: without loss of generality, this value 
can be assumed to vanish so that 

a{^) ~ ao^ + i/3o<^' + ■ ■ ■ (1.7) 

where defines the strength of the coupling between matter and scalar excitations 
at linear level, while [3q expresses a quadratic coupling where a massive body sources 
two scalars. The case (ao = 0, /3o = 0) represents GR, while the case (ao 7^ 0, /3o = 0) 
corresponds to another type of mono-scalar-tensor theory, that of Brans-Dicke [261 [27] : 
in this theory the coupling function between matter and the scalar field is completely 
hnear and characterized by a single parameter cjbd such that 2u-QYi + 3 = l/ctQ- 

As we will now see, the pair of parameters (ao? /^o) describes peculiar PN deviations 
of the class of mono-scalar-tensor theories under consideration with respect to GR. At 
Newtonian level, only ao contributes through a scalar analog of Newtonian potential: 
however, this effect is unobservable because it just redefines the gravitational coupling 
as an effective constant: G^^^ = Gn{1 + a^). It is only at next-to-leading order that 
scalar effects are no longer degenerate with the GR predictions. This could have been 
anticipated from eq. (11. 5p . where the PPN parameters are seen to modify only some 
terms at (9(1PN). If one indicates gravitons by curly lines and scalars by straight 
lines, the IPN order of the mono-scalar-tensor model is given by the diagrams of 
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, A B A B 

j- OTSnSSX) + J O O 

B C 




Figure 1.4: The IPN diagrams contributing to the modified EIH A^-body action (11. 5p 
in presence of a mono-scalar-tensor theory of gravity. Gravitons are represented by 
curly lines, while scalars correspond to straight lines. Figure taken from ref . |18] . 



Fig. II. 4[ where the left hand side column contains the usual GR terms and the velocity 
dependence of the vertices in the first line of diagrams has been left implicit (cfr. 
Fig. Em in Chap.E]). 

The series of diagrams of Fig. I1.4l has enabled DEF to calculate how and Pq enter 
the PPN parameters 7ppn, /3ppn introduced in eq. (11. 5p : 

^PPN = 1-— -2 , /3PPN = 1 + - 2 • (1-8) 

1 + "o 2 (1 + a^) 

The factor in 7ppn comes from the exchange of a scalar particle between two bodies: 
this is the case for the second diagram in the first line and for the third one in the second 
fine; the factor ao/^o^o in /^ppn comes from the purely scalar three-body process of the 
V-shaped diagram: here two bodies exchange one scalar with a third body that has 
a quadratic interaction weighted by (3q. Through eq. ( II. Sp . these scalar effects are 
linked to the terms that are multiplied by the PPN parameters in the modified EIH 
action (11. 5^ : therefore, experiments that probe the conservative dynamics of a system 
at IPN can be used to bound ao and /3o- In the case of the Solar System, these bounds 
are presented in Fig. II. 5| where, for symmetry reasons, only the upper ao plan is 
reported. The various tests to which the figure refers are: 

perihelion shift = measure of Mercury's perihelion advance [28] : 

LLR = Lunar Laser Ranging [29], which monitors the Earth- Moon distance to check 
if there is a difference between the Earth's and the Moon's accelerations towards 
the Sun; this due to the so-called Nordtvedt effect [30], which is absent in GR 
because of the strong equivalence principle (SEP): I will comment more on this 
effect in Sec. II. 2j 

VLBI = Very Long Baseline Interferometry [31j, an experimental technique that meas- 
ures the deflection of light in the curved spacetime of the Solar System; 

Cassini = measure of the time-delay variation from Earth to the Cassini spacecraft 
when this last was near solar conjunction [32] . 

As it is evident from the flgure, the most constraining experiment for what concerns 
the linear coupling constant is Cassini [32], which implies that ao must be smaller than 
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Figure 1.5: Solar-System constraints on the parameters (ao, f^o) defining a mono-scalar- 
tensor theory of gravity. The dark grey region is allowed by all the tests, which are 
explained in the text. Figure taken from ref. |18j . 



3 X 10^'^. On the contrary, the quadratic coupling /3o results basically unbound from 
Solar-System tests. The reason for this different behavior is due to the way and (3q 
enter the PPN parameters: /3ppN depends on the product between /3o and Oq, so that 
a small value of ao as constrained by 7ppn prevents (3o to be bound even in sign. 

Every point in the diagram corresponds to a theory of gravity: GR, which can be 
regarded as the very peculiar theory with vanishing values of the scalar parameters, 
sits in the origin of the axes. The full vertical axis /3o = corresponds to Brans-Dicke 
theory: because 2ubd + 3 = l/al, a bound on translates as cjbd > 4 x 10^ p2] . 

The horizontal axis corresponds to theories that DEF call "perturbatively equivalent 
to GR", which are characterized as follows. From the structure of the PPN parameters 
eq. (11. 8p and of the diagrams of Fig. II. 4[ one can foresee that the higher post-Newtonian 
orders of the scalar-tensor theory (ao,/3o) will involve at least two factors ao, so that 
schematically 

deviation from GR = x 

where the Aj's are constants built from ao, Po and the other parameters entering the 
coupling function a{ip) at higher orders. Because is already constrained to be small, 
one could expect the theory to be close to GR at any order. However, the scalar 
sector of the theory can exhibit non-perturbative effects [33] that make the square 
brackets of eq. (11. 9p become large enough to compensate even a vanishingly small Oq. 
These effects may occur in strong-field conditions and are analogous to the spontaneous 
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Figure 1.6: Heuristic argument to explain the phenomenon of spontaneous scalariza- 
tion. When (3q < and the compactness mGN/Rc^ of a body is large enough, it is 
energetically favorable for the object to create a local scalar field different from the 
background value. The body becomes thus strongly coupled to the scalar field. Figure 
taken from ref. 135). 



magnetization of ferromagnets [34]. In fact, if the so-called compactness mGN/R(? of 
a body is greater than a critical value, an object like a neutron star can become strongly 
coupled to the scalar field and acquire a charge with respect to it. A heuristic argument 
is enough to illustrate this peculiar phenomenon. As an example of a theory that is 
perturbatively equivalent to GR, one can take the limiting case of a model for which 
ao vanishes strictly, so that the coupling function is a parabola a((y9) = |/3oV^^- The 
behavior of the scalar field can be described by the usual 1/r fall-off outside the compact 
object and by the value that the field assumes at the center of the body. The energy 
of the scalar is made up of two contributions, one that comes from its kinetic term and 
one from its coupling with matter (see eq. (11. 6p for both): as a rough estimate of the 
energy configuration, one can write 

which, for negative values of /3o < 0, is the sum of a parabola and a Gaussian. Fig. 11.61 
shows plots of the function E{ipc) for three different values of the ratio m/R affecting 
the compactness: when this ratio is large enough, the function E{ipc) has the shape of a 
Mexican hat, which makes it energetically favorable for the star to create a nonvanishing 
scalar field (pc, and thereby a nonvanishing scalar charge given by a'(v?c) = Po'Pc- 

The heuristic argument has been confirmed by a full-fledged treatment of the physics 
involved in the phenomenon of spontaneous scalarization, which requires taking into 
account the coupled differential equations of the metric and the scalar field, plus using 
a realistic equation of state to describe nuclear matter inside a neutron star [33 1 13 ^ 136]. 
The main consequence is that the various PN-expanded predictions for the observables 
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Figure 1.7: Scalar charge ua versus baryonic mass m^, for the model a[}p) = — 3<y9^ 
(i.e., /3o = —6). The solid line corresponds to the maximum value of allowed by 
Solar-System tests and the dashed line to ao = 0. The dotted lines correspond to 
unstable configurations of the star. Figure taken from ref. |35] . 

depend on the scalar parameters in much the same way they did in the weak regime: 
the only difference is a replacement of the parameters (ao,/3o) by their strong- field 
versions, which are body-dependent 



As a consequence, the effective gravitational constant between two bodies A and B 
reads 



instead of G^^-^ = G]\f{l + Oq): the fact that there is now a body- dependence makes 
this change observable when one is dealing with different objects. In a similar way, 
the strong-field analogues of the PPN parameters 7ppn and /SppN (H.Sp involve a factor 
of a A for body A and a factor of for body B, instead of an ao for both bodies. 
Remarkably, the linear coupling is what plays the role of a scalar charge in the process of 
spontaneous scalarization: this takes place above a critical mass, whose value decreases 
as "— /3o" grows. On the other hand, if f3o is positive, strong field effects do not occur 
because the energy function E{{p) of eq. fll.lOp does not develop the essential Mexican- 
hat configuration. The behavior of the scalar charge is plotted in Fig. 11.71 for the 
particular model /3o = —6: it is seen that can be as high as 0.6, so that two factors 
of it can induce deviations from GR up to 36%. 

There is a major advantage in using systems made of compact objects like neutron 
stars as testing grounds for gravity with respect to the Solar System. In the former situ- 
ation, both the velocities and the gravitational fields at play are a factor of 10^ larger 
than our planetary playground: this gives access to higher orders in the PN expansion, 
notably to those which are odd in time and describe the dissipation of mechanical 
energy through the emission of gravitational radiation. Therefore, in the presence of 



a A = dlnrriA/dipo 
/3a = d a a/ d(po . 
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binary pulsars, one is testing gravity is in its strong-field and radiative regime. We 
will discuss the astrophysical laboratories represented by two-body systems compris- 
ing neutron stars in more detail in Sec. 11.31 where we present the phenomenological 
framework employed to extract physical information from them, i.e. the parametrized 
post-Keplerian approach. For the time being, it is enough to discuss how the emission 
of gravitational radiation changes in a scalar-tensor theory with respect to GR. This 
difference basically consists in the fact that in a scalar-tensor theory the system can lose 
its mechanical energy through an extra "decaying channel" represented by the helicity-0 
modes. The resulting change in the orbital period is a different function of time than 
in GR. The mechanical energy lost by the system at the orbital scale is received by an 
observer at infinity in the form of a fiux of GWs; in the case of a scalar-tensor theory, 
this fiux is schematically given by 

GWfiux = | Q^"^l^P"^% 



helicity-2 

fMonopole Dipole Quadrupole Q ( ^ \ K (\ 13) 



1 



hclicity-0 



where the subscript "helicity-2" refers to the GR prediction, while the curly brackets 
labeled "helicity-0" contain the extra contributions predicted in scalar-tensor theories. 
Relying on the different powers of 1/c, one would generically expect that the mono- 
polar and dipolar contributions of helicity-0 are much larger that the usual helicity-2 
quadrupole of gr|. A closer inspection of the terms reveals that this is not always the 
case. In fact, the structure of the scalar monopole is of the type 

Monopole _G \ d{mAaA) ^ 0(171303) ^ ^ f ^ ^ 1 ^ iH^) 
c c \ dt dt J J 

so that it reduces to 0{l/c^) if the stars are at equilibrium, dt{mAC(A) = 0, which is 
the case for all binary pulsars detected so far. On the other hand, a collapsing star 
would make this monopolar term huge. For what concerns the scalar dipole it has the 
form ^ 

Dipole _ G ( C/lniAmB \ , \2 , f ^ 



so that it can be efficiently excited only in presence of asymmetric systems, for which 
c^A 7^ this is the case of a pulsar-white dwarf binary: the pulsar's scalar charge 
aA may be of order unity, as we saw in Fig. II. 7[ whereas a white dwarf is a weakly 
self-gravitating body so its scalar charge is rather given by ao which is known to be 
small. 

Observations of binary pulsars and composite systems, like the pulsar-white dwarf 
binary, can be used to put upper bounds on the extra decaying channel represented by 



3. For what concerns the hehcity-0 quadrupole it can be neglected because it is roughly 0{al) with 
respect to the helicity-2 quadrupole of GR [37l [36] , 
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scalar waves. A similar analysis can be done also in the case of fiducial GW signals, like 
those awaited for detection at the LIGO [38], Virgo [39] and LISA [40] interferometers. 



A thorough investigation of the subject has been put forward by DEF and the most 
up-to-date constraints are reported in ref . [18] , from which Fig. 11.81 is taken. This 
comprehensive graphical panel has been obtained translating the bounds on the strong- 
field parameters aA, f^A as constraints on their weak-field analogues ao, (3o: these are 
indeed the parameters appearing in the action 11.61 through the coupling function (I1.7P . 
In such a way, one can compare the different probing power that various testing grounds 
have with respect to the scalar-tensor theory characterized by ao, /3o- The allowed 
regions lie below the various curves, except for the case of the closed curve "LIGO- 
Virgo NS-NS" for which the region below the curve is excluded. The grey region is 
consistent with all the tests: this region contains GR, which sits at — oo on the vertical 
axis because the scale of the plot is semi-logarithmic. 

From top to bottom, the new acronyms with respect to Fig. 11.51 are: 

SEP = Strong Equivalence Principle: this test refers to the bounds coming from a set 
of low-eccentricity pulsar-white dwarf binaries, which can best probe the effect 
being asymmetric systems; 

B1534+12 refers to the homonymous pulsar [IT] ; 

J0737-3039 is the double pulsar [M]; 

B1913+16 is Hulse- Taylor pulsar [T3] : 

J1141-6545 is a pulsar-white dwarf binary [12] : 

PSR-BH is the curve that corresponds to a possible observation of a pulsar-black-hole 
binary: it is supposed that this discovery is not made by means of GW detection 
but rather through the analysis of electromagnetic signals. 

The equation of state used to describe neutron stars is the polytropic described in 
ref. [33] because it can represent also the behavior of more realistic mo delsB 

From the figure it is interesting to remark that, for the purpose of constraining 
scalar-tensor theories, future GW detections seem not to be competitive with binary 
pulsar observations already at hand. Moreover, it is very surprising to note that the 
curve "PSR-BH" is more constraining than the one labeled "LISA NS-BH". However, 
both these conclusions strongly depend on the type of alternative theory chosen and on 
the parametrization used for describing its PN deviations from GR. I will comment on 
both these assumptions in what follows next, starting from the parametrization issue. 

From eq. (I1.7p one can see that the pair (ao, Po) parametrizes the first two orders 
of the coupling function between matter and the scalar field: as a consequence, the 
choice of the pair (ao,/3o) is adequate to confront scalar-tensor theories with GR only 
at the lowest radiative orders accessible with binary pulsars. In fact, these astrophysical 
laboratories have allowed one to probe the PN expansion only up to the 0{l/c^) terms 
in eq. fll.l3p . On the other hand, the aim of GW astronomy is to probe gravity at 
least up to (9(l/c^^)! In Sec. 11.41 I will describe a phenomenological approach to test 
the PN expansion up to this order; later on, in Chapters |4] and |5l I will describe an 



4. Private communication by Esposito-Farese. 
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Figure 1.8: Constraints on generic tensor-scalar theories imposed by tests of relativ- 
istic gravity in different dynamical regimes, notably: Solar-System experiments, classic 
binary-pulsar tests and future detections of inspiraling binaries with laser interferomet- 
ers (see text for explanation of individual acronyms and the significance of the various 
tests). The hatched region is allowed by all the tests. Figure taken from ref. [T8] . 



attempt, which is similar in spirit but stems from a field-theoretical approach. By 
means of these extended PPN frameworks one can have a clear view on how deeply 
various tests are able to probe gravity in its strong-field and radiative regime: as it 
will follow from the presentation of these frameworks, GW observations have a unique 
potential of detecting and constraining high-order gravity non-linearities. 

Concerning the specific choice of alternative theory, it should be stressed that the 
scalar field studied by DEF is not "cosmological" in the sense that it does not modify the 
dynamics of the universe as described by GR. On the other hand, in the last few years, 
cosmological models have been investigated where a non-GR behavior of the universe 
on Hubble scales can be ascribed to some scalar field: for example, one can characterize 
the DGP cosmological model [13] in terms of a scalar field called the Galileon [H]. In 
DGP, the luminosity distance has a dependence on the redshift that is not the same 
of GR and affects the propagation of GWs on cosmological scales. The space-based 
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interferometer LISA is being designed to detect coalescences up to redsliift 2; ~ 20 [15] , 
so it will probe the Galileon scalar-tensor theory in a way that is not possible with 
observations made with binary pulsars. Moreover, the frequency band of LISA will 
give access to sources in a mass range that is very different from that of binary pulsars: 
if one of the component object is a neutron star, the other one in the binary would be a 
IO^Mq black hole! This type of events are called extreme-mass-ratio inspirals and will 
allow one to map the space-time surrounding a black hole: in a such a way, it will be 
possible to test GR in a very peculiar way, through confirming the existence of black 
holes and ruling out alternative types of very compact objects [46]. Another example 
of the limitation of binary pulsar tests in bounding alternative theories is offered by 
Chern-Simons modification of GR [17|: if this is the correct theory of gravity, we could 
only perceive a deviation from GR by means of GW tests, as we will see in Sec. II. 4[ 

There is another point to be stressed about the conclusions on scalar-tensor theories 
derived by DEF. The analysis summarized by Fig. 11.81 assumed that the constraints 
put by GW tests on the parameters (ao,/3o) refer to detection of single events. This 
partially explains^ why the curve "PSR-BH" is more constraining than the one labeled 
"LISA NS-BH", i.e. why detecting a neutron star-black hole system by pulsar timing 
could give more information than can be done by analyzing the same system through 
the emitted GWs with the most sensitive space based interferometer. This comparison 
does not take into account one of the advantages of entering the era of GW astronomy, 
i.e. the fact that, already with advanced ground based interferometers, it will be 
possible to observe many coalescences per year: for the case of a neutron star-black 
hole system, the number of detections per year is likely to be ten but could reasonably 
be as high as three hundreds |48] (see Table ISTT] in Chap. [5]). This large statistics will 
bring crucial improvements in constraining parameters like (oq, Pq) for the following 
reasons. In fact, these parameters are peculiar of the theory of gravity and do not 
depend on the astrophysical system: therefore, they are extrinsic with respect to the 
emitting source. As a consequence, if events are detected, each one can be interpreted 
as providing an independent measurement of an extrinsic parameter, for example Pq: 
the total error on (3q would then be a fraction 1/a//V of the error coming from a single 
observation. Detection of many systems has a further advantage in measuring extrinsic 
parameters: as we have seen for tests with binary pulsars, there might be systems that 
are better suited for measuring one parameter instead of another. For example, the 
scalar charge of a neutron star is better measurable if one can detect a dipolar flux from 
an asymmetric system like a pulsar- white dwarf binary (see eq. (11.151) ). Moreover, the 
bound derived on a given parameter from one system might be used as a prior when 
testing a different parameter with another source: in this way, one should be able 
to disentangle the effects of various extrinsic parameters. A final note concerns the 
possibility to further increase the accuracy of GW experiments by cross-correlating the 
outputs of many interferometers, which enables one to dig signals that lie below the 
sensitivity of an individual instrument: this has proven crucial in obtaining the bound 
on the GW background of stochastic origin with the LIGO instruments [2] that I will 



5. The other reason has already been discussed before: it is the fact that one is hmiting the 
confrontation to the lowest radiative orders given by eq. (jl.l3|) . 
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discuss in Chap. |3l 



1.2 Probing the weak-field regime of gravity - 

The parametrized post-Newtonian formahsm 

In the previous section I have discussed scalar-tensor theories as an example of 
alternative theories of gravity. In the presentation we have seen that, at first PN order, 
the action governing the dynamics of scalar-tensor theories is similar to that of GR: 
eq. (11. 5p shows that the structure underlying both theories is the same, the differences 
being encoded in parameters that multiply the PN expansion terms of O {v"^ ^itlGn /r). 
This similarity is not peculiar to scalar-tensor theories, it is a feature of the PN regime of 
metric theories of gravity; for this reason, in eq. ( 11. 5p the parameters have been labeled 
with the index PPN, referring to a general framework that can be used to describe 
this regime, irrespective of the underlying theory. The use of parameters to describe 
the first PN limit of metric theories of gravity goes under the name of PPN approach. 
Among the proponents and major exponents of this framework is Clifford Will, who 
has treated the subject in two monographs on tests of gravity and confrontation among 
theories [271 [19]. I refer to these reviews for extensive discussions and relevant citations 
of the PPN; in this section I give a concise account of only those features that overlap 
with scalar-tensor theories and that are more useful for the rest of the thesis. 

The basic idea was formulated by Eddington [22J who wrote the Solar-System metric 
as the IPN limit of Schwarzschild metric in isotropic coordinates introducing some 
phenomenological parameters /3ppN and 7ppn in front of the different powers of the 
dimensionless ratio mGN/rc^: 

900 = -1 + 2^ - 2/3pp. (^) \o(^), (1.16a) 
9^, = 5., (l + 27PPN^) + O (^) (1.16b) 

9o^ = + o(^-^Y (l-16c) 

In this particular form of the metric the parameters /3ppN and 7ppn can be given a 
heuristic physical interpretation. From eq. (I1.16ap the parameter /3ppN can be regarded 
as a measure of the amount of non-linearity that a given theory of gravity puts into 
the goo component of the metric. On the other hand the parameter 7ppn measures the 
amount of space-time curvature produced by a body of mass m at a distance r from 
it: in fact, to IPN order, the spatial components of the Riemann tensor read p7j 

Rijkl = 37PPN — {rijnkSii + niUiSjk - riiUkSji - UjUiSik - '^SjkSu + ^5jfc5j7^ , 



where Ui = Xi/r. 
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The most general version of the PPN metric has a total of ten parameters: here 
I will only consider /3ppN and 7ppn because they are the only non-vanishing ones in 
scalar-tensor theories; within this choice, the metric reads [19] 



goo = -l + 2f/-2/3ppNf/' + (27ppN + 2 + a3)$i (1.17) 

-F2(37ppn - 2/3ppN + 1)$2 + 2<l>3 + 67ppn$4 + 0{e^) 
9o^ = -^(47PPN + 3)\/,-^W, + 0(e^/2) 

gij = (1 + 27ppn[/)(5,,- + (9(e2) , 

where e is the PN expansion parameter of order [v"^ ^tuGn /"r) and capital letters in- 
dicate metric potentials such as: 

U = ! , .. d?x' (1.18) 



X — x' 



in which p is the density of rest mass as measured in a local freely falling frame mo- 
mentarily comoving with the gravitating matter and = {dx^ /dt) is the coordinate 
velocity of the matter. 

The form (11.171) for the metric corresponds to the standard PPN gauge; this choice 
is only dictated by convenience and it relies on these assumptions: 

1. matter is modeled as a perfect-fluid; 

2. the spatial part of the metric is diagonal and isotropic; 

3. the metric is expressed in a quasi- Cartesian coordinate system at rest with respect 
to the Universe rest frame. 

With the explicit form (I1.17P one can obtain the equations of motion in the PPN metric, 
which are given by a set of coupled equations for matter and non-gravitational field 
variables in terms of other matter and non-gravitational field variables. For example, 
the geodesies of light rays can be written as [27] 

^ = (1 + 7ppn) m - 2h{n ■ VU)] (1.19) 
= -(l + 7PPN)t/, 

where Xp is the deviation of the photon's path from uniform straight line motion. 
These equations enable one to compute the deflection of light and the time delay of 
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an electromagnetic signal which are produced by a massive object; both these effects 
are characterized by the same factor (1 + 7ppn) that appears in eq. (ll.lQp . The best 
constraint on the parameter 7ppn comes from the measure of the time-delay variation 
from Earth to the Cassini spacecraft when the latter was near solar conjunction [32j 
and gives 7ppn ~ 2.3 x 10~^. 

In the case of massive bodies with spherical symmetry, the equations of motion 
have a quasi-Newtonian form 

a = (mp/m)Vt/ (1.20) 

where m is the inertial mass of the body and rrip its passive gravitational mass, i.e. the 
mass that determines the force on a body in a gravitational field: this mass given by 

nip = m — r]^Eg (1-21) 

where Eg is the negative of the gravitational self-energy of the body [E^ > 0) and 

?7n = 4/3ppN - 7PPN - 3 (1.22) 

is a parameter that measures the resulting violation of the massive-body equivalence 
principle, which is called the Nordtvedt effect. As we know from the discussion of 
Fig. 11.51 this effect is absent in GR but present in scalar-tensor theories. It refers 
to self-gravitating astronomical bodies, for which the ratio Eg/m may be significant: 
4.6 X 10~^° for the Earth and 0.2 x 10~^° for the Moon; therefore, if the Nordtvedt effect 
is present the Earth should fall toward the Sun with a slightly different acceleration 
than the Moon. This difference in accelerations affect the Earth-Moon distance that 
can be monitored with the lunar laser-ranging experiment (LLR): using a retroreflector 
left on the Moon by the most famous Apollo mission in August 1969, LLR has made 
regular measurements of the round-trip travel times of laser pulses sent from the Earth. 
The accuracies are at the level of 50 ps, which means that one is sensitive to errors 
of 1 cm in distance. At this level of precision, it is important to make sure that one 
is really probing the strong version of equivalence principle, the one valid for massive 
self-gravitating bodies. In fact, it could be that the Earth-Moon system exhibits a 
violation of the weak equivalence principle due to the different chemical compositions 
of the bodies. To this end, the world-leader group Eot-Wash has tested laboratory- 
sized bodies whose chemical compositions mimic that of the Earth and Moon. In 
this way, one could constrain the difference in accelerations due to a violation of the 
weak equivalence principle to 1.4 parts in 10^^ [19], which bounds the difference in 
accelerations due to a violation of the strong equivalence principle at the level of about 
2 parts in 10^^. For what concerns the Nordtvedt parameter, the resulting bound is 
|?7n| = (4.4 ± 4.5) X 10~^, which in turns constrains the parameter /3ppN to be smaller 
than 2.3 x 10"^ [29]. 

In Chapters m and [5] I will present a field-theoretical approach that allows a different 
interpretation of /3ppN- To conclude this section, in Fig. 11.91 1 report the constraints 
on the PPN parameters 7ppn and /3ppN imposed by Cassini, LLR other Solar-System 
experiments. 
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Figure 1.9: Constraints on the PPN parameters 7ppn and /Sppn imposed by Solar- 
System experiments. "Mars radar ranging" refers to measures of radar echo delay 
between the Earth and Mars as a probe of the propagation of light in the curved space- 
time of the Solar System. For the other tests see the discussion following eq. (11. Sp . The 
region allowed by all the tests is the intersection of the almost vertical strip labeled 
"Lunar Laser Ranging" and the horizontal bold line labeled "Cassini". GR sits in the 
point (/3ppN = 1,7ppn = !)• Figure taken from ref. [18]. 
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1.3 Probing the strong-field/radiative regime of 
gravity - Part I 

The parametrized post-Keplerian formahsm 

In this section I describe a phenomenological approach to testing gravity first pro- 
posed by Damour in ref. [50] and then extended by Damour and Taylor in ref. jSlj. 
The presentation follows the one recently given by Damour at the SIGRAV lectures [9]. 

The equations of motion of N compact bodies in GR can be derived from the 
following action 

S=Y^^ f d^+'x^R{g)-J2^aC I ^-g,,{z^)dzUzl, (1.23) 

a=l 

where the first piece is the usual Einstein-Hilbert action that accounts for the pure 
gravitational sector and the second term describes the coupling of a point-like source 
to gravity; the space dimension D is kept as an arbitrary complex number until the 
end of the calculation when it is set to its physical value D = 3: this serves to regular- 
ize the divergencies stemming from the point particle approximation. The equations 
derived from this action can be solved applying the weak-field PM expansion that we 
discussed in Sec. 11.11 Using this method, Damour and collaborators have derived the 
equations of motion of two compact bodies at the 2.5PN order [521 [531 iSl]: this level 
of accuracy is sufficient for describing the onset of dissipative contributions in binary 
pulsars' dynamics according to 

= ^loi^a - Zb) + C'^ A\^{Za - Zb, Va, Vb) 
+ A\^{Za - Zb, Va, Vb, Sa, Sb) 

+ c-^ Al,{za - Zb, Va - Vb) + 0{c-^) (1.24) 

where: 

Za denotes the coordinate of the pulsar a, Va its velocity and Sa its spin; 

Al^Q = —nibGN^zl — zl)/\za — Zb\^ denotes the Newtonian acceleration, i.e. the OPN 
order in the expansion; 

y4^2 is the term of order (9(f/c)^, i.e. the IPN correction; 

y4^4 stands for the 2PN modification which comprises spin-orbit effects; 

y4Jj5 is the term that had a crucial role in the discovery of gravitational radiation; it is 
the first term in the expansion which is odd in time and describes the decrease 
of mechanical energy of the system by an amount which, on average, equals the 
energy lost in a fiux of GWs at infinity. For this reason, the detection through 
radio pulses of the presence of in binary pulsar dynamics constituted a direct 
proof of the existence of GWs. 
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In order to investigate the pulsar dynamics through experiments, one needs to link the 
equations of motion (ll.24p to observational effects through two steps: 

1. solve the equations of motion fll.24p so to get the coordinate positions Zx and 
as explicit functions of the coordinate time t; 

2. relate the coordinate motion Za[t) to the pulsar observables, i.e. to the times of 
arrival of electromagnetic pulses on Earth. 

The first step has been accomplished in ref. [55], where the relativistic two-body mo- 
tion has been written in a simple quasi-Keplerian fashion. The second step has been 
accomplished through the derivation of a relativistic timing formula |56l [57] , a multi- 
parameter mathematical function that encodes many different physical phenomena: 

- dispersion effects; 

- travel time across the Solar System; 

- gravitational time delay due to the Sun and the planets; 

- time dilation effects between the time measured on the Earth and the Solar- 
System-barycenter time; 

- variations in the travel time between the binary pulsar and the Solar-System 
barycenter like those due to relative and proper motion; 

- time delays happening within the binary system. 

For the purpose of discussing the relativistic timing formula, one can concentrate on 
the time delays that take place within the binary system. These are the "aberration" 
effect and the delays named after Einstein, Roemer and Shapiro; I will discuss them in 
what follows. 

In describing pulsar timing, it is convenient to use a "multi-chart" approach which 
owes its name to the adoption of + 1 separate coordinate systems to account for the 
motion of bodies (see ref. [9] for more details and full list of citations). The space- 
time inside and around each body a is described by local coordinate charts denoted 
by X", where a = 0, 1, 2, 3 and a = 1, 2, . . . , A^. The spacetime outside the A^ "tubes" 
that contain the bodies is described by one global coordinate chart denoted by x'^ with 
H = 0,1, 2, 3. The different coordinate systems are related among themselves through 
expansions of the form 



where zj^ denotes the global coordinates of the center of mass of the pulsar, Ta the local 
(proper) time of the pulsar frame and where, for instance. 



By means of the multi-chart decomposition, it is possible to combine the information 
contained in several expansions: in the global frame one can use a weak-field expansion 




1 



(1.25) 




(1.26) 



(1.27) 
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where the indices (0), (1) in the exponents correspond to powers of Newton's constant; 
in the local frames one can use expansions of the type 

G^p{X2) = G%{X:- rria) + ^^^(XJ; m„ m,) + ■ ■ • , (1.28) 

where G^^^{X;ma) denotes the (possibly strong-field) metric generated by an isolated 
body of mass m^. The separate expansions (11.271) and (ll.28p are then matched in some 
overlapping domain of common validity of the type 

T^a Gn/c^ ^ Ra ^ \x — Za\ <^ d \xa — Xh] with b ^ a. (1-29) 

With this setting in hand, it is possible to relate the time delays that happen within 
the binary pulsar to quantities which are observable in the Earth system. 

In the rest frame attached to the pulsar a (X° = cTa,X*) the pulsar phenomenon 
can be modelled by the secularly changing rotation of a beam of radio waves: 

*.^/«.ra.T..n.T.4n.T.^ + .... (1.30) 

where is the longitude around the spin axis. Rest frame quantities like Ta and 
have to be translated into the global coordinates that describe the dynamics of the 
binary system; this is done by using the link (ll.25p . Among other results, one finds 
that a radio beam emitted in the proper direction in the local frame appears to 
propagate, in the global frame, in the coordinate direction given by 

n' = N' + --N' —— + 0i^] : (1.31) 
c c V c / 

this change in direction is called the aberration effect. 

The Einstein time delay stems from the link between the pulsar proper time Ta and 
the global coordinate time t = x^/c = z^/c used in the orbital motion; to first PN order 
{0{Gj\f, f ^)), this time link is given by 



i.e. by a sum of the special relativistic and general relativistic time dilation effects: the 
combined effect of this terms clearly deserves the name of Einstein time delay [57] . 

The Roemer and Shapiro time delays rather originate when one computes the global 
time needed to reach the Solar-System barycenter by a light beam emitted by the pulsar; 
this is achieved by imposing that this light beam follows a null geodesic according to: 

= rfs^ = g^,u{x^) dx'' dx" ~ - (^1 - dt^ ^ ^ 5") ' ^^'^^^ 

where U = maGN/\x — Za\ + mbGN/\x — Zb\ is the Newtonian potential within the 
binary system. Indicating by te the time of emission and by ta the time of arrival, the 
interval between the two reads 

t^-t,= ['• M . 1 M.| + i /*" + ^\ \M . (1.34) 

Ju CJt. Jt. \ X-Za X-ZbJ 
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The result of the first integral is (1/c) \zssB{ta) — Za(te)\, i.e. the light crossing time 
between the Solar-System barycenter (SSB) and the pulsar: the fact that the pulsar 
moves on an orbit Za{te) gives rise to the so-called Roemer time delay. The second 
integral rather accounts for the fact that the beam is propagating in the curved space- 
time generated by the companion: this effect is called the Shapiro time delay. 

Once expressed as functions that relate the pulsar proper time to the Solar-System- 
barycenter time, the delays discussed so far will contribute to the GR timing formula 
as [57| 

tssB - to = D-'[T + A^(T) + Aij(T) + A^j(r) + AsiT)] , (1.35) 

where, here and henceforth, T denotes the pulsar proper time, D is an overall Doppler 
factor and the terms denoted by Aj's encode the contributions from the various physical 
effects: A for the aberration, E for Einstein, R for Roemer and S for Shapiro. The 
structure of the timing formula fll.35p can be conveniently expressed in a condensed 
parametrized fashion 

tssB -to = F [T^; {p^}; {p^^}; {q^^}] , (1.36) 

where: 

T/v is the pulsar proper time corresponding to the A^*'^ pulse, which is related to the 
integer N by an equation of the form 

Npuise = const + UpT +^UpT^ + ■ ■ ■ , (1-37) 

in which Up and Up are some spindown parameters; 
{p^} is the set of Keplerian parameters, like the orbital period Pb ; 
jpPXj -g q£ separately measurable post-Keplerianpaxameteis, like the derivative 

of the orbital period Pf, ; 

{q^^} is the set of non separately measurable post- Keplerian parameters, like the 
overall Doppler factor D ; being non separately measurable, these parameters 
can be absorbed into changes of the other parameters and set to some arbitrary 
values. 

The relativistic timing formula of eq. f ll.35p was derived by Damour and Deruelle [57] 
for GR but its validity is more general: remarkably, the mathematical form of its para- 
metrized expression (11.361) was found to be applicable also to scalar-tensor theories. As 
we saw in Sec. II. H these theories predict that binary pulsars observables differ markedly 
from their behavior in GR. Notwithstanding these differences, the parametrized form 
of the timing formula (I1.36P holds in the case of scalar-tensor theories too, thus en- 
abling one to treat binary pulsars observations by means of a single phenomenological 
approach: this is the so-called parametrized post-Keplerian (PPK) framework ^501 [51] 
and consists in the following steps. By a least-square fit of the observed arrival times 
to the parametrized formula (ll.36p . one can extract from raw data the best fit values 
of both the Keplerian parameters and a subset of post-Keplerian (PK) parameters: as 
we will see below, three PK parameters could be measured for the Hulse- Taylor pulsar 
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Figure 1.10: Plot of the parameters 7 and u for the Hulse- Taylor pulsar as functions of 
the constituent masses. The errors on the parameters are smaller than the thickness of 
the respective curves. The meeting point of the two curves defines the measured values 
for the pulsar and its companion. Figure taken from ref. [58] . 



while, in the case of the double pulsar, one has obtained six PK parameters and a 
measure of the mass ratio. The parameter extraction is phenomenological in that it 
does not rely on the validity of a specific theory of gravity. However, each theory of 
gravity makes its peculiar predictions about how the PK parameters are related to the 
Keplerian ones and to the masses of the binary system. It is through these links that 
one can test the consistency of a theory of gravity with experiments. The masses are 
a priori unknown and have to be inferred from the measured values of the parameters. 
Because every PK parameter p is a function of the binary masses and m^, it is 
possible to invert this function to obtain a relation of the type = fP[mh], i.e. a 
link between the values of and those of which are allowed by the parameter p. 
These links correspond to curves in the plane of the binary masses: the graphs of two 
such curves overlap, within errors, in a common region, which gives the values of the 
binary masses together with their statistical errors. The situation is shown in Fig. I1.10[ 
where two PK parameters are plotted for the case of Hulse- Taylor binary pulsar . The 
parameter 7 measures the amplitude of the Einstein time delay, A^;: it has the di- 
mension of time and should not be confused with the dimensionless PPN parameter 
7ppn probed by Solar-System experiments. The parameter u measures the periastron 
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advance per orbit. Both 7 and u are determined so precisely that the thickness of the 
respective curves is larger than the errors on the parameters; the inferred values for 
the masses are indeed very accurate [58]: Mp = (1.4408 ± O.OOO3)M0 for the pulsar, 
Mc = (1.3873 ± O.OOO3)M0 for the companion. 

Once two curves (parameters) have been used to fix the values of the masses, a third 
one should consistently overlap in the same common region of the masses plane: the 
theoretical prediction for the third parameter can be regarded as a check of the con- 
sistency between the theory under examination and experiments. As a third parameter 
let us take the time derivative of the orbital period Pf, , which is a direct consequence 
of the term A5 ~ 0{v/c)^ in the equations of motion for a binary system ( 11. 24^ : 
therefore, any test involving the observable Pb will be a mixed type of test of the 
radiative/strong-field regime. On the other hand, the two parameters u and 7 do not 
involve radiative effects. In the context of GR, the predictions for these parameters 
have been worked out in refs. [561 EZl [59]. Considering the fractional periastron ad- 
vance per orbit k = A9/2n = {Cj)/n = (u) Pb/27i, where n = 2ti / Pb indicates the 
orbital frequency, the GR predictions for /t, 7 and Pb read 

k'^^ma.mb) = Pi , (1.38) 

1 — 

7^^(m„ mb) = - Xb{l + Xb) I3l , (1.39) 

n 

109-^ 1 _L 73 2 I 37 4 

where e is the eccentricity of the orbit and 

M = ma + mb (1.41) 

Xa = nia/M ■ Xb = rrib/M ; X^ + Xb = 1 (1.42) 

,„(M).(^l^)". ,1.43) 

It is interesting to look at the scaling of /3o so as to compare the PN order of the 
various observables. By means of Kepler's law, = m/a^ where a is the semi-major 
axis of the orbit; this enables to write /3q ^ (M n)"^^^ = m/a, which, a part for a factor 
of Gat, is of order e . Therefore, k and 7 are seen to be 0{e) corrections to the Keplerian 
motion, while Pb is (9(e)^/^. 

In scalar-tensor theories the PK parameters k, 7 and Pb will not look the same 
as eqs. fll.38p - fll.40p : first of all, they will depend on the coupling constants of the 
theory; secondly, they will have a different functional dependence on the masses and 
on Keplerian parameters like the eccentricity and the orbital frequency. The expressions 
for the PK parameters k, 7 and Pb in scalar-tensor theories are reported in ref. [9] (to 
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Figure 1.11: A confrontation of the consistency tests reahzed by means of post- 
Keplerian parameters as measured in the Hulse- Taylor pulsar, with niA the pulsar's 
mass, ruB the companion's. The panel on the left shows that GR passes the test with 
flying colors, the panel on the right indicates that a scalar-tensor theory with Pq = —6 
is ruled out. The widths of the lines are larger than la error bars. Figure taken from 
ref. M. 
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where C^J^ is the effective gravitational constant of eq. (11.12^ . which comprises scalar 
effects, and where for the period decay we have indicated only the change in the quad- 
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Figure 1.12: Orbital decay of PSR B1913+16. The data points indicate the observed 
change in the epoch of periastron with date while the parabola illustrates the theoretic- 
ally expected change in epoch for a system emitting gravitational radiation, according 
to GR. Figure taken from ref. |60] . 



rupole contribution of helicity-2. Taking again Hulse- Taylor binary pulsar, a confront- 
ation of the consistency tests realized by means of these PK parameters is reported in 
Fig. 11.111 the panel on the left shows that GR passes the test with flying colors, the 
panel on the right indicates that a scalar-tensor theory with Pq = —6 is ruled out. 

Concerning Hulse- Taylor pulsar, it is interesting to report a test of GR realized by 
means of the periastron precession. In a way similar to the one used to determine the 
time of the A^-th pulse emitted by the pulsar (ll.37p . one can write a relation for the 
time of periastron passage Tp 

Np,ri = nTp + ^hT^ (1.47) 

where higher derivatives of the orbital frequency n have been discarded because of their 
smallness in Hulse- Taylor binary. The time of periastron passage Tp will then differ 



32 



1.3 The parametrized post-Keplerian formalism 



from the ratio Nperi/n as 




(1.48) 




n 



which is a parabola with coefficient —h/2n < 0; using the GR prediction for h, this 
parabola lies on top of experimental data with the remarkable agreement shown in 



An even more spectacular laboratory for testing GR through PK parameters is 
represented by the double pulsar J0737— 3039, whose probing power for scalar-tensor 
theories has been reported in Fig. 11.81 The peculiarity of this system relies in the fact 
that both the component neutron stars are pulsating in the direction of the Earth, thus 
making it unique so far. Besides being oriented favorably, this binary is so relativistic 
that one could measure more PK parameters than in the Hulse- Taylor case and in much 
less time. Notably, on top of the triple (7, u, Ph) it has been possible to determine three 
other parameters. The ffist two are the shape and range of the Shapiro time delay: 
with an obvious choice of symbols, they are indicated by s and r. The third parameter 
is the rate of spin-precession Qso caused by relativistic spin-orbit coupling. Moreover, 
the fact that both neutron stars are detectable as pulsars has enabled one to measure 
the orbit of both objects around the common center of mass: in other words, if i is the 
inclination angle of the orbit with respect to the line of sight, the projected semi- axes 
xj = ajsini is known for both objects. The ratio of the two semi-axis is identical 
to the inverse ratio of masses, which is theory independent. Therefore, by measuring 
the orbits of the two pulsars relative to the centre of mass, one obtains a precise 
measurement of the mass ratio parameter R = 171^/^3- The tests corresponding to 
all these parameters are reported in Fig. I1.13[ where the curves in the masses plane 
overlap in the region allowed by the geometrical constraint sini < 1. As shown by the 
figure, GR passes all the tests corresponding to the measured parameters. 



Fig.ini 



Theoretical and phenomenological frameworks for testing gravity 



33 



T 1 1 ^— 1 , 1 T^— 1 , T 1 , 1 p-p 1 , r 




Mass A (MgyJ 

Figure 1.13: Mass-mass diagram of the double pulsar system summarizing the meas- 
ured PK parameters in combination with the derived mass ratio R = mA/TriB = 
Xb/Xa (solid red line) and constraints given by the mass functions of the binary sys- 
tem. The latter are indicated by the coloured regions which mark areas in the diagram 
that are excluded by the Keplerian mass functions of the two pulsars and the condition 
that sinz < 1. Further constraints are shown as pairs of lines enclosing permitted re- 
gions as predicted by GR. These are: the measurement of the advance of periastron Cj 
(diagonal dashed line); the measurement of the gravitational redshift/time dilation 
parameter 7 (dot-dash line); the measurement of the Shapiro parameters r (solid green 
line) and s (dotted green line); the measurement of the orbital decay Ph (dot-dot-dot- 
dash line) and the rate of spin precession of B, VLso (almost vertical dashed line). The 
inset is an enlarged view of the small square which encompasses the intersection of the 
tightest constraints. The permitted regions are those between the pairs of parallel lines 
and we see that an area exists which is compatible with all constraints, delineated by 
the solid blue region. Figure taken from ref. [6T] . 
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1.4 Probing the strong-field/radiative regime of 
gravity - Part II 

Extending the parametrized post-Newtonian 
framework to the gravitational-wave phasing 
formula 

In the PPN framework a phenomenological attitude is applied towards the first PN 
order in the expansion of the gravitational field: all metric theories share a common 
structure that is conveniently represented by means of parameters. The same philo- 
sophy can be applied to the mathematical expressions that describe the GW waveform. 
In recent years, indeed, an extension of the PPN framework has been proposed where 
the coefficients to be bound by experiments are those of the phasing formula of a GW 
signal. This line of investigation has been pursued in a series of three papers from 
Arun, Iyer, Sathyaprakash and collaborators [151 UHl EI]- I describe this framework 
in the present section, while in Chap. H] I will present an example of a study which is 
similar in spirit but has a field-theoretical interpretation. 

Because interferometers are most sensitive to the phase evolution of a GW signal, 
one usually neglects the PN corrections in the amplitude with respect to those of 
the phase [62] and works in the so-called restricted post-Newtonian scheme. This 
is one of the assumptions that can be made for both simplifying the treatment and 
demonstrating the feasibility of the tests themselves. Other simplifications stem from 
neglecting the spin of the objects and the eccentricity of the orbits. 

To obtain the phase, one can proceed as follows. From eq. fll.24p we know that the 
equations of motion for a binary system contain terms that are odd in time and therefore 
describe dissipative effects in the dynamics; notably, the first term that corresponds to 
the loss of the system's mechanical energy is ~ (v/c)^: once averaged over time, this 
loss equals the leading order fiux of GWs at infinity. As a consequence, the resulting 
change in the orbital phase of the binary system can be computed from the energy 
balance equation 

_^^!e_llZ£).^W, (1.49) 

where E is mechanical energy of the system and J-" the averaged fiux of emitted grav- 
itational radiation. The time arguments of the two quantities indicate that the energy 
carried by GWs at a time t and a distance r from the system are balanced by the 
system's loss at retarded time t — r/c. Therefore, eq. (11.491) expresses the conserva- 
tion of energy and a similar relation holds for the angular momentum. As we will see 
below, when one considers PN orders higher than (v/c)^, non-linear phenomena affect 
the propagation of GWs in a very non-trivial way: for example, the back-scatter of ra- 
diation produces the so-called tails of the wave that have support also inside the light 
cone and therefore propagate with an effective speed smaller than c. In this context 
it is not obvious a priori that the flux at time t is exactly compensated by a loss at 
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retarded time t — rjc. In the case of GR, eq. f ll.49p has been checked to be vahd up to 
relative 1.5PN order [63], while for alternative theories, like the scalar-tensor theories 
studied by DEF, the relation (11.49^ is usually assumed. 

The PN expressions for the mechanical energy and the GW flux are presently known, 
in the case of GR, up to order or 3.5PN [M] and are given by 



1 ^ 

k=0 

32 9 1 n V ^ ^ I. 1712 



T = V J-,^^ - ^ ln(.;) (1.50) 

^-^ 105 



k=0 



where u = 17111712/ M"^ is the symmetric mass ratio of the system and the explicit 
expression of the PN series coefficients Ek and J-^ as functions of u can be found in 
ref. [65]. Expressing E and J-" as functions of the velocity v, the energy balance equation 
transforms into the following coupled ordinary differential equations 

dv dE/dt —J-' 



dt dE/dv E'{v) 

-r = ^ = ITT (1-51) 
dt M ^ ^ 

where E'{v) = dE/dv and (p is the orbital phase. Finally, the GR prediction for 
the phase (pcw of gravitational radiation is obtained, at leading order, by two times 
the orbital phase, fGw{t) = 20(t): for example, in the extreme- mass-ratio limit its 
expression reads 

* ^ ,^dE/dv' 



(PGw{t)= druGwir) = ——— dv'v' , (1.52) 

J to (jN^-y^H Jv{t) J' \y ) 

where ojgw is the GW frequency and Mu the mass of the heaviest object in the system. 

From eq. (11.511) . we can see that the terms of PN phasing formula are combina- 
tions of the PN coefficients of E and T (ll.50p : therefore, each PN term in the phase 
encodes non-linearities that are intrinsic to GR. One such peculiar non-linearity is the 
aforementioned phenomenon of GW tails [661 EH EH] , which arises starting from order 
1.5PN in the phase. This effect is due to the fact that waves which are sourced by 
the quadrupole moment of the system can scatter off the curved (Schwarzschild) back- 
ground generated by the source: the scattering effectively delays these waves, which 
therefore propagate also inside the light cone (a Feynman diagram describing this effect 



is reported in Fig. 2.21(a) of Chap. [2]). Other viable theories of gravity have their own 
predictions for the PN coefficients of the phasing formula, in analogy with the altern- 
ative predictions for the PK parameters that I have discussed in Sec. 11.31 Through 
the measure of GWs, an accurate scrutiny of the PN coefficients will be possible; as 
in the case of the PPN framework, this will serve a two-fold purpose: check the con- 
sistency of GR itself and provide bounds on the parameters characterizing alternative 
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theories. Therefore, it is relevant to consider how to extend the PPN framework to the 
strong- field/radiative regime probed by means of GWs. Such an extension has been 
put forward in refs. [15l [161 IE] and consists in regarding the PN coefficients of the GR 
phasing formula as parameters to be estimated. To illustrate the procedure, we report 
the expressions for the waveforms that are valid in GR and then present the possible 
tests on its PN coefficients. 

In the restricted PN approximation, the response of an interferometric antenna to 
the incident radiation from a source at a luminosity distance Dl is 

h{t) = ^ [7:MF{t)f" cos [^Gw{t)] , (1.53) 

where Ai = u'^^^M is the so-called chirp mass of the system and F{t) is the instantan- 
eous frequency of radiation given by 

^^^^ = 2^^^- 

For the purpose of testing the PN phasing coefficients it is useful to work with the 
Fourier transform of the signal h{f) = J^^h{t) exp{27iift)dt. Using the stationary 
phase approximation it has been shown [69l [70] that one can write the waveform in the 
frequency domain as 

Hf) = Af-'/'exp [^vl>(/) +^^] , (1.54) 
with the Fourier amplitude A and phase "^{f) defined by 



5DL7r2/3 V 24 ' 

7 

^if) = 27r/t, - + + ln/]/(^-^)/3 . (1.55) 

j=0 



In the expression for the phase, tc and $c are integration constants that represent 
respectively a fiducial epoch of merger and the phase of the signal at that epoch. 
Moreover, the ^/'-coefficients of the PN expansion of the Fourier phase are given by 



(1.56) 
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and the a-coefficients are [7T] 

3715 55 

ao = 1 , ai = , a2 = — — + —i' , as = -167r , (1.57) 
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where the constant appearing in the expression for is Euler's constant, •je — 0.577. 

There is a total of nine PN parameters: seven are the coefficients of 0{v"') terms 
for n = 0,2,3,4,5,6,7 and two are coefficients of 0{v"'lnv) terms for n = 5,6. We 
stress the absence of tpi, the coefficient multiplying the term v^: this feature is peculiar 
to GR as we will discuss also in Sec. 15.3.31 On the other hand, the coefficient ip^ gives 
a contribution to the phase (I1.55P that is constant as it multiplies /°: for this reason, 
the coefficient ip^ cannot be used as a test parameter and only redefines the phase of 
coalescence. 

The PN coefficients of eqs. (I1.56P and (I1.57P have been quoted entirely to show 
explicitly that, in the absence of spin, GR predicts that they depend only on the masses 
of the systemQ through the parameters M and u. If the PN limit of GR is the correct 
description of the binary dynamics, the parameters of eqs. (I1.56P and (I1.57P must be 
consistent with each other within their respective error bars, much in the same way 
as the PK parameters inferred by binary pulsars observations (see Sec. ll.3p H. When 
talking about the set of the phasing parameters, we will collectively denote them by 
ipj^s, i.e. we will not use a different notation for the ipji^s; similarly, the error bars will 
be denoted by Aipj. 



6. An obvious exception to this is represented by the coefficients which are constants, hke a^. 

7. As suggested in ref. |15) . there are also non-perturbative schemes that one could implement 
to check the consistency of GR using GWs: for example, one could test the full numerical relativity 
predictions. Another possibility is offered by Pade-approximants [72] and the effective one-body 
approach |73l 1741 [75] : here, by means of re-summation techniques, one builds model waveforms that 
include both the inspiral and the merge. Implementation of these schemes to test gravity has not 
appeared in the literature yet. 
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Notwithstanding the proven consistency of GR, the presence of the statistical errors 
leaves room for alternative theories. With respect to GR, the phasing coefficients of 
alternative theories have a different behavior: either their functional dependence on 
the constituents masses is not the same as eqs. (11.56^ and (11.57^ . or they depend on 
other types of parameters. We already met an example of the latter category when we 
considered the scalar-tensor theory of DEF in Sec. 11.11 in fact, in eqs. f ll.13p -f ll.15p we 
saw that the parameter a a-, describing the linear coupling of a scalar field to matter 
in the strong-field regime, causes a change in the period decay of a binary pulsar at 
orders (f/c)^ and (f/c)^. In the language of the present section, these changes would 
bring two extra ■j/'/s of "negative" order with respect to the Oivjcf GR quadrupole: 
'0-2 for the monopolar term of Oiyjc)^ and for the dipolar term of 0(w/c)^. In 
ref. [15] the focus was rather towards a more phenomenological alternative in which 
the graviton has a mass. This possibility had been previously entertained by Will [76] , 
who begins by writing the velocity a graviton would have in a local inertial frame 

2 2 4 

where is the rest mass, E the energy. This velocity difference could be tested 
by monitoring a source that emits both gravitons and photons: an example could 
be supernova located at a distance D ~ zc/Hq, with z « 1 the redshift and Hq 
the Hubble constant. With the appropriate detectors on Earth, we would measure a 
difference in the times of arrival of the two types of radiation 



where the index "a" refers to the arrival times, the index "e" refers to the emission 
times. Such a time difference is related to the velocity difference by 



where a typical value has been chosen for D. However, for the purpose of detecting 
the presence of m^, this technique is not the optimal one: in fact eq. fll.59p contains 
the unknown Atg which is model dependent. A better strategy is to make use of 
gravitational radiation only, as suggested by Will [76]. In terms of the frequency /, the 
energy reads E = hf, where h is Planck's constant; eq. fll.58p can then be re-written 
as 



At = Ata - (1 + z)Ate 



(1.59) 




(1.60) 




(1.61) 



where we have introduced the graviton Compton wavelength and assumed that the 
frequency is such that hf » mgC^. Through eq. fll.6ip the velocity difference is 
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related to the graviton frequency, which lends itself to be tested by detecting GW from 
compact binaries. In fact, during the observation of an inspiral, the frequency chirps: 
it enters the sensitivity band of an interferometer at its lower end, for example ~ 50Hz 
in the case of ground based instruments where low-frequency seismic noise is relevant, 
and exits the band when the binary reaches the last stable orbit before merger, around 
~ 500Hz for first generation instruments. According to eq. fll.6ip . this frequency sweep 
will cause Vg to increase towards c, which in turn modifies At according to eq. ( ]1.60p . 
with a shrinking effect. From the point of view of GW interferometry, this shrink will 
reduce the amount of time that a given number of cycles spends in the frequency band 
of the instrument. In order to detect this effect by means of interferometers. Will 
computes how the time change is reflected in the Fourier transform of the GW phase. 
Within a 1.5 PN accuracy, the result is that a finite graviton Compton wavelength 
introduces an extra term at IPN order, i.e. it causes a shift 



128z/ vr^DM , , 

To study and constrain the effect of terms like eq. (11.621) . the strategy of refs. [151 [161 
[T7] is analogous to what is done for PK parameters: two ipj's are used to plot two curves 
of the type 7712(^1, ipj ^Aipj) in order to fix the masses of the emitting binary (within 
errors); afterwards, to test the theory, one draws (at least) a third curve for a test 
parameter ipT to see if it overlaps in a region of the mass-mass plane that is common 
with the first two parameters used. The difference with binary pulsars observation 
is that in the present case all the information about the emitting sources has to be 
extracted from the phase. In GW astronomy in fact, signal detection and parameter 
extraction are realized through matched-filtering, a technique where the output time 
series of the interferometer is convoluted with a bank of theoretical waveforms called 
templates. Every template depends on specific values of the ipj^^ so that, maximizing 
over the signal-to- noise ratio (SNR), one can identify the template and the ip/s that 
matched the best: this theoretical waveform can then be considered as representative 
of the actual signal (within a certain confidence level). 

A first possibility is to fit a GW signal through a template where all the ipj^s 
are treated as independent of one another and try to measure each of the PN coeffi- 
cients |T5]. This type of tests requires SNRs as high as 1000, which could be the case 
of a supermassive black-hole binaries observed by the space interferometer LISA for 
one year before merger. In order to study the feasibility of the tests, one can consider 
GW events in such a way as to minimize the errors in measuring the PN coefficients, 
which results from a compromise among many quantities: SNR, binary total mass and 
luminosity distance of the source from the experiment. Moreover, for the purpose of 
estimating phasing coefficients, one can model LISA as a single interferometer and neg- 
lect the orbital motion of its spacecrafts. These assumptions are justified by the fact 
that the angles describing the position of the source and the orientation of LISA are 
parameters that modify the amplitude of the signal and not its phase: for this reason, 
there is no correlation among the angles and the physical parameters contained in the 
phase. Within this context, ref. [15] obtains the errors reported in Fig. 11.14] where the 
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Figure 1.14: Errors in the measurement of the PN coefficients ipi's of the GW phasing 
formula f ll.SSp with LISA. The coefficients 4'5i and ipei of eq. fll.SSp are here absorbed in 
ip5 and ipGy respectively. LISA is modeled as a single interferometer with no amplitude 
modulation (see text for explanation). The sources are optimally oriented supermassive 
black- hole binaries located at a luminosity distance of = 3 Gpc. Figure taken from 
ref. m. 



coefficients V'sz and ipei are absorbed in -05 and V'e, respectively. Looking at the figure, 
one can see that the best candidate source for a test with LISA is a black hole binary 
with total mass of M = 2 x lO^'^M©: this is indeed the system chosen in ref. |15]. The 
resulting test is represented in Fig. 11.151 where the consistency among the parameters, 
comprising of their la error bands, is evident: the curves m2{mi, tpj^Atpj) have indeed 
a non- vanishing intersection in the point Mi = M2 = 1 x 1O^'^M0. 

With one-year observations of coalescences in the mass range (2 x 10^-2 x 10'')Mq, 
ref. [15] concluded that LISA tests could distinguish massive gravity theories from 
GR at the IPN level, provided the graviton Compton wavelength is not larger than 
(5.5 X lO^'' - 3.8 X 10^^)km, i.e. (2 x 10^ - 1 x 10^)pc. It is interesting to convert this bound 
in a constrain on the graviton mass. In units h = c = 1 one has 1 = ~ 200 MeVfm"^, 
where fm=fermi=10~^^m: with a conversion factor of lm=(2 x 10~^ev)~^, the highest 
bound on Xg gives rUg < 10~^^eV. 

As we have seen, a test where one tries to measure the PN coefficients all at the 
same time is only possible for a very specific choice of the constituents masses and for 
very high SNR. In reality, the parameters are correlated among them: this covariance 
enhances the errors in parameter estimation and dilutes the accessible information. 

A more effective test is the one proposed in ref. [16j , where the consistency check is 
run on one coefficient at a time. The coefficients tpo and 1/J2 play the role of a base with 
respect to which one expresses all the other parameters except the one on which the 
test is realized, which is taken as independent too. The choice of tpo and 1^2 as a base 
is motivated by the fact that, being the least suppressed in the PN expansion, these 
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Figure 1.15: Plot showing the regions in the M1-M2 plane that correspond to la 
uncertainties in the PN parameters ip^s. The source is a 2 x black hole binary 

at luminosity distance of = IGpc: it is imagined to be observed for a year by 
LISA as described by a single interferometer with no amplitude modulation (see text 
for explanation). Figure taken from ref. [77] . 



parameters are the best measured; moreover they are the only ones which will not be 
modified by spin when one decides to introduce this degree of freedom in the treatment 
too. Therefore, test parameter ipx is chosen in the set ipT = 4'3,4'4:,4'5h4'6,4'6hi^7 and 
the GW signal is fitted with a template which depends on ipo, ip2 and ipT- 

From Fig. 11.161 one can see that the choice of a restricted parameter set improves 
the accuracy in estimating the coefficients: all the regions spanned by the parameters 
in the mi-m2 plane are smaller than those of Fig. 11.151 most notably for the higher PN 
order coefficients. It is remarkable that, also in the present case, the coefficient ip4 is 
the worst determined: even the highest PN orders perform much better. This is mainly 
due to the structure of the phasing function (11.551) . where a PN coefficient of order k 
brings a correction oc f''^^ to the overall factor f~^^^. This is the reason why could 
be discarded from the very start of the study. Because of the phasing structure, terms 
which are close to ip^ in the expansion will be slowly varying with frequency. However, 
■ipe multiplies the term f~^^^^, while -ipi multiplies the term f~^^^, which is a decreasing 
function of the frequency: therefore, the contribution of 1^4 will be close to and then 
very much affected by the covariance with ip^,. 

The results presented so far refer to tests conducted with a possible GW signal 
detected by the space interferometer LISA. Ref. [16] considered also tests realized 
using ground-based interferometers, of second and third generation, like Advanced 
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Figure 1.16: Plot showing the regions in the mi-m2 plane that correspond to 
la uncertainties in a test parameter chosen one at a time form the set ipT = 
ip3,ip4,ip^i,iljQ,'iljQi,ipY. The source is a supermassive black hole binary with masses 
(10^, 10^)Mq, located at a redshift of z = 1; it is supposed that the source is observed 
for a year by a single interferometer of LISA type (see text for explanation). Note that 
the la uncertainty in ip^ is smaller than the thickness of the line. Comparison with 
Fig. 11.151 shows that using a single test parameter at a time enhances the estimation 
precision. Figure taken from ref . [77] . 



LIGO [78j and the Einstein Telescope (ET) |79j. The main conclusion is that already 
a second generation detector can provide tests which are sufficiently accurate to start 
probing GR through GWs. The existence of such a potential has been confirmed by 
a very recent study [17], which extended the investigation to the use of so-called full 
waveforms (FWF) as opposed to the restricted ones (RWF). Investigating the mass 
range (11 — 11O)M0 in case of Advanced LIGO, it is found that one could measure ip3 
with a fractional accuracy better than 6% and ip^i with a fractional accuracy better 
than 23%, corresponding to an improvement varying from 3 to 100 with respect to 
the use of RWF. Most importantly, ref. \17\ analyzes in detail the efficiency of the test 
realized with the reduced set of coefficients as compared to the full set. In the latter 
case, if, for example, the 1.5PN coefficient did not overlap with the regions allowed by 
the other parameters, one could directly ascribe this to a failure of GR at 1.5PN order. 
This correspondence could not be so evident when using the reduced set of coefficients. 
Let us stick with ip^ as the testing parameter and suppose, for example, that the correct 
theory is Chern-Simons gravity [17], which starts deviating from GR at 2PN order. The 
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template used for fitting the signal will have three independent parameters: ipo, ijj2 and 
ips, meaning that ipo and 1^2 are parametrizing all the coefficients of order higher than 
k = 2, except for tp^. Employing such a parametrization in presence of the deviation 
implies that one is describing non-GR effects by means of GR terms: this can affect 
the estimation of the test parameter ip^, even if the 1.5 PN order is not the one from 
which the deviation actually starts. With these caveats, it is possible to look at the 
consistency checks of Fig. 11.171 which refers to the third generation interferometer ET: 
from the various plots one can see that each test parameter investigated consistently 
overlaps with the basis pair ('0o,'?/'2)- 

Before ending this section, it should be mentioned that very recently, another phe- 
nomenological framework has been proposed to discuss deviations of the GW phasing 
formula from the GR prediction. In ref. [80], Yunes and Pretorius have put forward 
what they call the parametrized post-Einstein (PPE) approach, where they parallel the 
PPN spirit: while in the latter one parametrizes the metric in a model- indipendent way, 
in the former one parametrizes the full waveform, i.e. both the phase and the amp- 
litude. This is achieved using the existing knowledge about various alternative theories 
of gravity such as Brans-Dicke [26], massive graviton effects [76] and Chern-Simons 
theory [81] to write down a generic parametrized waveform in the Fourier domain. 
Apart from extending the use of parameters to the amplitude too, another difference 
with refs. [151 UHl [H] is that the PPE also considered the contributions from the merger 
and ringdown phases of the binary evolution beyond the inspiral. In both approaches 
spin degrees of freedom have not been included for the moment. 

It is in this active domain of research that the parametrized framework presented in 
this thesis sets itself; in Chapters H] and [5] I will discuss a field-theoretical extension of 
the PPN that stems from the EFT approach to inspiraling binaries: as we will see, this 
is another convenient attitude to build parametrized waveforms for GW astronomy. 
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Figure 1.17: Plots showing the regions in the mi-m2 plane that correspond to la 
uncertainties in the basis parameters ipQ, ip2 and in the test parameter chosen one at a 
time form the set ipx = ip3,ip4,ip5i,ijjQ,ipQiy'^7- The source is a black- hole binary with 
masses (2,2O)M0, located at a luminosity distance of Di = 300 Mpc; it is supposed 
that the source is observed by the Einstein Telescope (ET), a third-generation ground 
based interferometer. The basis parameters ipo and iIj2 are those from which one can 
measure the masses of the two black holes. Each test parameter, together with its 
error, spans a given region in the mi-m2-plane: if GR is the correct theory of gravity, 
then the three parameters ipo, "4^2 and ipT should have a non-empty intersection in the 
mi-m2 plane. A smaller region leads to a stronger test. Notice that all panels have 
the same scaling except the top middle panel in which Y axis has been scaled by a 
factor of 10. With respect to Fig. I1.16[ the waveforms used for the tests presented here 
comprise amplitude corrections. Figure taken from ref. [17j. 



Chapter 2 



The effective field theory approach 
to inspiraling compact binaries 

2.1 An invitation 

During the inspiral phase of a binary system the constituents are separated enough 
for their dynamics to be non-relativistic, i.e. the leading order interaction between 
the two bodies is Newtonian gravity and corrections to it can be organized in the PN 
approximation of GR. 

Being perturbative in nature, the PN expansion can also be interpreted from a 
field theory standpoint and reformulated in terms of Feynman diagrams. As I have 
discussed in Sec. II. ![ a first application of these instruments to the two-body problem 
in classical gravity dates back to 1960 p^; more recently ref. [82] by DEF has started a 
series of works in which the PN limit of scalar-tensor theories of gravity is conveniently 
studied by means of diagrammatics. In the last few years the use of these instruments 
in classical gravity has revived thanks to the introduction of an EFT approach [6]. In 
this context the non-relativistic dynamics of the two-body system is studied in a way 
analogous to bound states in QED and QCD for such a reason the EFT framework 
of ref. [6] has been dubbed non relativistic General Relativity (NRGR). 

The use of an EFT framework in the two-body problem of GR is motivated by the 
fact that the binary dynamics is a typical situation in which the disparity of physical 
scales allows one to separate the relevant degrees of freedom in subsets with decoupled 
dynamics. For example, because the extension of the radiating source is much smaller 
than the wavelength of gravitational perturbations, the modes corresponding to the in- 
ternal dynamics of the compact object can be integrated out and the source is effectively 
described as a point particle. 

At high enough order the point particle approximation breaks down and causes 
divergences to appear. However this happens in any field theory coupled to point 
sources and only means that the theory needs to be supplemented by a more complete 
(fundamental) model. Quoting ref. [6], the NRGR attitude towards the use of the 
point particle description and the related problem of divergences can be summarized 
as follows: "Rather than trying to resolve the point particle singularities by using a 
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match 




Figure 2.1: Graphical representation of the step-hke construction of a tower of EFT 
theories for the two-body problem in GR. From top to bottom, the characteristic values 
of the length scale i are: the size of the compact object Rs, the orbital radius r and 
the radiation wavelength A. Figure courtesy of Chad Galley. 



specific model of the short distance physics, in an EFT framework we systematically 
parametrize our ignorance of this structure by including in the effective point particle 
Lagrangian the most general set of operators consistent with the symmetries of the well 
understood, long wavelength physics. In our case, this long wavelength physics is GR 
and the symmetry that constrains the dynamics is just general coordinate invariance.". 
This parametrization is a model-independent description of the compact source in terms 
of a point particle. In fact, the EFT is "tuned" to reproduce the low energy regime 
of the full theory by means of the so-called matching, a procedure which consists of 
two stages. In the first, one calculates a convenient observable in the full theory and 
expands the result in the low energy limit. This form of the result is then compared with 
the prediction of the EFT which depends on the coefficients of the effective operators: 
by equating the two results one can then fix these coefficients. 

Making use of the action for gravity as dictated by GR and the sum of two point 
particle actions is an adequate description of the two-body system as far as the objects 
are isolated, i.e. until the focus is on length scales smaller than the orbital radius. 
At this scale another EFT has to be built to describe the coupling of radiation with 
a composite, two-body source. This step-like approach is conveniently represented by 
Fig. 12.11 which graphically describes the construction of a tower of EFT theories. It 
should be stressed that these theories are built one at a time and one for each of the 
characteristic values of the length scale i, which, with reference to the figure, are: the 
size of the compact object Rs, the orbital radius r and the radiation wavelength A. It 
is this separation of the intricacies of the binary problem that makes NRGR an ideal 
framework to study gravitational radiation: NRGR is effective in the sense that it is op- 
timized to streamline the perturbative calculations of the PN regime. Moreover, NRGR 
makes a systematic use of Feynman diagrams, which can be dealt with by automatic 
calculation on a computer. By means of these vantage points, it has been possible to 
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obtain previously un-computed spin contributions to the source multipoles at next-to- 
leading order this provides the last missing ingredient required to determine the 
phase evolution to 3PN accuracy including all spin effects. 

Besides being a convenient tool for calculations, Feynman diagrams can be used 
to recast the non-linearities of GR in a field-theoretical language: by means of such a 
reformulation, it is possible to have a different point of view on the physical information 
that can be extracted from experiments of relativistic gravity. This is the standpoint 
which is taken in the present thesis and it is exposed in Chapters H] and (5) there I 
show how an EFT approach allows one to investigate the extent to which experiments 
can constrain the values of the non-Abelian vertices of GR. This represents a specific 
field-theoretical interpretation of observations but has a general relevance because the 
analysis is model-indipendent. In this chapter I set the stage for these investigations 
by giving an introduction to NRGR based on refs. [SI [83]; along the treatment, I find it 
convenient to reproduce some formulae of ref. |83] as they stand, notably those where 
Feynman diagrams are drawn inside equations thanks to the feyn package for ETf^X. 
With respect to the field-theoretical investigations of DEF that I have sketched in 
Sec. II. H here I content myself with GR and I do not take extra fields into account: 
in such a way, I can keep the discussion self-contained and describe the original spirit 
of the approach. An example of the application of the EFT to the case of one extra 
scalar field is the subject of Chap. |3], where I report a study of the renormalization of 
the energy-momentum tensor for point-like and string-like sources. 



2.2 Effective field theories 

The EFT treatment of binary systems pioneered in ref. [6] is motivated by the 
presence of three length scales in the problem (see Fig. \2.2\i : 

Rs = size of compact objects, 
r = orbital radius, 

A = wavelength of emitted radiation (2.1) 

which turn out to be interconnected by means of the relative velocity v. For a gravit- 
ationally bound system, in fact, one can use the virial theorem to obtain the estimate 

= ^-^' (2-2) 

where we have introduced the Schwarzschild radius Rs = 20^771; moreover, because 
the binary's motion is circular with frequency Uq, one can write 

r 

V = Uor - (2.3) 
A 

in which use has been made of the fact that the GW frequency uqw is proportional to 
the orbital one, ojgw — <^o- 
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^ ^ - - Rs 



Figure 2.2: A schematic description of an inspirahng compact binary where the relevant 
length scales are shown. The disparity of these scales makes EFT techniques very well 
suited to describe binary systems. Figure courtesy of Chad Galley. 



Because of the link (12. 2p between the velocity of the sources and their gravitatio- 
nal-field strengths, the length scales in (12. ip are indeed connected one another through 
the PN expansion parameter v. this shows that high PN orders encode physics from 
different scales at the same time. Through this link it is also evident that the scales 
respect the following hierarchy: Rs << r << A. This type of problems can be conveni- 
ently dealt with by means of an EFT where the physics which is relevant at a definite 
scale is treated separately and "once for all". 

In particle physics, EFT schemes are typically used to study the impact of some 
ultra-violet (UV) scale A on the physics at low energy characterized by the infrared (IR) 
scale 9, with 9 « A. For the sake of concreteness, let us consider a theory described 
by the action functional S[(j), $], with collectively denoting light (or massless) degrees 
of freedom and $ standing for some heavy fields with masses M$ ~ A . If the focus 
of the investigation is an energy regime E 9, it is not necessary to resort to the full 
theory iS[0, $]: because the fields $ are too heavy to be excited at energies E ^ 9, only 
the fields will have a non-trivial dynamics; therefore, it is convenient to solve for the 
fields $ and plug the solution back in the action (cfr. the Fokker action in Sec. II. ip . 
Formally, this is achieved by integrating out the heavy modes through the path integral 



where the (local) operators Oi{x) are, in principle, infinitely many and the coefficients 
of the expansion Cj, called Wilson coefficients, are energy-dependent. This effective 
action does not represent a "new" theory but rather approximates the full one in some 
regime; as a consequence the operators Oi{x) must respect the original symmetries. 
For practical purposes, only a finite number of such operators contribute to eq. (12. 5p . 
The effective action is supposed to reproduce the predictions of the full theory for the 
low energy dynamics up to some corrections, i.e. eq. (12. 5p is a perturbative series in 
a small parameter e = 9 /A. This is related to the fact that, because every operator 
Oi has some mass dimensions Aj, the value of its Wilson coefficient at an energy scale 




(2.4) 



and results in an effective action of the type 




(2.5) 
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H = A can be shown to be0 



= A) 



AA.-4 



(2.6) 



where ~ therefore, the higher the mass dimension Aj is with respect to 

4, the lower the contribution of the corresponding operator is to the effective action; 
operators with Aj < 4 have their coefficients renormalized, i.e. redefined in order to 
take the effect of the cutoff into account. As a consequence, the impact of the high- 
energy physics on the IR dynamics is very simple: all the UV dependence appears 
directly in the coefficients of the effective action, which makes it possible to determine 
how a low energy observable depends on A. More specifically, this decoupling between 
the two scales allows one to determine which operators are needed at a given order in 
the perturbative expansion by using a generalized version of dimensional analysis called 
power counting. Let us see how this works in practice. Because accounts for the low 
energy dynamics of the theory, its four-momentum is of the order of the IR scale, i.e. 

9: the same scaling will characterize derivatives 9^ ~ ^, while for coordinates one 
has ~ 9^^. With these rules, one can power-count every term in the action (12. 5p . 
notably the scalings of different operators Oi with respect to each other. The leading 
order term in the perturbative series (12. 5p is the kinetic term, which scales as 



where the overall mass dimension is zero because the action is a scalar under Lorentz 
transformations. Eq. (12. 7p tells us that the scaling of the field itself is ~ 6'. Therefore, 
an interaction term like c^{dfj_(pd^(l)Y will be power-counted as {9 / A)^ and, when integ- 
rated to form an action, it will give a contribution which is suppressed with respect to 
the leading order term (12. 7p by [9/ A)^. Power-counting terms in the action enables one 
to determine which operators are needed to compute physical observables at a given 
perturbative order. 

Building an effective action is convenient in two typical situations, depending on 
whether the full theory is known or not. An example of the first case is represented 
by the electroweak interactions: if one is interested in energies E -C Tnw,Zi it makes 
sense to integrate out the VF, Z bosons because they are too heavy to have a dynamics. 
The result is Fermi theory of weak interactions, which makes it simpler to calculate 
amplitudes like the one depicted in Fig. 12. 3[ where the weak process receives a QCD 
correction [H]. Even when the full theory is not completely known the decoupling 
allows one to treat the low energy dynamics. The structure of the effective action is 
dictated by the symmetries that survive at low energies: by writing down the most 
general set of operators consistent with these symmetries, one is accounting for the 
UV physics in a completely model independent way. For instance, the Standard Model 
of particle physics is believed to be an EFT below a scale A ~ 1 TeV. The indirect 
manifestation of this scale has already been constrained with the LEP experiment at 




(2.7) 



1. Remember we are using units in which c = ^ = 1 so the scahngs are naturaUy expressed in terms 
of mass dimensions. 
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g 




(a) (b) 



Figure 2.3: QCD correction to a typical four-Fermi process in the Standard Model, 
calculated in (a) the full theory with propagating W bosons and (b) the effective 
Fermi theory of weak interactions. Graph (b) reproduces graph (a) up to corrections 
suppressed by powers of E'^/m'^ ^ 1. Figure taken from ref. [83] . 



CERN and with the SLD experiment at SLAG [85], while the direct probe of the TeV 
scale is awaited from the LHG at GERN. 

A final case where an EFT approach is useful is represented by a theory that is 
known but strongly interacting. The example that mainly concerns the present thesis 
is GR: this can be regarded as an effective theory of quantum gravity with predictive 
power below the strong coupling scale Mpi ~ 10^^ GeV. The EFT interpretation of 
GR is reviewed in detail in refs. [86l [871 IHH] and allows to address some issues that 
emerge within a quantum field theory reformulation of Einstein's theory. For example, 
in an EFT framework, the fact that GR is non-renormalizable does not constitute a 
problem: the UV divergences that start arising at first quantum order are dealt with by 
renormalizing the parameters of higher derivative terms in the action; perturbatively, 
only a finite number of terms is required for each order. The specific EFT of ref. [6] 
is designed to study the classical dynamics of binary systems so that, as we will see, 
NRGR does not need to care about quantum corrections. From the following sections, 
it will be clear that even in a classical context, it is useful to resort to quantum field 
theory tools like Feynman diagrams: in analogy with Fig. 12. 3[ the NRGR description 
of non-relativistic compact binaries can be schematically described as in Fig. 12.41 

The considerations spelled out so far are not exhaustive of EFT features but they 
suffice to move to the discussion of the two-body problem: we will come back to the 
EFT ingredients directly when they are at play. 



2.3 Non-Relativistic General Relativity 

In the dynamics of a binary system there are three separated length scales (12. ip . 
so the EFT treatment will make use of the path integral (12. 4p in two stages: the first 
for integrating out the physics pertaining the internal structure of the compact object, 
the second for that of the orbital dynamics. This is what we will see explicitly in the 
following, starting with the scale represented by the Schwarzschild radius Rs- 
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(a) 



(b) 



Figure 2.4: In analogy with Fig. 12. 3[ this picture is a very schematic representation 
of how non-relativistic compact binaries are treated in an effective field theory, where 
the perturbative PN expansion is expressed in terms of Feynman diagrams. Panel (a) 
represents the orbital motion of a binary pulsar and the gravitational radiation emitted 
by the system, panel (b) is the Feynman diagram which accounts for the lowest order 
interaction between the binary constituents, i.e. the Newtonian potential. 

2.3.1 The effective theory for isolated compact objects 

The first length scale that has to be integrated out is the size of the compact object. 
According to the EFT philosophy, one can build an effective action by identifying the 
degrees of freedom which are relevant at the scale of interest: including the most general 
operators consistent with the symmetries provides a model-independent treatment of 
a system. To describe GW observables, the relevant degrees of freedom are the long 
wavelength gravitational perturbations and the world-lines of effective point particles. 
In a first approximation one can neglect the spins of the objects: as we saw in eq. (11.24p . 
spin effects in the equations of motion do not show up until 2PN order. 

The symmetries that govern the dynamics of these degrees of freedom are 

1. general coordinate invariance, — )■ x'^(x); 

2. worldline reparametrization invariance (RPI), c — )■ o'ic); 

3. invariance under 5*0(3) transformations, which guarantees that the compact ob- 
ject is perfectly spherical and that it has no permanent moments relative to its 
own rest frame. 

Such EFT is the appropriate one for describing non-rotating (Schwarzschild) BHs in- 
teracting with external gravitational fields. The effective action consistent with these 
criteria is then 



Seff[x'',g^,u] = SEH[g] + Spp[x,g], 
where the action for gravity is the usual Einstein-Hilbert action|fl 





(2.9) 



with Mp^ = 327rGiv and R{x) the Ricci scalar, and the point particle term is the usual 
minimal coupling 




2. With respect to Sec. 11.11 here we adopt the "mostly minus" convention for the metric 
77py=diag(+,--,-). 
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where m is the particle mass. Both actions can be supplemented with terms con- 
taining higher order powers of the curvature. For the gravitational action (12. 9 p these 
terms are suppressed by powers of Mpi, the scale where quantum gravity is expected: 
these contributions are completely negligible in the classical regime of binary dynam- 
ics. Concerning the point particle action (I2.10p the situation is different. Higher order 
operators that contain the Ricci tensor R^^ turn out to be redundant because they can 
be removed by a field redefinition (see the appendix of ref. [83] for details). A second 
possibility is to have higher powers of the Riemann tensor R/^uap, notably through its 
components of electric and magnetic type parity Efj^i, and B^i,, respectively, which are 
given by 

E^u = R^au^x'^x^, (2.11) 

Therefore, one can re-write the point particle action ( ]2.10p as 

Spp = -m j dr + CE j drE^^E^"" + cb J dTB^.B^"" + ■ ■ ■ : (2.12) 

The operators constructed from E^^^, and Bf^^, are the first in an infinite series of terms 
that systematically encode finite-size effects, which are due to the physical extension 
of the source, in our case a black hole. One way to see this is to compute the matching 
coefficients ce,b- the explicit matching calculation [6l [83] shows that both ce and 
Cb scale as Mp^Rg, so they approach zero quite fast for a vanishing Rs- Therefore, 
the operators built from and B^^, do not contribute to the binary dynamics of 
Schwarzschild black holes until 5PN: this gives a proof of the effacement principle of GR 
in the EFT context. Another way of seeing how the non-minimal terms are associated 
to finite-size effects is to look at the equations of motion that derive from eq. (I2.12p : 
ahead of the calculation, one can see that the particle will not move on a geodesic of 
the background gravitational field: in fact, this is what one gets by extremizing the 
minimal action (I2.10p . Indeed, geodesic deviation implies stretching by tidal forces, 
which occurs when one considers the motion of extended objects in a gravitational 
field. As a consequence, one can rely on the minimal point particle action fl2.10p up 
to the very high 5PN order for what concerns the dynamics of two non-rotating black 
holes. Schematically, this is described by Fig. 12.51 which points at the first EFT to be 
constructed in the tower of theories of Fig. 12.11 

To study the non-relativistic dynamics of a binary system, one can expand both 
the gravitational action (12. 9p and the point particle action (I2.10p by defining the gra- 
vitational field as a small perturbation over flat space-time in the following wayEI 

Mpi 

The perturbative expansion can be very conveniently dealt with by means of Feynman 
diagrams, i.e. by associating drawing signs to the various fields and interaction vertices. 



3. This choice is not the only one possible, cfr. eq. p.lOp of Chap. [3l 
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Figure 2.5: A schematic representation of what the EFT for isolated compact objects 
is about: the extended astrophysical object can be effectively regarded as a point- 
particle probe of the long-wavelength gravitational perturbations. Figure courtesy of 
Chad Galley. 



We have already met an example of a symbolic representation dictionary in Sec. II. ![ 
where the rules for diagrammatic computations in the context of scalar-tensor theories 
have been reported in Fig. 11.21 In the case of Einstein gravity, the Feynman rules for 
perturbative calculations can be found in the classic lectures by Veltman [89] and the 
more recent review by Donoghue [86]. It is custom to assign curly or wavy lines to the 
gauge bosons of a theory, like photons and gravitons, and continuous lines to fermions 
in the matter sector. For NRGR one can stick to these conventions. However, as we 
will see, the use of continuous lines for matter sources should not be taken as meaning 
that these are quantum degrees of freedom. 

Using the expansion f l2.13p for the purely gravitational sector (12. 9p one has 



2M|i 



{dhf + 



hjdhf h\dhf 
Mp, ^ 



+ ■ ■ 



PI 



-1 




+ ■ ■ 



(2.14) 



where gravitational self-interactions are represented by Feynman vertices containing 
many graviton lines. The first term in the expansion (12.141) represents free propagation 
of a graviton from a space-time point to another; in a compact fashion, this can be 
written as 

{h^y{x)Ki3{y)) = Df{x - y)P^y-ap , (2.15) 
where the Feynman propagator Df{x — y) is given by 



Df{x - y) 



{27tY F + ie 



-ik-{x—y) 



(2.16) 



and the tensorial structure reads 



(2.17) 



For what concerns the interaction of the gravitational field with the point particle. 
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plugging the definition (12.13P in eq. f l2.10p gives 

1/2 



m 



Mpi 



M, 



PI 



m dr 



m 



2Mi 



PI 



m 



'PI 



8M^ / driJi^yX^x ) + ■ ■ 



(2.18) 



where df'^ = rj^^dx^dx^ and x^ = dx^/df. In order to draw the Feynman diagrams 
corresponding to these vertices it is important to remark that the matter sources to 
which gravity couples are macroscopic objects. As a consequence, the binary constitu- 
ents are non-relativist ic particles endowed with typical three-momentum of the order 
p ~ mw, with V the orbital velocity; on the other hand, gravitons will rather have 
three-momentum k ~ ^, where h is the Planck constant and we use c = 1. When a 
compact object emits a single graviton, momentum is effectively not conserved and the 
non-relativistic particle recoils of a fractional amount roughly given by 



M 

IpI 



|k| 
IpI 



h 

L 



where L ~ mvr is the angular momentum of the system: it is clear that for macroscopic 
systems such quantity is utterly negligible. To summarize, in NRGR one describes the 
massive constituents of binary systems as static, background sources of gravitons. In 
this context, the Feynman diagrams corresponding to the vertices of eq. fl2.18p are given 
by 



L 1 




+ 






i 1 



+ 



(2.19) 



where: 

the continuous lines represent matter sources and they are straight to indicate that 
there is no recoil in the emission of gravitons; 

the arrows on the continuous lines describe the flow of time of the particle world-line; 

each curly line represents a factor h^j^x'^x" (see eq. fl2.18p ). 

Composing these matter-gravity vertices with those of the purely gravitational sector 
(12.141) . one ends up with diagrams of the type of Fig. 12.61 that describe the interaction 
between two sources. By calculating the corresponding amplitudes, one is integrating 





Figure 2.6: The flrst few diagrams contributing to Seff{xa) in Lorentz covariant per- 
turbation theory. Figure taken from ref. |83] . 
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out the gravitational perturbations h^i, of eq. f l2.13p . The result is an effective action 
where the only degrees of freedom which are present are the particles' world-lines: this 
corresponds to the functional Seff{xa) defined via a path integral of the type of eq. (12. 4p 

exp[iSeff{xa)] = J Dh^^{x)eyip[iSEHih) + iSpp{h,Xa)], (2.20) 

where the particle world-lines a;^(r) that source the gravitational perturbations h^^, are 
held fixed. Once the effective action ( 12.20p has been obtained, it is possible to calculate 
physical observables. Notably, taking the real part of eq. ( 12.20P generates the coupled 
equations of motion for the two-body system and consequently the mechanical binding 
energy. On the other hand, the imaginary part of the effective action measures the 
total number of gravitons emitted by a fixed two-particle configuration over an 

arbitarily large time T — )■ oo 

ilmS.„W = i/dBdn^. (2.21) 

where dT is the differential rate for graviton emission from the binary system. Although 
the graviton number is not a well defined observable classically, the classical power 
spectrum can be obtained by integrating dP = EdV and reads 

As discussed in Sec. 11.41 once we have E and P, we can derive the GW phase v^gw 
seen by the detector; for example, in the mass ratio limit the phase reads 

^VGw{t) = f drcocwir) = T dv'v''^§l^ , (2.23) 

J to (jnM i^[y ) 

where ojgw is the GW frequency, M the mass of the heavier body in the system and 
where, with respect to eq. fll.52p . we wrote the GW flux as P instead of F. 

An example of amplitude computation 

To see how the computation of scattering processes like those of Fig. 12.61 proceeds, 
here I will take as an example the case of the one-graviton exchange diagram in 
panel (a). The gravity-matter interaction is given by the linear term in eq. (I2.18P : 
for the a-th source this reads 



2M] 

where I have specified the space-time dependence of the field /i^jy and where I remind 
that the four- velocity of the particle satisfies x'^x'^ r]^^ = 1 . Including an analogous 
term for particle b gives a factor of hai3{xb)x'j^x^ and therefore the propagator of the 
gravitational field f l2.15p : 

)Pf.u;aP . (2.24) 
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The tensorial part of eq. fl2.24p brings the following contraction between the four- 
velocities of the particles: 

^a^a P^iJ^;al3 ^b^b ~ l^^a^a [Va^iVlSu + 'nau'np^i ~ "n^ivTlap] ^b^b (2.25) 
= ^ [{ia ■ ib) + (ia ■ Xb) - {Xa " Xa){Xb " Xfe)] 

= ^ [2{Xa ■ XbY - 1] 

Combining all the terms, the amplitude for the one-graviton exchange diagram in 
Fig. 12.6( a) gives the following contribution to the effective action 




(2.26) 



where the numerical factor of 1/8 comes from the two matter-gravity vertices together 
with the propagator, while the sign difference with respect to eq. f l2.25p is due to the 
fact that every Feynman rule is defined with a factor of i. As it stands, eq. f l2.26p 
cannot be assigned a unique power of v to assess the perturbative order at which it 
contributes to the effective action and, then, to GW observables. This is because the 
perturbative series of Fig. 12.61 is fully covariant and, as such, is not optimal to perform 
calculations in the limit of small three-velocity. As we will see below, for a term like 
the one-graviton exchange it is easy to do the expansion for small velocity; on the other 
hand, in a diagram like the one in Fig. 12.6( c). it would be more cumbersome to keep 
track of all the necessary terms at a given order in v. Since in the end one is interested 
in computing GW observables up to a fixed order in the velocity, it is crucial to develop 
a set of rules that assign a unique power of v to each diagram in the theory. Such a 
set of rules constitutes the power counting, a characteristic ingredient of every EFT as 
we discussed in Sec. 12.21 For NRGR this is presented in the following sections. 

2.3.2 Power counting and mode decomposition for gravita- 
tional perturbations 

To establish a power counting scheme, it is very useful to recognize that one can 
decompose the gravitational perturbations as follows: on one hand there is a back- 
ground whose range of variability is the wavelength of gravitational radiation A; on 
the other, there is a set of perturbations that have support only on scales below the 
orbital radius r. In other words, one can distinguish between hard modes and soft 
modes. A similar decomposition of gravitational perturbations is adopted to distin- 
guish GWs from a background field and to define the energy-momentum tensor carried 
by GWs [90] 0. Because gravitational radiation is characterized by the scale A which 

4. A complete treatment of the so-called short-wave approximation can be found in Sees. 1.4 and 
1.5 of ref. \90\, while the "Further reading" at the end of Chap. 1 contains the relevant citations. 
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is large compared to r, one can take A as the IR scale of the full theory and r as the 
UV one. This is reminiscent of the general discussion on EFTs that I presented in 
Sec. 12.21 Because A and r are two separate scales, the physics characterizing each one 
is decoupled from the other: in these situations an EFT approach is very convenient. 
On these grounds, one can think of GW observables as regarding the IR sector of the 
theory, while the gravitational dynamics below the orbital radius is a short distance 
phenomenon and corresponds to the UV region. The high energy modes are responsible 
for binding the binary system through forces of Newtonian type: for this reason, the 
UV modes are referred to as potential modes. Because retardation effects are negligible 
on the orbital scale, Newtonian dynamics is instantaneous: from a field-theoretical per- 
spective, potential modes are then off-shell and can only correspond to internal lines in 
a Fenyman diagram. The soft modes which have support on the IR scale A correspond 
to GWs and are called radiation modes: according to the particle physics parlance of 
EFT, physical radiation is described by on-s/ie// modes that are represented by external 
lines in a Feynman diagram. 

The separation of scales made explicit by the EFT framework is physically equi- 
valent to the matched asymptotic expansions adopted in the standard PN approach, 
which have been introduced in refs. pH |92]. Indeed, the consistency of the scale fac- 
torization proposed in NRGR has been recently verified up to order 3PN in the GW 
flux [93]. 

To make contact with the standard PN treatment, one can observe that poten- 
tial modes are those obeying a Poisson equation, so they correspond to quasi-static 
gravitational fields whose profile is set by the matter sources. On the other hand, radi- 
ation modes are the gravitational perturbations that satisfy a wave equation, i.e. they 
represent GWs. 

In order to illustrate more vividly both the meaning and the consequences of the 
mode decomposition, now I will go back to the one-graviton exchange amplitude fl2.26p 
and use it to compute the effective action in the non-relativistic limit: with this in hand, 
I will calculate gravitational observables to lowest order along the ways discussed after 
eq. ( I2.20p . Let us start the evaluation of eq. fl2.26p from the following term 



this integral can be split up into contributions from the following two regions of mo- 
menta: 



Space-like momenta of the form 

F->(fc°~^,k~V (2.28) 

In units c = 1, the typical three-velocity v is much smaller than one, therefore 
/c^ = k^k^ 7^ and these modes can never be on-shell; as a consequence, they 
never contribute to the emitted power through the imaginary part of the effective 
action: indeed, eq. fl2.27p can have an imaginary when the ie term in the propag- 
ator becomes important, i.e. when k'^ = 0. On the other hand, these modes will 
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determine the equations of motion through the real part of the effective action. 
We recognize these modes as being the potential gravitons. 

• Null momenta whose scaling is given by 

F^(A;°~^,k~^). (2.29) 

These modes satisfy k"^ ~ 0, so they can give rise to the imaginary part of 
the effective action and, therefore, to the radiation that propagates out to the 
detector. It is straightforward to identify these modes as the radiation modes. 

Since only potential modes contribute to the real part of the effective action, we may 
calculate the non-relativistic limit of this quantity by expanding the propagator in 
eq. (12:271) as 



1 1 



1 + ^ + 



~ [1 + Oiv')] ; (2.30) 



the leading order term k ^ gives the 1/r dependence of Newtonian potential once 
Fourier-transformed according to 



Ak-x 1 1 



■5(a;°) . (2.31) 



(27r)4 k2 47r|x 

Let me stress that we are computing the effective action with observables in mind: 
on one hand we will take the real part of the action to obtain the energy through the 
equations of motion, on the other we will take the imaginary part as an intermediate 
step towards the power emitted in GWs. To compute the equations of motion we 
have to consider another process with respect to the one-graviton exchange. This is 
because of the virial theorem, by which terms 0{G]^) contribute at the same order 
as those 0{v'^). As a consequence, the Newtonian contribution to the one-graviton 
exchange fl2.26p coming from eq. f l2.3ip is of the same order as the kinetic energy 
of the particles, which is represented by the diagrams in Fig. 12.71 To evaluate these 
diagrams, one has to expand the point particle action in the velocity, rather than in the 
gravitational field as we did in eq. fl2.18p : to first non trivial order in f , this expansion 
gives 

dTa = [r],,dxV]'/' = [{dxy - v^] ^ dx' \l - ^v^j . (2.32) 



With the expansions of the propagator and of the proper time, we can write the 
non-relativistic limit of the real part of the action as 

^eSefAxa] = \Y.j d^''^^^ - E / dx'^f^ + ■■■, (2.33) 

a a,b 

where we have used Gn = (327rM|i)^^. As it could have been awaited, this is just the 
action for classical non-relativistic particles interacting through a Newtonian potential. 
Note that the second term contains divergent self-energy contributions whenever a = b 
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0{v 



ib) 



Figure 2.7: Feynman diagrams describing free propagation of matter sources, i.e. the 
kinetic energy of the particles at 0{v'^). Because of the virial theorem, these diagrams 
are of the same order as the Newtonian potential, which is 0{Gn)- 



in the sum. Diagrammatically, they are described by a potential graviton which is 
emitted and reabsorbed by the same source. These divergences can be absorbed by 
renormalizing the masses of the particles. Formally, they can just be set to zero, by 
evaluating the momentum integral in dimensional regularization using the formula 



1 T{D/2-a) /x 



(27r 



(k2)" (47r)«/2 r(a) 



2 \ a-D/2 



(2.34) 



and taking the limit x — before setting D = 3. 

To calculate the imaginary part of the effective action one can use 



Im 



fc2 + ie 



-ITX- 



2|k 



■Kk^ = |k|) 



(2.35) 



which guarantees that only on-shell particles contribute to the radiated power. Then, 
the imaginary part of eq. fl2.26p gives 



1 



IGM^J (27r)3 2|k| 



(2.36) 



fcO = |k| 



This results in a differential power given by 



dP 



1 Gn^2 



dnd\k\ TAtt^ 



(2.37) 



A:0=|k| 



where 



so that k^ = /cq = and one obtains that the system does not source radiation through 
its mass monopole. The expression (12.371) coincides with what one would find both 
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from the classical calculation of the energy-momentum tensor at an asymptotically 
large distance from the source and from the tree-level amplitude for single graviton 
emission by the binary source. The latter alternative will be used later on in Sec. 12.3.41 

Coming back to the mode decomposition of the gravitational field, we split it up 
in potential modes H^y{x) propagating on the slowly- varying background field h^^{x), 
which represents the long-wavelength radiation: 

Ku{x) = H^,{x) + h^,{x). (2.38) 

For what concerns Feynman diagrams, we will draw the radiation modes hf^u with wavy 
or curly lines and the potential modes Hfj^^ with dashed lines. As a consequence of the 
scaling of their momenta eqs. f l2.28p and f l2.29p we have the following power counting 
for the two sets of modes 

doH^u ~ -H,, , d,H^, ~ -H^, (2.39) 



doh^i, ~ -hf,i, , dih^i, ~ -h^^ . (2.40) 

In order to treat all derivatives acting on fields on the same footing and assign them 
a common scaling ~ v/r, it is convenient to re- write H^j^i, in terms of its Fourier 
transform 

^e^''-^i^k^,(xo). (2.41) 

With this redefiniton, spatial derivatives on iJ^^, are replaced by factors of k multiplying 
Hkf,u{xo) and the only derivative left is doH]^f,u{xo) ~ {v/r)H]^f,^{xo). 

Having split the gravitational perturbations h into two sets of modes, the path 
integral of eq. fl2.20p now needs to be worked out in two stages: one over H, another 
over h. This would look like one is doubling the calculations to be done; rather on 
the contrary, this is one of the advantages of the EFT approach, as we defined it in 
the introduction: the EFT "allows one to separate the relevant degrees of freedom in 
subsets with decoupled dynamics". Dividing the calculation of the path integral fl2.20p 
in two stages is the aim of the EFT approach to binary systems: in fact, this choice 
eases the treatment of the intricacies which are inherent to the problem. 

In a first stage, one performs the path integral over the hard modes H]^^^ 

^iSNR(h,Xa) _ J £)JJ^^^^^O^^iSEH{h+H)+Spp{h+H,Xa) ^2 42) 

where h^i, is treated as a non-dynamical background. This is a formal way of doing 
the matching to the long-distance EFT containing hni, and the particle world-lines: in 
practice, it will be obtained through the calculation of Feynman diagrams that contain 
potential modes as internal lines. This action, which is derived in eq. (12.841) . constitutes 
the theory valid at the orbital radius: because of the matching, it contains explicitly 
the short distance scale r in the coefficients of operators. 
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Figure 2.8: The self-energy diagram whose imaginary part gives the radiated power 
(see text for discussion). 



In order to compute observables, the EFT prescription requires obtaining an effect- 
ive action that depends only on the particle world-lines. Therefore, given the action 
SNRih, x), the second stage in the NRGR approach is to perform the path integral over 
the radiation modes h,,i, 



Rather than solving this path integral, to obtain the action Seff{xa) one can resort to 
computing the diagrams that have no external graviton lines like the one in Fig. 12.81 
where the double line describes the binary system as a composite source; the orbital 
scale is un-resolved because one has integrated it out by solving for the potential modes. 
Feynman diagrams are easier to calculate in this theory than in the full theory, i.e. the 
one in which gravitational perturbations are not split into potential and radiation 
modes: in fact, momentum integrals are characterized by a single scale v/r. One 
can use these diagrams to compute the GW observables of eqs. fl2.33p and (12.361) : for 
example, from the imaginary part of the diagram in Fig. 12.81 follows the leading order 
quadrupole radiation formula 



where the brackets denote a time average, dots stand for time derivatives and Qij is 
the quadrupole moment of the system that sets the strength of the interaction with 
radiation gravitons in the vertices. The time evolution of Qij{t) can be obtained by 
solving the equations of motion for the world-lines Xa(t), which follow from the real 



Before moving to calculating SiyR(h,Xa), it should be noted that eq. (I2.42p does 
not contain, as it stands, the gauge fixing terms necessary to make sense of the path 
integral; this is because the explicit form of such terms is not needed in the subsequent 
discussion. Nevertheless, it is worth stressing that the most convenient way to do the 
gauge fixing is the so-called background field method [91], where the gauge is chosen 
in such a way as to preserve the invariance under diffeomorphisms that transform the 
background metric h^^. If such a gauge fixing scheme is chosen, the action for h^i^ is 
guaranteed to be gauge invariant. This requirement places strong restrictions on the 
form of the EFT that describes radiation. 




(2.43) 



(2.44) 



part of SNR{h, x). 
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2.3.3 The conservative dynamics at first post-Newtonian or- 
der - The Einstein-Ilbert-Hoffmann Lagrangian in the 
effective field theory approach 

With the aim of calculating SNR{h, Xa) as a sum of Feynman diagrams, we now 
proceed to deriving the power counting rules for the interaction terms in eqs. (12 .Mp 
and fl2.18p involving only potential modes, i.e. radiation will not be considered for the 
moment and we postpone its treatment to next section. 



Power counting and scalings for diagrams with only potential modes 

The scaling of -ffk/ii/ in terms of kinematic variables can be obtained by looking at 
its propagator, i.e. at terms in the gravitational action that are quadratic in H\^^jj\ 



= — / dt 



2 J (27r) 



2 



(2.45) 



where = -f^^k- Because potential modes have do ~ vdi, the terms in the second 
hne of this equation are suppressed relative to the first line by f ^ and can be treated 
as perturbative corrections. Therefore, to leading order, the propagator for ifk/ii/ reads 

{TH^^,{x')H^^^{0)) = -^{2'Kf6^{k + q)5(x°)P^,;„^ , (2.46) 

where T stands for time-ordering and the tensor structure is given by eq. ( 12.17p . Be- 
cause of the factor k~^, eq. ( 12.46P is essentially the Fourier transform of the Newton po- 
tential (cfr. with eq. (I2.3ip ). which will be then expressed by a one-graviton-exchange 
like Fig. 12.6( a). Since the tensor structure (I2.17P contains factors of r^^^,, one has 
P^iu■,a|3 ~ C(l) and the scaling of H^^^^ is determined by the other terms in eq. (I2.46p 

Hl^^(^\(^'\Q-'^r% (2.47) 

SO that 

i^k^. -r^y^, (2.48) 

where 5^(k) has been power-counted as (k)~^. In terms of H]^^^, the linear interaction 
between gravity and the matter source of eq. (12.180 reads 



2M, 



PI 



j dfh,, x^x" = j da;°^e^'^-^(^°)ifk^, x^x" (2.49) 
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Therefore, to lowest order in velocity, the matter-gravity interaction vertex scales as 



ra 



m 



-1/2. 



M, 



PI 



(2.50) 



To complete this expression we need the scaling of m/Mpi. this is fixed by means of 
the virial theorem, 

GnTji m?' 



mv r = Lv, 



(2.51) 



where we have used Gat ~ and we have introduced the orbital angular mo- 

mentum L = mvr. The interaction in eq. (12.50p then scales as 



m 



2M, 



PI 



(2^ 



^ik-x(xO) o- rl/2 



(2.52) 



the exchange diagram with two insertions of this term corresponds to the Newton 
potential between the point particles, its scaling would then be as one power of the 
angular momentum 



(2.53) 



This simple result has an important implication: it shows that, once translated in 
NRGR, the perturbative series of PN corrections will be of the type: 



L{l + v^ + v^---) 



(2.54) 



where I have used only even powers of the velocity because the conservative dynamics 
I am dealing with now is time symmetric. Having assessed the scaling of Newtonian 
potential, it is straightforward to see that a correction of order IPN comes from the 
sub- leading terms in the quadratic action for H^^^^ fl2.45p . In fact these perturbations 
give a propagator which reads 



(i/kM.(^°)^qa/3(0))« = -(27r)353(k + q)i-l5(xO)P^,;, 



(2.55) 



where the subscript is to distinguish this propagator from the one in eq. fl2.46p and 
where 



k4 k4 k4 k2 



(2.56) 
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1 » 



2 ^ ! 



Figure 2.9: The Feynman diagram stemming from the sub-leading contributions to 
the propagator: in NRGR, this 0{v^) correction corresponds to the insertion of an 
operator, denoted by ®. 



which is indeed suppressed by with respect to eq. f l2.46p . The diagram for the 
corresponding one-graviton exchange is reported in Fig. 12.91 where the correction to 
the propagator is described by an operator insertion denoted by the symbol (g). 

The propagator is just the first term in the gravitational action: in eq. fl2.14p we 
wrote its analytical and diagrammatic expansions as 



- 2M^i J d^Xy/gR{x) ^ J 



X 



{dhf + 



h{dhf h^dhf 



Ml 



+ 



PI 



+ ■ ■ 



PI 



(2.57) 



where it is understood that each term in the series is suppressed with respect to the 
previous one: now, we can substantiate this statement by looking at the scalings of 
potential graviton self-interactions. In terms of -ffk^j^, the cubic term of the gravitational 
action (12.57^ has the structure 



(2.58) 



where we have used 



d^xexp |i ^Xl'^"^ } " {^^'' 



(2.59) 



The novelty of eq. f l2.58p with respect to the scalings worked out so far is in the factor k^, 
which is due to the fact that the Ricci scalar has two derivatives acting on the metric. 
Therefore, the full scaling of the interaction term f l2.58p is 



1 



>- ~ 

M, 



X 



PI 



X 



(2.60) 
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Figure 2.10: The so-called seagull diagram, which owes its name to the two-graviton 
vertex sourced by matter. In NRGR the seagull diagram represents one of the correc- 
tions to the Newtonian potential of 0{mGN/r) with respect to it. 



which gives the following scaling of the corresponding exchange between two particles 



Vl 



Vl 



(2.61) 



i.e. a term in the two-body potential that is suppressed by f ^ relative to the leading 
order Newtonian exchange diagram. 

The three-graviton potential f l2.6ip is not the only term which is suppressed by 
f ^ with respect to Newtonian potential because of gravity non-linearitiesQ there is 
indeed another diagram which is shown in Fig. 12.101 and comes from the quadratic 
term {hxxY in the matter-gravity action fl2.18p . To lowest order in the velocity, this 
term reads HqqHqq (cfr. eq. f l2.49p ). so the quadratic interaction term and its scaling 
read 



5, 



m 



VP 



Lv (r 
~ X - I X 



dx' / e^(''+^)-(-°)/fkoo(x°)i/qoo(x°^ 

Jk,q 

2 



(^) X (^'^) 



X 



(2.62) 



m \r . 

Working at order e ~ {mGN/r) ~ v'^, one should also consider terms which are 
suppressed with respect to Newtonian potential because of velocity corrections, as it 
is the case for matter-gravity vertices. In eq. fl2.49p we have written down the terms 
coming from the velocity expansion of H]^^,^ x'^x'^ , the lowest order one being ifkoo 
which scales as L^/^ , as determined in eq. fl2.52p . Analogously, one can see that 

H^oi v' ~ L'/\' (2.63) 



(2.64) 



Another vertex scaling as L^l'^v^ comes from expanding the proper time of the a-th 
particle as in eq. ( 12.32P that we report here for convenience 



dTa = dx'^^/T 



V ~ 

a — 



1- -Va 

2 " 



(2.65) 



5. Remember that the PN expansion parameter is e ^ {mGN/r) ^ . 
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Figure 2.11: The Feynman diagrams that, in NRGR, represent the next-to-leading 
order contributions to the conservative dynamics, i.e. corrections to the Newtonian 
potential which are 0{mGN/r), upper row, and O(w^), lower row, with respect to it 
(cfr. Fig. 11.41 in Sec. II. ip . Diagrams taken from ref. [6]. 

taking the f ^ suppression from here, the polarization that couples to a matter source 
can only be i^koo (cfr. eq. (I2.49P ). The resulting vertex completes the point particle 
action fl2.49p to O('y^), to this PN level the linear coupling of a potential mode to a 
world- line reads 

= I V i^i^e^^-''(^")i/kM. i^i" (2.66) 

ziWpi J 



m 



2Mpi J (27r) 

If we combine the various diagrams and interaction vertices we have treated so far, 
we can build five diagrams for two-body potentials that are suppressed by v"^ with 
respect to the Newtonian one. These diagrams are presented in Fig. 12.111 as shown in 
ref. [6] , the sum of the corresponding amplitudes reproduces the famous Einstein-Infeld- 
Hoffman Lagrangian |2T], i.e. the IPN corrections to the two-body non-relativistic 
motion 



L„2 




G%mim2{mi + 1712) 

(2.67) 



2IX12P 

3(v^ + V2) - 7vi ■ V2 



2 , 2\ _ (Vl • X12)(V2 ■ X12 



X12 



where X12 = xi — X2 is the orbital radius and: 

the first term is the leading order relativistic correction to the particles kinetic energies: 
it is reported in round brackets because it does not come from the diagrams of 
Fig. I2.11( rather, it comes from diagrams of the type of Fig. 12. 7[ with the v- 
expansion pushed up to order t>*^; 

the second term comes from the gravity self-interactions of eqs. (12.6 ip and fl2.62p : 
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the terms in square brackets arise from diagrams where the one-graviton exchange is 
accompanied by velocity-suppressed interactions; these are: the u ^ correction of 
the propagator (I2.55P and the velocity-dependent vertices of eq. (I2.63p . 

This constitutes a first step toward the derivation of the path integral over potential 
modes fl2.42p . The EIH Lagrangian f l2.67p contains the gravitational forces between 
the non-relativistic particles which arise from diagrams with no external factors of the 
radiation field. The functional integral in eq. f l2.42p also generates couplings of matter 
to radiation from diagrams with one or more external radiation gravitons which will 
be the subject of next section. 

Before moving to the description of radiation in NRGR, let us take a deeper look 
at scalings. Even if to deal with GWs we only need a classical field theory, it is 
interesting to see what the scalings in powers of h would be. Let us then assess the 
scaling of a process where two world-lines exchange gravitons through a quantum loop: 
by comparing this process with the two-body potentials worked out in Sec. 12.3. ![ here 
we will justify why we did not need to keep this kind of diagrams into account]^. 
The confrontation among diagrams is reported in Fig. I2.12[ where, by means of curly 
lines, we show Newtonian interaction together with the three-graviton potential and 
the quantum loop one. With reference to Fig. 12.121 the following scaling laws can be 
associated to the different contributions: 



m 



PI 



[dtd'kY[6{t)6'^'\k)k-'] 



m 



PI 



[dtd'kY[6{t)6'^'\k)k-'\ 



M, 



-dt6'^'^\k) {d^k)' 
PI 



m 



PI 



[dtd'kY [Sit)S^'\k)k-'Y ^dt6^'\k) (d'k) 



PI 



where the terms in the first square brackets come from the matter-gravity coupling, 
those in the second come from the graviton propagators and finally the terms in the 
third square brackets come from the three-graviton vertices. 

To restore factors of Planck's constant h, one can apply a rule that relates the 
number X of internal graviton lines (graviton propagators) to the number V of vertices 
and the number C of graviton loops 



£ =X- V + 1; 



(2.68) 



then, taking into account that each internal line brings a power of h and each interaction 
vertex a from the interaction Lagrangian, the total scaling for diagrams where the 
only external lines are massive particles is h'-'~^. According to this rule the third 



6. This argument is also discussed in Chap. [3] because contained in the pubhcation of mine that is 
reported there. However, for a more complete treatment of NRGR it nicely fits this chapter, too. 
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2.3 The effective theory for radiation 




(a) Newton 



(b) IPN 



(c) Quantum loop 



Figure 2.12: Contributions to the gravitational scattering amplitude of two massive ob- 
jects. In this picture curly lines indicate gravitational perturbations over flat space-time 
and do not specifically refer to radiation gravitons of the mode decomposition fl2.38p . 
From left to right, the diagrams represent respectively the leading Newtonian approxim- 
ation, a classical contribution to the IPN order and a negligible quantum loop process. 



diagram of Fig. 12.121 involves one more power of h than the first two. The diagram 
with a graviton loop is then suppressed with respect to the Newtonian contribution, 
apart from some powers of v, by a factor h/L <^ 1, whereas the second diagram in 
Fig. 12.121 is a IPN contribution which does not involve any power of h. Equivalently 
one can notice that there is no kinetic term in the Lagrangian for the matter source 
since the massive object is not a propagating degree of freedom, so the IPN diagram is 
not a loop one. As I will discuss in Sec. 13.31 these scaling arguments remain unchanged 
when one adds other fields/particles, like a scalar, and/or when one introduces another 
mass scale [95], provided that the virial relation (12.21) correctly accounts for the leading 
interaction. 

2.3.4 The effective theory for binary systems - Incorporating 
radiation in the treatment 

In the previous section we have started the calculation of the path integral over 
potential modes (I2.42p : the aim is to get to a theory which is valid at scales larger then 
the orbital radius and only contains radiation gravitons. The first step consisted in 
integrating out potential modes from the two-body potentials. To complete the theory 
below the orbital radius we now have to add radiation modes and treat their couplings 
with both potential modes and matter. These couplings are indeed generated by the 
functional integral in eq. (12.421) and it is only by integrating potential modes out of 
them too that we will obtain a theory sensitive to radiation only. 

There is no big conceptual issue in incorporating radiation, as it can be done fol- 
lowing much of the same steps we described in the previous section. For example, one 
can start again from the propagator to assess the scaling of radiation gravitons h: in a 
suitable gauge [6J, this is given by 




(2.69) 
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(/,p) (pO + feO,p + k) 

Figure 2.13: A three-graviton vertex where a radiation graviton (curly line) interacts 
with two potential gravitons (dashed lines). In order to have a definite power counting, 
the radiation graviton should not impart any momentum to the potential graviton. 
This exemplifies the necessity of multipole-expanding the radiation modes as done for 
QCD in ref. [96] (see text for discussion). Figure taken from ref. |83j . 



where the tensorial structure P^u;ai3 is the same as in eq. f l2.17p . Because h represents 
radiation, its momentum will scale as the inverse of the wavelength, so ~ v/r; 
eq. f l2.69p then implies that radiation modes should scale as 

V ~ - . (2.70) 

r 

This rule allows one to power count terms in the action containing the radiation field; 
however, this is not enough in order to obtain an EFT for radiation that has manifest 
velocity power counting. Let us see why through the following example. Consider the 
interaction vertex depicted in Fig. I2.13[ which involves both potential and radiation 
modes. Following the direction of the momentum p^^, after the interaction the potential 
mode has a three-momentum p + k so that its propagator will scale as 

Fig.m^ ^ =^[l-2p-k + ---] ; (2.71) 
(p + k)^ p^ 

because the individual momenta scale as |p| ~ 1/r and |k| ~ v/r, this propagator 
contains an infinite number of powers of v and cannot be properly power-counted in 
the velocity. To ensure that this does not happen, it is necessary to arrange that 
radiation gravitons do not impart momentum to the potential modes, as first pointed 
out in the context of non-relativistic gauge theories in ref. [96]. This is achieved by 
multipole-expanding the radiation field at the level of the action by means of 

V(x, x°) = V(X, x°) + 5x*9i V(X, x°) (2.72) 

+ ^5x^5x^'9i9,V(X,x°) + --- 

where X is an arbitrary point, which in the case of binary systems can be very conveni- 
ently chosen as the center of mass defined by Xcm = J2a^a^a/ ^a^a- The reason 
why eq. (12. 72^ avoids any transfer of momentum from radiation to potential modes, 
thus preserving the velocity power counting, is the following. The expansion f l2.72p 
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Figure 2.14: Feynman diagrams describing the couphng (12.73^ of a radiation graviton 
to the particle word-hnes. 

is a redefinition of the radiation field, which now is only a function of time because 
the space argument in /i(X, x°) is a fixed point in coordinate space (that can be taken 
as the origin of the coordinate system). This means that the couphngs of the radi- 
ation field to either potential modes or point particles will not be Fourier expanded for 
what concerns h(K,x^) (cfr. eq. f l2.82p ): in this way, terms like fl2.7ip will not appear; 
rather, the factors 5x*9j ~ in eq. ( I2.72p will bring the necessary ingredients to form 
the multipoles of the system. 

Let us start by treating the coupling of the radiation field to the particle word-lines: 

- E ^ / V(x,x°)(i;X) ^ (2-73) 

- E ^ / {^oo(x,x°) + 2/io.(x,x°)< + /i,,(x,x°)« + ^t;2/ioo(x,x°)| , 

where we have performed the velocity expansion up to (9(f^) like we did for potential 
modes in eq. fl2.66p . In diagrammatic terms the action f l2.73p is described by Fig. 12.141 
According to eq. fl2.72p . every term in this action has to be multipole-expanded: we 
display these expansions explicitly in what follows comprising terms up to order v"^. 
For the polarization hoo, we have 



+ iXa-3.A„„(0,i-°) 

+ ^Sx^'Sx,'d,aihm(0, i») + ^vlhniO, x°)\ , (2.74) 

where we have separated the various terms in lines according to their scaling in the 
velocity and where we have made explicit the choice X = Xcm = for the origin of the 
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coordinate system. The leading order of eq. fl2.74p is the coupling of hoo to the mass 
monopole of the system M = '^^ma- In eq. ( I2.37P we saw that the mass monopole 
does not lead to the emission of radiation, so this coupling does not contribute to the 
power: indeed, hoo is not a physical degree of freedom in GR; moreover, one can show 
that the mass monopole is a conserved quantity, so it cannot be a source of radiation. 
The second term on the right hand side of eq. fl2.74p is of order v with respect to the 
leading order because of the term 5xa*9j; it can be transformed into 

MX^^9,/ioo(0,xO) (2.75) 

using the definition of the center of mass coordinate to write J2a^a'Ki — ^-^cmi 
because we could choose the origin of the system to be in the center of mass, we have 
= 0, so that this coupling of radiation to matter is zero. The terms in the last 
line of eq. fl2.74p are of order f ^ and cannot be made to vanish. 
For the polarization hoi, we have 



2Mpi 
m, 



2Mpi 



dx°2|/ios(0,x°)< 

+ 5iLjdJios{Q,xl)vl) . (2.76) 



Here the first term is of order v and can be re-written as 

2/io.(0,x°)P^^ = (2.77) 

where we have defined Pcm = ^a'^^'^a total linear momentum of the system 

and used the fact that it is zero in the center of mass frame. This term would have 
contributed to radiation at order v as the second term of eq. (I2.74p : the fact that both 
these couplings vanish corresponds to the absence of dipole radiation in GR. For what 
concerns the second term on the right hand side of eq. fl2.76p . it is of order and 
cannot be made to vanish. 

Finally, concerning the term proportional to the polarization h^g in eq. (12.731) . we 
have 



E 



_fno_ 
2Mpi 



2Mpi 



dxUKs{0,xlyX , (2.78) 



which is of order and does not vanish. 

Collecting the non- vanishing terms of order from eqs. fl2.74p . fl2.76p and (I2.78p . 
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(a) 



(b) 



Figure 2.15: Feynman diagrams representing the radiation sector at first non trivial 
order in GR. See text for discussion. Figure taken from ref. [6]. 



we have the 0{v'^) Lagrangian expressing the coupling of radiation to point particles 

+ Wa,^o.(0,x°)< + /i..(0,x°)t;xl • (2.79) 



As remarked in ref. [6], this Lagrangian has two problems: it is not gauge invariant un- 
der infinitesimal coordinate transformations and seems to predict that the un-physical 
modes /iqo and Jiqs can be sourced by the system. These two problems are related. The 
solution comes from considering the couplings of radiation to potential modes. These 
contributions are represented by the second and third diagram of Fig. I2.15[ which 
encodes the complete radiative Lagrangian at order v"^. Indeed, by means of power 
counting, one can check that 



'H 



/ n 



~ V 



5/2 



V ~ — ^ 



(2.80) 



so that every graph in Fig. 12.151 scales as ^/Lv^I'^ . Because the monopole contribution 
scales as y/Lv^^"^, the diagrams in Fig. 12.151 are suppressed by w^: in fact, as we will 
see below, the third diagram gives rise to the coupling of radiation with the system's 
quadrupole, which is the term of order f ^ in the multipole expansion. 

Let us now go through the diagrams of Fig. 12.151 more in detail. The first graph 
corresponds to eq. fl2.79p and it has been discussed above. 

The second graph comes from a two-graviton coupling of the point particle similar 
to the one that gives rise to the seagull diagram fl2.62p . As in that case, the lowest 
order contribution to a vertex {h x i)^ corresponds to taking the polarization 00 for both 
modes: the only difference is that here we have a potential and a radiation graviton. 
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so that only the first is integrated over. The explicit result reads [6]: 



Fig. EMh) 



2Mi 



PI 



Xi - X2 



00 



(2.81) 



which we will comment about when we sum all the diagrams to obtain the final ex- 
pression. 

For the diagram in Fig. l2.15T c). we need the vertex hH'^, whose expression reads [6] 



^hH'^ — k^^oo 



4 



t/xO 



(2.82) 



+hi 



- 5i,k^ ( -T^^ri^-k„. + ^i^ki^-k 



where the integral over momentum k has been suppressed and where the dependence 
(Xcm,a;°) of hfj,iy has been left implicit. Given this term, it is easy to show that 



-k^/ioo + -k'^hii - kikjhij 



To solve the last of the momentum integrals, one needs a generalized version of the 
Fourier transform fl2.3ip 



-ik-xkjkj 1 



57r X 



Sin 



(2.83) 



using as well as the equations of motion at leading order, the final result is 

Fig. 12.15(c) = ^— [ dx^ 3G'Afmim2 _ ]_y^jy^^-^^^ 
mpi J 2 xi-X2 2^ 



Adding together the results of eqs. fl2.79p and fl2.84p to what one obtains by sup- 
plementing eq. (I2.8ip with its mirror image under exchange of the particle label, one 
finds the complete Lagrangian at order ^/Lv^I'^ 



Ly2 [h] 



2Mi 



00 



PI 



1 X - 2 G]\fmim2 
- > ma^r„ — ; - 



Xl -X2I 



2Mi 



PI 



-(-ijkl-'kdjh, 



i'l-Oi 



2M, 



OiOj ) 



(2.84) 



where: 
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Figure 2.16: A schematic representation of what the EFT for binary systems is about: 
at length scales larger than the orbital radius, and smaller than the typical wavelength 
of radiation, the two-body system can be regarded as a point particle with suitable 
gravitational interactions. Figure courtesy of Chad Galley. 

The first term is the coupling of Jiqq to the Newtonian energy of the two-particle 
system; this term can be regarded as a correction to the mass monopole of the 
source given by the kinetic and gravitational energy 

a a ^ ^ ' ' 

The second term is a coupling of the graviton JiQi to the total mechanical angular 
momentum of the system, L = Xla-^a ^ f^a^a- Both the mass monopole and 
the angular momentum are conserved at this order in the velocity expansion: 
therefore, /loo and do not represent physical contributions to radiation. 

The last term is the coupling of the source moment m^x^x^ to the (linearized) 
Riemann tensor of the radiation field, which reads 

RoiOj = ^ {dohj + didjhoo - d^dih^j - dodjhoi) ; (2.86) 

one can show that Rqq = Roioi = for on-shell graviton matrix elements, so that 
radiation only couples to the traceless quadrupole moment of the source 

Q'' = J2^a(«-l^%) . (2.87) 

The derivation of the radiation Lagrangian is an interesting example of the systematic 
approach of EFT. In fact, because of power counting, the radiation Lagrangian at 
(9(f^) cannot just contain terms from the multipole expansion of the point particles 
couplings to radiation, i.e. it cannot be given just by eq. ( 12.79P : at the same order in 
the velocity one finds contributions from diagrams with graviton self-interactions. This 
is consistent with the strong equivalence principle, according to which gravity couples 
in the same way to all sources of energy-momentum, including the gravitational field 
generated by the source. 

What eq. fl2.84p represents is the matching of the short- distance theory valid below 
the orbital scale r to the long-wavelength theory valid at scales £ ~ r . This matching 
constitutes the second step in the construction of NRGR as a tower of effective theories 
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Figure 2.17: Feynman diagram representing the one-graviton emission amphtude from 
the binary system. At length scales larger of the order of the orbital radius, and smaller 
than the typical wavelength of radiation, the two-body system can be regarded as a 
point particle with suitable gravitational interactions (see eq. ( I2.84p ). 




Figure 2.18: The self-energy diagram whose imaginary part gives the radiated power. 



that we discussed in the introduction of this chapter and is described in Fig. 12.161 as it 
is evident from this figure, the orbital separation between the constituent objects is not 
resolved at this level and the binary system itself is now described as an effective point 
particle. As anticipated, the scale r is only present in the Wilson coefficients of the 
operators of the theory f l2.84p . The Feynman diagrams corresponding to this theory 
are of the type of Fig. I2.17[ where we have used a double-line notation to indicate that 
this is a vertex in the EFT above the scale £ ~ r. The strength of the interaction 
vertex is one of the couplings of eq. (12.841) . 

Once we have calculated the action SNR[h, (xa)], the EFT prescription requires that, 
in order to compute observables, we integrate out the radiation modes and obtain a 
theory where the only degrees of freedom are the particle world-lines. As discussed at 
the end of Sec. 12.3.21 this action can be obtained by summing the Feynman diagrams 
that have no external graviton lines like the one of Fig. 12.81 that we report here for 
convenience. Now we know that the interaction vertices in this diagram are those of 
eq. (I2.84p . among which only the last one is a physical radiative coupling: taking the 
imaginary part of a diagram with two such couplings in the vertices, one can then 
derive the leading order quadrupole radiation formula. An alternative to taking the 
imaginary of diagrams like the one in Fig. 12.181 is to compute the square of the one- 
graviton emission amplitudes, whose corresponding Feynman diagram is the one in 
Fig. 12.171 This is indeed what has been found in eq. (I2.37P by direct calculation and 
replaces de facto the explicit computation of the effective action Seff{xa)- 

To end this section I would like to emphasize that, for the purpose of characterizing 
a GW signal, it is not enough to consider the diagrams of Fig. I2.15[ i.e. the dissipative 
dynamics of the system. In fact, future GW detections at interferometers will probe the 
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non-hnear dynamics of GR in both its conservative and its radiative regime through 
monitoring the phase of the waveform: as I have aheady shown more than once by 
nowEI, a typical expression for the phase is the one corresponding to an extreme-mass- 
ratio 

^^Gwit) = fdrujawir) = ^ T dv'v''^§l^ , (2.88) 
J to LtnM i^[v ) 

where it is evident that the phase receives contributions from both the energy func- 
tion E and the power emitted P. As a consequence, tests of gravity which involve 
the orbital decay, but are not limited to it, should take both regimes into account. 
This is indeed the case of the parametrized post-Keplerian test (7,a;,P), that I have 
discussed in Sec. II. 3[ where the orbital decay is measured together with parameters 
that do not belong the dissipative dynamics. Another example is represented by the 
phenomenological tests that I have investigated personally. These studies have been 
inspired by the NRGR approach and are discussed in the following two chapters. 



2.3.5 Further examples of diagrams corresponding to known 
effects 

In order to further elucidate the rationale behind the decomposition of gravitational 
perturbations in high frequency modes and a slowly varying background (12.38^ . it is 
instructive to present here some Feynman diagrams that correspond to physical effects 
pertaining the different regimes of conservative and radiative dynamics. An interesting 
case is the one related with spin phenomena. The inclusion of spin degrees of freedom in 
NRGR has been addressed in refs. [971198]: here, the conservative dynamics is studied at 
leading order and the corresponding spin-orbit and spin-spin potentials are reproduced. 
Spin effects are incorporated by adding the world-line degrees of freedom A;^(A), where 
A is the affine parameter of the world-line and A;[ is the boost that transforms the 
locally flat frame, labelled by a, to the co-rotating frame labelled by J. These frames 
are defined by vierbeins, which are given by and = e'^Af for the locally fiat and 
co-rotating frame respectively and verify e'J^e'lg^'^ = 7]°"^. Then one can introduce the 
generalized angular velocity given by VL^^ = e'^'^{Dej/dX). Finally, the spin is defined 
as the tensor Sf^u that is the conjugate momentum to Qfj,^ so that the form of the 
world-line action reads ^8] 

S = -J2[l + / \sl'^\jdx}j , (2.89) 

where the sum extends over the binary constituents and S^"^ = S'^'^e^e;^. From this 
action one can derive the Mathisson-Papapetrou equations of motion for spinning test 
particles [99| 1100] . The spin-gravity coupling in f l2.89p can be rewritten by introducing 
the Ricci rotation coefficients u^'' = elDf^e"-" [97] 

Sspin-grav = f S Lab^'jf'u^dX , (2.90) 



7. See the discussions about eq. p.23|) in Sec. 12.3.1] and eq. (|1.52l) in Sec. 11.41 
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v° iP- 
(a) (b) 

Figure 2.19: Panel (a): Feynman diagrams representing the spin-orbit interaction at 
leading order. Panel (b): Feynman diagram representing the spin-spin interaction at 
leading order. A grey blob represents an insertion of the spin operator in the point 
particle action; the factors of correspond to the suppression with respect to the 
leading order coupling of a potential mode to a world-line (see eq. f l2.50p ). Figures 
taken from ref. [97] . 



with Sf^ = S^^e'^^e^,, the spin in the locally flat frame. Expanding f l2.90p in the weak 
gravity limit one obtains the Feynman rules for NRGR with spin |97i |98j; for example, 
at lowest order, the spin couphng of a potential graviton to a matter source reads 

l'spin = 2M^i^'^^''^'-^ ■ (2-91) 

By making use of similar rules for higher order terms, one can derive the conservative 
dynamics for a binary system in the case of spinning bodies. At leading order the spin- 
orbit and spin-spin interactions are represented by the Feynman diagrams of Fig. 12.191 
as one can see from the figure, with respect to Newtonian potential, the spin-orbit 
interaction is 0{v^) while the spin-spin is 0{v^), i.e. these interactions constitute 
higher PN corrections than the EIH Lagrangian of Fig. 12. lit this explains why we 
did not take spin into account so far. From Sec. 11.31 we know that it has recently 
become possible to measure and constrain the general relativistic spin-orbit interaction 
by pulsar timing, even if with a slightly worse accuracy than other post-Keplerian 
parameters (cfr. Fig. 11.131) . 

As discussed at the end of the previous section, future GW detections at inter- 
ferometers will probe the non-linear dynamics of GR in both its conservative and its 
radiative regimes through monitoring the phase of the waveform. Concerning the dis- 
sipative dynamics, the presence of spin affects the radiative multipole moments, like 
the mass quadrupole lij which equals Qij at leading order. At next-to-leading order, 
lij receives contributions from the spin of both objects: these terms have been recently 
calculated in ref. [8] and the corresponding diagrams are represented in Fig. 12.201 The 
confrontation of Fig. 12.191 with Fig. 12.201 constitutes another manifestation of the dif- 
ferent roles that potential and radiation modes have in NRGR. 

Another interesting example of the complementarity of H and h is provided by phe- 
nomena which are still to be investigated in NRGR: the tail of the radiation, discussed 
in ref. [66], and the non-linear memory effect, predicted in refs. [1011 1102"] . These are 



78 



2.3 The effective theory for radiation 




Si 
-9- 




S2 

a) 



b) 



Figure 2.20: NRGR diagrams describing some of the spin contributions to the quad- 
rupole moment at next-to-leading order. The spin-graviton vertex of particle i is still 
represented with a grey blob in these diagrams it is labeled by S^. Figure taken from 



peculiar non-linearities of GR and their detection would then be an invaluable confirm- 
ation of the validity of the theory. The possibility of testing GR through detection of 
the tail effect has been considered in ref. |103] while the importance of the non-linear 
memory has been recently investigated in ref. [104] . 

The tails of GWs result from the non-linear interaction between the quadrupole 
radiation generated by an isolated system with total mass-energy M and the static 
monopole field associated with M. Their contributions to the field at large distances 
from the system include a particular effect of modulation of the phase in the Fourier do- 
main, which has M as a prefactor and depends on the frequency ujgw as "ugw Incucvi/"- 

On the other hand, the non-linear memory effect is a slowly-growing, non-oscilla- 
tory contribution to the GW amplitude. It originates from the fact that, in a non- 
linear theory like GR, GWs can themselves constitute the source of GWs. In an ideal 
interferometer a GW with memory causes a permanent displacement of the test masses 
that persists after the wave has passed. Remarkably, the non-linear memory affects the 
signal amplitude starting at leading order, i.e. the same as the Newtonian quadrupole; 
however, due to the non-oscillatory behavior, detecting the memory at interferometers 
is cumbersome. 

In both the tail and memory phenomena, gravity non-linearities affect the propaga- 
tion of the emitted waves. Therefore, in NRGR both effects will be represented by 
diagrams where the orbital scale has been integrated out and the outermost field line 
is a radiation graviton h resulting from a non-linear interaction with other gravitons. 
The situation is described in Fig. 12.211 where the arrows on the double lines indicate 
the flow of time for the composite binary source. In the case of the non-linear memory, 
panel (b), the outgoing graviton results from the interaction between a wave emitted 
at time ti and a wave emitted at time ^2 > ti, so the three-graviton vertex involves 
only gravitons of radiative type. For what concerns the tail instead, panel (a), the GW 
produced at a time ^2 back-scatters off the pre-existent background gravitational field§ 
for this reason, one of the gravitons in the three- vertex is indicated with a dashed line. 



ref. |8]. 



8. This long-range static field is due to the two-body system so it is described by a potential mode 
that is not the same as the one integrated out to take care of the orbital dynamics. 
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(a) 



(b) 



Figure 2.21: NRGR description of the non-hnearities affecting the vacuum propagation 
of a GW: panel (a) reports the tail of the radiation, panel (b) the non-linear memory. 
Figures stemming from a discussion with Ira Rothstein at Carnegie Mellon University. 

2.4 Final remarks 

NRGR is adequate for describing the emission of gravitational radiation from bin- 
aries of comparable masses. In this case the perturbative expansion is a series in the 
weak-field/slow- motion parameter e of the usual PN approximation. Comparable-mass 
inspirals are the relevant binary sources in the frequency band of first-generation grav- 
itational antennas; in the future, the space-based interferometer LISA will also be 
sensitive to binaries which are very asymmetric in the component masses, the so-called 
extreme- mass-ratio inspirals (EMRIs). A representative example of these systems is 
constituted by a stellar-mass BH inspiraling a super-massive BH so that their mass ra- 
tio V = mim2/{mi + 7712)"^ is typically of the order of 10~^ and naturally lends itself to 
be used as a perturbative expansion parameter different from the PN e~f^~mGAr/r. 
The light object can be thought of as a test particle moving in the gravitational back- 
ground produced by the super-massive BH so the standard picture in this context is 
BH perturbation theory. Analogously, the EFT approach tailored to describe EMRIs 
will take advantage of the smallness of which is not the case for NRGR. I do not 
present this EFT here since I did not make use of it for the investigations presented 
in this thesis. The interested reader can find more details about it in the works of 
refs. |105[ I106j by Chad Galley who took up its studying in recent years for accurately 
modeling EMRIs in view of LISA. 

To end this chapter, it should be mentioned that the NRGR treatment of the 
dynamics of binary systems has some analogies with the approaches that were tra- 
ditionally used before its introduction, notably the one from Blanchet, Damour and 
collaborators (see the review |107] for a complete treatment and full list of citations). 
In this approach space-time is divided in zones and the matching is done between two 
different expansions which are possible for the gravitational field according to the zone. 
With reference to Fig. I2.22[ the procedure can be sketched as follows. The zone com- 
prising the binary system and its neighbors until a length scale i = Ri is the so-called 
near zone: outside this zone the metric admits the weak-field post-Minkowskian ex- 
pansion which is controlled by Newton's constant G^. The slow- motion PN expansion 
is rather possible below a length scale i = R2. This scale is of the order of the radiation 
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Far zone 




Figure 2.22: A pictorial description of the approach from Blanchet and Damour (see 
the review |107j for a complete treatment and full list of citations). Figure courtesy of 
Chad Galley. 



wavelength, so it is bigger than Ri: there exists a region where both approximations 
are valid which makes it possible to match the two expansions. 

For what concerns the differences between the NRGR framework and the standard 
approach one could remark what follows. First of all, NRGR operates at the level 
of the action, which is a scalar under Lorentz transformations; on the contrary, the 
standard approach works with the equations of motion and has to worry about vector 
and rank-two tensor indices. The standard approach then matches metric components, 
while NRGR matches actions. The former keeps higher order terms until the end when 
they drop out; the latter can rather take advantage of an explicit power counting that, 
ahead of any calculation, enables one to select the interaction terms needed at the 
desired perturbative order. 



Chapter 3 



Classical energy-momentum tensor 
renormalization via effective field 
theory methods 

In this chapter I present the work I have conducted with Riccardo Sturani, who has 
been part of my group at the University of Geneva during my PhD studies. In this 
work |108] we applied the EFT of ref. [6] to study point-hke and string-like sources in 
the context of scalar-tensor theories of gravity. Within this framework we computed the 
classical renormalization of the energy-momentum tensor at next-to-leading order: for 
what concerns helicity-2 gravitons, this means up to first PN order, whereas regarding 
helicity-0 modes this amounts to the (non- derivative) trilinear terms in the interac- 
tion. In such a way, we could write down the corrections to the standard (Newtonian) 
gravitational potential and to the extra-scalar potential for both types of sources we 
investigated. In the case of point particles, we obtained the usual results concerning 
the extra scalar ip, which we endowed with a mass m^: at momenta q higher than the 
mass, the effective potential mediated by ip has a logarithmic profile, rather than the 
behavior typical of IPN terms in Einstein gravity; at low momenta, q <C m^, the 
Yukawa suppression takes place as usual. In the case of one- dimensional objects, we 
calculated the renormalization of the tension for cosmic strings. Before us, this issue 
had been investigated in ref. |109j . by Dabholkar and Harvey (DH), and in ref. |110j . by 
Buonanno and Damour (BD). These two works conducted different analysis, leading to 
apparently conflicting results for the string-tension renormalization. Even if this dis- 
crepancy had already been explained in ref. |110j , we found it interesting to re-analyize 
the whole subject in the context of NRGR. Now, to enter the subject in more detail, I 
will partially follow the introduction given in BD. 

The importance of cosmic strings is due to the fact that the existence of these 
one-dimensional topological defects is a generic feature of grand unified theories and 
of cosmological models motivated by string theory. It is believed that cosmic strings 
form abundantly during phase transitions in the early universe and that they represent 
relevant sources of GWs in the form of a stochastic background 1112] . In this 
context, a significant scientific result of first generation ground-based interferometers 
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has been accomplished in 2009: in fact, the LIGO- Virgo collaboration could put a direct 
upper bound on GWs of stochastic origin 12] in the same frequency region of integral 
constraints previously obtained from Big Bang nucleosynthesis |112l [3] and from the 
cosmic microwave background (CMB) [4]. The situation is described in Fig. 13. II taken 
from ref. p]: we refer to this work for complete discussion and full list of references. 

Loops of oscillating cosmic strings can produce other types of radiation than the 
gravitational one, according to the fields they couple to: for example, one can consider a 
dilaton scalar $ and an axion field B^i, at the same time as the gravitational field g^y. 
All these fields back- react on the string source, thus leading to a self-interaction of 
the string with itself. In the limit of an infinitely thin string, the self-interaction is 
divergent, much in the same way as the self-interaction of a point-like electron with its 
own electromagnetic field. In this last case, Dirac noticed that the divergence could be 
cured by renormalizing the mass, i.e. by assuming that the bare mass of the electron 
depends on a cutoff radius 5 [113] : 

m{5)=mR- — , (3.1) 

where m{5) is the UV divergent bare mass of the electron, mpi its renormalized mass 
and e its electric charge. Notably, this cutoff dependence is consistent with the as- 
sumption that m{5) represents the total mass-energy of the particle plus that of the 
electromagnetic field contained in a sphere of radius 5; in other words: 



m{52) - m(5i) = + £x T°° (3.2) 

J6i 



where T^ieidi^) = -E(r)^/(87r) = e^/(87rr^), with r the radial coordinate. If one gener- 
alizes this argument to the case of a string in four dimensions, the linearly-divergent 
electron mass translates as a logarithmically-divergent string tension 

/i(5)=/iR + Clog(^^^ , (3.3) 

where C is a renormalization coefficient that we will discuss next and is an arbitrary 
length which has to be introduced because of the logarithmic dependence of /i on 5. 
The coefficient C receives contributions form each fundamental field with which the 
string interacts: if one considers the case of the gravitational field g^^, the dilaton 
scalar $ and the axion field B^^^, C is given by 

C = Cg + C^ + CB. (3.4) 

In the spirit of Dirac's argument, the individual renormalization coefficients Cj's would 
have to be determined by the energy of the corresponding field. This is the method 
used by DH and us: it gives the correct vanishing value for C in the case of fundamental 
strings, for which a peculiar link exists among the couplings of the various fields to the 
string (see the discussion below eq. (I3.32p ). On the other hand, the individual values 
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of the Cj's obtained by DH and us do not match the ones calculated by means of 
other approaches in the literature, notably that of BD. This is due to the fact that BD 
calculated the renormalization of the string tension by means of a different quantity, as 
explained by BD and also by us in ref. |108j . While DH and we computed the energy- 
momentum tensor for a string, BD derived its effective action. In the latter case, by 
taking the coincidence limit in 4-D coordinate space, BD found divergent contributions 
to the total interaction energy which are localized on the source. These terms are not 
obtained within the approach based on the energy-momentum computation. However, 
the effect of the divergent pieces is to contribute to the renormalization of the bare 
unobservable tension of the string: therefore, on physical grounds, the results of the 
two approaches are equivalent. 

In the EFT terminology that we used, this difference can be grasped more clearly 
once expressed in terms of Feynman diagrams. DH and we calculated a quantity 
that corresponds to Fig. 13. 6[ where there is an external line corresponding to the 
gravitational field; on the contrary, BD obtained the amplitude of the self-interaction 
process depicted in Fig. 13.71 A more complete discussion is contained in our work that 
I report next. In this publication a different notation is adopted with respect to the 
rest of the thesis for what concerns fields in Feynman diagrams: here, curly lines are 
used to represent the scalar degree of freedom of the gravitational field (see eq. (I3.10p 
for the metric parametrization) ; whereas a double line here describes a string source, 
instead of an unresolved binary system as in Chapters [2] and |H 
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Figure 3.1: Comparison of different measurements and models of stochastic back- 
grounds of GW. The upper bound derived in ref. [2] is the one pointed at by the arrow 
and refers to LIGO's [38] fifth scientific run S5: at 95% confidence level, it is given 
by Qq < 6.9 X 10~^ and applies in the frequency band 41.5-169.25 Hz. This bound 
should be compared with the previous LIGO S4 result and with the projected sensit- 
ivity of Advanced LIGO [78]. The bound due to BBN [112[ [3] is an indirect integral 
bound: it applies over the frequency band denoted by the corresponding horizontal 
line. A similar integral bound (over the range 10^^^ — lO^'' Hz) can be placed using 
CMB and matter power spectra [1]. Projected sensitivities of the satellite-based CMB 
experiment Planck |114] and the space-based GW interferometer LISA [30] are also 
shown. The pulsar bound |115] is based on the fluctuations in the pulse arrival times 
of millisecond pulsars and applies at frequencies around 10~^ Hz. Measurements of the 
CMB at large angular scales constrain the possible redshift of CMB photons due to the 
stochastic backgrounds of GW: this limits the amplitude of the backgrounds at largest 
wavelengths (smallest frequencies) [3]. Examples of inflationary, cosmic strings and 
pre-Big-Bang models are also shown: however, the amplitude and the spectral shape 
can vary significantly as a function of model parameters. Figure taken from ref. [2]: 
refer to it for complete discussion and full list of references. 



Classical energy-momentum tensor renormalization 



85 



General Relativity and Gravitation 42, 2491-2509, (2010) 

Classical energy- momentum tensor renormalization 
via effective field theory methods 

Umberto Cannella and Riccardo Sturani 



abstract 
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Goldberger and Rothstein, to study point-like and string-like sources in the context 
of scalar-tensor theories of gravity. Within this framework we compute the classical 
energy-momentum tensor renormalization to first Post-Newtonian order or, in the case 
of extra scalar fields, up to the (non-derivative) trilinear interaction terms: this allows 
to write down the corrections to the standard (Newtonian) gravitational potential and 
to the extra-scalar potential. In the case of one-dimensional extended sources we give 
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one paper and by Buonanno and Damour in another, already discussed in the latter. 
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3.1 Introduction 

We consider in this work the classical renormalization of the energy-momentum 
tensor (EMT) of fundamental particles and strings due to their interaction with long 
range fundamental fields, including standard gravity. The gravitational self-energy of 
a massive body for instance, arises because of gravitons' self-interactions, it can be 
described as an effective renormalization of the massive body EMT and it is fully 
classical having its analog in Newtonian physics. Such self-interactions, even if they 
involve point-like particles, are not divergent when gravity is present, as on general 
grounds General Relativity imposes a lower limit on the size of massive objects: their 
Schwartzchild radii. 

In the case of one- dimensional extended objects like strings, no horizon analog is 
present and no fundamental lower limit can be imposed on their size: classical contri- 
butions to the EMT due to self-interactions of gravity can (and do indeed) diverge in 
this case. Letting the source size shrink to zero and keeping fixed other physical para- 
meters like mass and charge (and eventually neglecting gravity), usually one encounters 
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infinities, or equivalently, physical quantities depending critically on the source size. 
Dirac [113] emphasized that the cutoff dependence of the energy of the electromag- 
netic field sourced by an electron can be absorbed by an analog dependence of the 
bare electron mass, to provide a finite, physically observable invariant mass. However 
the usual way to calculate mass renormalization is by considering the virtual process 
of emission and reabsorption of a massless field, like for mass renormalization of the 
electron in standard electrodynamics, rather than a renormalization of the EMT, i.e. 
of the particle coupling to gravity, as we are going to do here. The above mentioned 
virtual processes are usually considered in the context of quantum field theory, but 
they show their effects also classically, when heavy, non-dynamical, non-propagating 
sources are considered, as we will show. 

In order to compute these quantities we make use of the formalism introduced in 
[HI [83], which is an effective field theory (EFT) method borrowed from particle physics, 
where it originated from studying non-relativistic bound state problems in the context 
of quantum electro- and cromo-dynamics [71 1116] : for this reason, it has been coined 
Non Relativistic General Relativity (NRGR) (see also [12j for the first application of 
field theory techniques to gravity problems). Here we apply NRGR in the framework 
of scalar-tensor theories of gravity for computing next-to-leading order corrections to 
the EMT renormalization, which in turn define, via the usual Einstein equations, the 
profile of the graviton generated by the sources. 

An example of such a renormalization has been worked out in [117] for point 
particles in the GR case and by [109[ 1110] for string-like sources coupled to an extra 
scalar, the dilaton, and an anti-symmetric tensor, the axion. See also [118[ I119[ 1120] for 
the string sources interacting with axionic and gravitational fields. We find particularly 
worth of interest the different analysis performed in [109[ 1110] , leading to apparently 
conflicting results for the string-tension renormalization. The explanation of the dis- 
crepancy is actually given already in [110] , but here we re-analyize it with the fresh 
insight available thanks to NRGR. 

The plan of the paper is as follows. In sec. I3.2l we summarize the basic ingredients of 
NRGR and set the notation for the case at study. In sec. I3.3l we apply EFT methods to 
a model where a scalar and the standard graviton field mediate long range interactions, 
to compute the effective EMT of a massive body. In sec. 13.41 we present the analogous 
computation for a one-dimensional-extended object in four dimensions. Finally we 
draw our conclusions in sec. 13.51 

3.2 Effective field theory 

We start by describing the basis of NRGR: in doing so we closely follow the thorough 
presentation given in [6], to which we refer for more details, with the exception of the 
metric signature, as we adopt the "mostly plus" convention: ri^i, = (— , +, +, +). 

In order to be able to exploit the manifest velocity-power counting, which is at 
the heart of PN expansion, we must first identify the relevant physical scales at stake. 
If, for simplicity, we restrict to binary systems of equal mass objects it is enough to 
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introduce one mass scale m and two parameters of the relative motion, namely the 
separation r and the velocity v. It turns out that, up to the very last stages of the 
inspiral, the evolution of the system can be modelled to sufficiently high accuracy by 
non-relativistic dynamics, i.e. the leading order potential between the two bodies is 
the Newtonian one. The virial theorem then allows to relate the three afore-mentioned 
quantities according to 

2 Gjv^ _s 

V ~ (3.5) 

r 

(where Ga? is the ordinary gravitational constant) and tells that an expansion in the 
(square of the) typical three- velocity of the binary is at the same time an expansion in 
the strength of the gravitational field. 

The compact objects being macroscopic, they can be considered fully non-relati- 
vistic {v << c) so that from a field theoretical point of view, and with scaling arguments 
in mind, the binary constituents are non-relativistic particles endowed with typical four- 
momentum of the order ~ (£■ ~ mv^, p ~ mv) (boldface characters are used to 
denote 3- vectors). Concerning the motion of the bodies subject to mutual gravitational 
potential, it is convenient to consider only the potential gravitons, i.e. those responsible 
for binding the system as they mediate instantaneous interactions: their characteristic 
four-momentum will thus be of the order 

V 1 
~ (^° ~ k ~ -) 

so that these modes are always off-shell {k^k^ ^ 0). 

When a compact object emits a single graviton, momentum is effectively not conserved 
and the non-relativistic particle recoils of a fractional amount roughly given by 

\6p\ Ikl h 

IpI ~ IpI ~ ^ ' 

where L ~ mvr is the angular momentum of the system: it is clear that for macroscopic 
systems such quantity is negligibly small. To summarize, an EFT approach describes 
massive compact objects in binary systems as non-dynamical, background sources of 
point-like type: quantitatively this corresponds to having particle world-lines interact- 
ing with gravitons. The action we consider is then given by 

S = Seh + Spp , (3.6) 

where the first term is the usual Einstein-Hilbert action 

Seh = 2M|, J d'x./^ R{g) , (3.7) 

with the Planck mass defined (non canonically) as Mpf = 327rGiv — 1-2 x lO^^GeV, 
and the second term is the point particle action 



S^ = -n.jir = -r.j^^J::X, (3i 
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in which (yf^j, is the metric field that we write as gf^u = fltif + ■ To make the graviton 
kinetic term invertible, one should also include a gauge fixing term like 



9f 



-M 



PI 



(3.9) 



with T^ = d''h^,-l/2d^K. 

We now parametrize the metric following |121] . instead of [6j, as 



(3.10) 



where = 0,..,3 and i,j = 1,2,3. We define 7*-^ as the inverse matrix of ■jij, so 
that 7'-' = and a* = 7*%^ . It is also useful to introduce = •jij — 5ij (so that 

= to first order) and <j = QjS^^ . Then, to quadratic order, the following action 
for non-canonically normalized fields is obtained 



_|_ c 

EH I quadratic '"^9/ 



(3.11) 



X 



The non-relativistic parametrization of the metric f l3.10p allows to write down all 
the terms that do not involve time derivatives in a simple way 



EH static 



2M\ 



PI 



dt rf^xi 



-7 



i?(7) - 2d,ipd,ip^'^ + -e^'^F.^F^f 



, (3.12) 



where F^j = diUj — djai is the usual field strength tensor. 
The canonically normalized fields cTjj, 0, Ai can be defined as 

Ai = V2Mpi ai . 



(3.13) 



The only interaction term we will need, as it will be explained, is the cubic one (T0^ 
given by 



S 



1 



EH\ 



2M. 



PI 



dt c/^x 



dicpdjcf) [ 6ikSji - ^Sij6ki ] (Tki 



The world-line coupling to the graviton thus reads 
S^^-mjdr^-mj^^;^^ 

= -m I dt e'l'/^^^^'^^h 



(3.14) 



(3.15) 



A, 



V2Mpi 

-m J dt e^/(2v^*^«) (^l - 1^2 ^ 



f,-V24>/Mpi^..yiyj 

V + . . . 



2^/2 M PI V2Mpi 
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(a) Newton (b) IPN (c) Quantum loop 

Figure 3.2: Contributions to the scattering amplitude of two massive objects. From 
left to right the diagrams represent respectively the leading Newtonian approximation, 
a classical contribution to the IPN order and a negligible quantum 1-loop diagram. 



The propagators we use are given by the following non-relativistic expressions, as we 
are treating the time derivatives in the kinetic terms as perturbative contributions, 

I ' K 

A,{t,k)A ,{t',k') = {27rf6it-t')6^^\k-k'):^6,, 

</.(t,k)0(t',kO = (2n)'6{t-t')6^'\k-k'):^ 
I I k^ 



where 



Pij,ki = - {SikSji + SiiSjk — 2SijSki) 



(3.17) 



As far as we are only concerned in scaling we can set k ~ 1/r, t ~ r/v and, by 
virtue of the virial theorem (13. 5p . m/Mpi ~ \/Lv. We can then immediately estimate 
what are the scalings of the contributions to the scattering amplitude of two massive 
objects: each of the three diagrams reported in fig. l3.2[ for instance, contributes to such 
process. By assigning a factor [^^^dt d^k] to a graviton-worldline coupling not involving 

velocity, a factor [(5(t)5'^^)(k)k~^] for each propagator, and a factor [-^dt (5(^)(k) (rf^k)^] 
for a three-graviton vertex, the following scaling laws can be associated to the different 
contributions of fig. | 
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3.3 Effective EMT in scalar gravity: the point particle case 



Even if we are actually dealing with a classical field theory, it is interesting to give a 
look at the scalings in powers of h. To restore ^'s one can apply the usual rule that 
relates the number X of internal graviton lines (graviton propagators) to the number V 
of vertices and the number £ of graviton loops 



then, taking into account that each internal line brings a power of h and each interaction 
vertex a, from the interaction Lagrangian, the total scaling for diagrams where the 
only external lines are massive particles is h'-'^^. According to this rule the third 
diagram of fig. 13.21 involves one more power of h than the first two. The diagram with 
a graviton loop is then suppressed with respect to the Newtonian contribution, apart 
from some powers of v, by a factor h/L <^ 1, whereas the second diagram in fig. 13.21 
is a IPN contribution which does not involve any power of h. Equivalently one can 
notice that since the massive object is not propagating (there is no kinetic term in 
the Lagrangian for such a source), the IPN diagram is not a loop one. These scaling 
arguments remain unchanged when other particles are added, like a scalar field, and/or 
another mass scale is introduced [95], as we will discuss in sec. 13.31 provided that the 
virial relation (13. 5 p correctly accounts for the leading interaction. 

3.3 Effective energy- momentum tensor in 
scalar gravity: the point particle case 

The usual way to obtain an effective action F out of a fundamental action Sfund 
is by integrating out the degrees of freedom we do not want to propagate to infinity 
according to the formal rule 



where $ denotes the generic field to integrate out. 

In practice this non-perturbative integration is replaced by a perturbative compu- 
tation, performed with the aid of Feynman diagrams like those of fig. 13.21 which shows 
some contributions to the effective action of two particles interacting gravitationally. 
At lowest order (Newtonian interaction) the diagram in fig. 13.2( a) represents the term 
responsible for the Newtonian 1/r potential between two massive objects. Stripping 
away one of the two external lines in this diagram an amplitude for the coupling of a 
single particle to a graviton is obtained: this amplitude is linear in the external graviton 
wave-function and defines the effective EMT of the particle. Thus, at Newtonian level 
the two diagrams in fig. 13.31 give the following contributions to the effective action 



£ = T- V + 1; 



(3.18) 




(3.19) 




(3.20) 
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Figure 3.3: Feynman diagrams describing the gravitational contributions to the ef- 
fective energy-momentum tensor of a particle at Newtonian level according to the 
parametrization f l3.10p used for the metric. 

where Xp is the three-vector of the position of the source particle and use has been 
made of the Newtonian value of the EMT defined as usual as 

-2 6S 



(3.21) 



Note that the contribution from v'^a^ is vanishing as the aij part of the metric field 
does not couple directly to a static massive source for which Tij{x) = 0, Tqq^x) = 
mS^^\-x — Xp). 

The second diagram in fig. l3.2l is a representative contribution of the IPN corrections 
to the Newtonian potential between two particles. Stripping away again one of the two 
external particle lines the diagram showed in fig. 13.41 is obtained, whose contribution 
to the effective action at next-to-leading order is 

rl^) =-L_ [ d^xaijix)T'^'\x) 



M, 



PI 



Mpi J (27r) 

^Mli J {2Tif {2'Kf k2 (k - q)H ' 

/■ dt 7^a.,(t, -q) (-5^^q + ^ ) e^^^^^ , (3.22) 



2iOM|,J (27r)3^'^'''' ^ " ^ ' q 

where q = ^/q ■ q and we have used eqs. fl3.5ip . The analogous quantity for vanishes 
as there is no 0^ vertex, see eq. (13.121) . Incidentally, we note that the EMT obtained 
from eq. (13.221) is transverse, consistently with the request that the effective EMT has 
to be conserved order by order (see |122j for an interesting discussion of scalar gravity 
at interacting level). 

Another check of the correctness of our result can be obtained by reconstructing 
the metric out of this effective EMT. The linearized equations of motion for gravity 
give 

0(t,k) ' 



k2 50(t,k) 



(3.23) 
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Figure 3.4: Feynman diagram describing the gravitational contribution to the effective 
energy-momentum tensor of a particle at first post-Newtonian order according to the 
parametrization used for the metric fl3.10p . 



which, using the first of eqs. fl3.53p . allows to compute the metric component ip accord- 
ing to 

, , (i){x) m f d^k e'^'^^^'^p^ , 
y^(^) = ^ r^.. = -^TJT / TTT^ 772 = (3-24) 



(3.25) 



where Gn has been reinstated in the final result and r = |x — Xp|. Analogously, for 
one has 

11 6S 



which, again using eqs. (I3.53p . leads to 



^ij (tj 2;) -Pj 



ij;kl 



m 
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PI 



kkh 



<t>=0=<;ki 



p-ik-(x-Xp) 



k k^ 



Given the metric parametrization fl3.10p we obtain 

\2 



900 = - 

901 = 



^ 2GNm _ ^{Gjrmy 



9ij 



2GN'm (Gjvm)M {GNmf x'x^ 

1 + — + .2 ) ^-j + —;r2 



(3.26) 



(3.27) 



which is the Schwarzschild metric to IPN order in the harmonic gauge, see [117] . 

Let us now consider an extra degree of freedom with respect to ordinary gravity, 
that is a massive scalar field ip whose action is given by 



1 



(3.28) 
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y 



Figure 3.5: Feynman diagram representing the self-interaction contribution of the 
massive scalar field ip to the energy-momentum tensor of a particle at next-to-leading 
order. 



where a cubic self-interaction has been allowed. The interaction with the gravitational 
field (Tjj, embodied by the trilinear term ipipa, can be derived from the kinetic term, 
namely 

S^\^^^ = ^J dtd'^d^^d,^ [a^^ - \5^^a^ . (3.29) 

There are no trilinear terms such as (pipip or (fxpip because of the specific metric para- 
metrization we chose fl3.10p . The field ip is assumed to couple to matter in a metric 
type in analogy with f l3.15p : 

S'^^ = _^e"^/(2V2A/pO j ^ 

for some dimensionless parameter a. Therefore the tree-level coupling of to matter 
at lowest order is very similar to the diagram on the left of fig. 13.31 

p{o) ^ am 
2V2M1 

At next-to-leading order we have two possible contributions. The first comes from a 
diagram like that of fig. 13. 4 1 where the two 0's are replaced with two ip^s: the amplitude 
is almost the same as eq. (13.221) . apart from an extra factor a^. The second contribution 
comes from the cubic ip self-interaction, depicted in the diagram of fig. 13. 5t 

Note that at high momentum transfer [q ^ m^) the integrand goes as q^^, whereas in 
the gravity case (13.22^ we had T^j{q) (x q: this difference leads to an effective potential 
due to the if) mediation which has a logarithmic profile, rather than the behavior 
typical of IPN terms in Einstein gravity derived in ^95j; at low momenta (g ^ m^) the 
Yukawa suppression takes place as usual. 



V'(t,Xp(t))rft. 



(3.30) 
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3.4 Effective energy- momentum tensor in 
scalar gravity: the string case 

In the case of a one-dimensional extended source we consider the Nambu-Goto 
string with action Sg given by 

Ss=iJi j e"*/(^*'^0 cirda / do.x^^dpx" e""^ B^, drda , (3.32) 

is 2y/2Mpi is 

where 7 = det7a/3, with 7^^ = daX^d^x^ g^^, x^ are coordinates in the 4- dimensional 
space, a and r are the coordinates on the world-sheet E spanned by the string in its 
temporal evolution. Such an action describes a fundamental string interacting with 
gravity via a string tension /x, with a scalar field $ through a coupling a^/ {■\/2Mpi) 
and with the antisymmetric tensor B^j^y through the coupling l3^/{2\/2Mpi). In this 
notation a supersymmetric string corresponds to a = /3 = 1. 

The convention for indeces is the following: denote the two directions parallel 

to the world-sheet while /x, z/, . . . are generic 4-dimensional indeces, then Latin letters 
... denote 3-space indeces and we will use a, 6 or c to denote the (two) spatial 
dimensions orthogonal to the string. 
The action Sf determining the dynamics of the fields is 



(3.33) 



where H^^p = d^Byp + dpB^y + dyBpp_. The only new propagator we will need with 
respect to the point-particle study is 

Bp,{t,\i)B p„{t\ k') = i {ripp'n,, - r]^^7].p) (27r)' 5(1 - t')5^''\k - k')^ • (3.34) 

Analogously to diagrams in fig. 13. 3t the effective action for the linear coupling to the 
string source of the fields 0, aij, $ and Bp^ is T^o) = rj*^ + r^] + tJ*^ + with 

rj*^ = / 0(Too + T,,5*^) d^x = 0, 

= / QjT'^ d'^x = -YT- / (^u{x{T,a))dTda , 

^ -^^ (3.35) 

= /" $ (Too - Ti.S'^) d^x = [ ^(x(t, a)) drda , 

= ^^1^ I d^x^dpx^e-^Bpy drda = ^ Bo,{x{r, a)) drda , 
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Figure 3.6: Diagrams reproducing the coupling to </> (curly line) and to 0"^ (long- 
dashed), or the effective energy- momentum tensor, of a string at next to lowest order 
in interaction. The diagram on the left vanishes (see discussion in the text). 



where use has been made of the explicit parametrization of a static string: = r, 
= a, and of the definition (13.211) for the string EMT T^^ giving 



T^, = diag(/i,-/i,0,0)(5(2)( 



X 



(3.36) 



Following the same reasoning as in sec. 13.31 the contributions to the renormalization 
of the effective EMT due to the dilaton and the antisymmetric tensor interaction can 
be computed, see fig. 13. 6[ We thus restrict to those trilinear interaction terms involving 
a graviton field, either a or a a, as an external line (in a completely analogous way 
the renormalization of the $ and B^^, coupling could be computed). We then have: 



—3— / dtd^^l- 
2MpiJ [2 

1 

+ 2 



ditdjt I 5"S'"' - -S''5"" 1 a,. 



(3.37) 



where we have defined the tensor 



and we have specified the antisymmetric tensor polarization indices to "01" , as this 
is the only polarization involved in this interaction, and omitted rewriting the terms 
coming from the pure gravity sector, i.e. and 0^0", because they read the same as 
in (fSUl). 

The diagram on the left in fig. 13.61 is actually vanishing because no (f) can attach 
directly to the string and no trilinear term with only one (f) is present in the action (13.331) . 
as it can be seen from fl3.12p or (I3.37p : this implies that the relation Tqo = —Tij6^^ holds 
also at next-to-leading order. We are thus left with the diagram on the right in fig. 13.61 
where the particles propagating in the internal dashed lines can be either two dilatons 
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or two antisimmetric tensors or two gravitons of the type aij. The contribution to T^J^- 
from the diagram involving two dilatons is 



2^,2 



Mpi 



PI 



dr 



(2^ 



(27r)' P (A; - g)^ 



Jab 



g2 



(27r)^ 



+ ( C^ - 2 ) 



(3.38) 



with C a divergent quantity, coming from the last integration in the first line, whose 
value can be read from eq. fl3.54p 



C = lim — 



l + -(7-2 + log [gV(47r)]+o(e)) 



(3.39) 



here dimensional regularization has been used, as this entry of the effective EMT is 
expected to be (logarithmically) UV divergent, see e.g. |109t IllOj . Note that the 
divergent constant only enters the Tn component of the effective EMT. 



For the B^^ interaction a simifiar result is obtained 
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(3.40) 



The contribution to the IPN effective action due to purely gravitational process, 
i.e. by the diagram on the right of fig. I3.6l with three cr's, can be computed by making 
use of the three graviton point function: 



{aii{ki)aii{k2)(Jij{q)) = -^6'-^\ki + k2 + q)q'^Sii6ji , 



(3.41) 



which has been obtained thanks to the Feyncalc tools |123j for Mathematica; the result 
is 
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AMI, 



j dr^,e-^^^H^^6''a.,{r,q) J 
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'Daii{T, q) 
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where D is a divergent constant, again entering the Tn component only, given by 



(3.42) 



D = lim - 

e-S>0 e 



l + -(7 + ln [q'/iAn)]) 



(3.43) 
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Figure 3.7: Feynman diagram representing the string tension renormalization as com- 
puted in |110] . The internal wavy line stands for all possible fields interacting with the 
string: dilaton, antisymmetric tensor and graviton of type a. 

The conserved effective EMT is thus given by the sum of the three contributions 
just calculated and reads 



together with Too = — Tii and Toj = . The coordinate space counterpart of (I3.44p is 
reported in the Appendix. 

We note that in the directions orthogonal to the string the EMT is still vanishing 
for cP' = thus preserving the no-force condition valid for supersymmetric strings 
of the same type (charge). The divergent part of the entry Tu is also vanishing in 
the supersymmetric case due to a cancellation among the different terms: therefore, 
the superstring tension, given by Tu, does not receive divergent contribution. This 
confirms the result of Dabholkar and Harvey |109j obtained through the analysis of the 
EMT's on the (linearized) GR solution around a string. 

In |110j Buonanno and Damour also found a non-renormalization, but via a different 
cancellation. The authors of |110j analyzed a physical quantity which is described by 
a diagram of the type depicted in fig. 13.71 where it is understood that each of the fields 
interacting with the string can propagate in the internal line. We now take a closer 
look at the different contributions to this process. Letting a aij propagate in the wavy 
line of fig. 13.71 yields a vanishing result given that the amplitude for such a process has 
the following behavior 



as it can be explicitly checked from eq. f l3.17p . This diagram vanishes for the same 
reason why two straight, static, parallel strings do not exert a force on each other: the 
amplitude for one graviton exchange between two such strings is proportional to the 
same vanishing quantity Pii;ii. The dilaton contribution to the amplitude of fig. 13.71 is 



(3.44) 




fig. m^r., oc r^aijakiT'^' (x fx^Pii,u = 



(3.45) 




(3.46) 
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whereas to find the effect of the antisymmetric tensor it is enough to replace "a^" 
with "—(3"^" in eq. ( 13.460 . as can be checked using ( 13.340 and ( 13.35^ . These three amp- 
litudes, condensed in the representation of fig. 13.71 have a close correspondence with 
what is found in [110] and show that the contributions to the superstring renormal- 
ization are different when calculated by looking at the self-energy as in |110j other 
than through the (effective) EMT as in [109] and in the present work; nonetheless, the 
non-renormalization property of superstrings is preserved in both approaches. 

The source of the discrepancy is explained in [llOj where it is observed that the dif- 
ference in the two ways of computing the renormalization of the string tension amounts 
to a (divergent) source-localized term, "as the interaction-energy cannot be unambigu- 
ously localized only in the field, there are also interaction-energy contributions which 
are localized in the sources" which are missed in one approach but accounted for in the 
other. Moreover, the contribution of the antisymmetric tensor to the string tension 
renormalization turns out to be the same with the two methods because this coupling 
to the string is metric-independent, so it does not contribute to the total EMT given 
by T'^'^ = 2g~^/'^5S / 5g^u. Of course the physical result cannot depend on the details of 
the calculation method: indeed the source-localized contribution just renormalizes the 
hare tension of the string and does not give physical effects. As observed in [110] , this 
constrasts Dirac's argument [113] about the connection bewteen the renormalization 
of a point charge and its divergent field self-energy. 

Therefore, we support the explanation of the discrepancy given by Buonanno and 
Damour jllOj and provide a computation of the renormalization of the EMT with a 
completely different technique than in Dabholkar and Harvey [109] . confirming their 
result. 

Following the track of the EFT methodes we employed, one could also compute the 
renormalization of the couplings of $ and B^y. For the dilaton coupling the relevant 
diagrams are two, both of the type fig. 13.61 with a $ as outer wavy line and either two 
5,,,,'s or a $ and dashed inner lines. For the antisymmetric tensor case, the 

external B^y can be attached to either a cxjj and a B^y or to a $ and a B^y. All the 
above mentioned trilinear vertices have the same dependence on external momentum 
as the gravity case. 

One final remark is needed about result (I3.44p . A tensor Tab{x) is conserved if 
T^^'^(x) = which, in Fourier space, translates naively to 



d^TM = -iq'^TM NO! (3.47) 

Clearly, with an EMT of the form (l3^ . for ^ (5'^ the right hand side of eq. (I3^ 
does not vanish. This happens because Tab{q) is not square integrable, thus it is not 
ensured that the derivative operation and the Fourier transform commute with each 
other, and indeed they do not in this case, see Appendix for details. 
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3.5 Conclusions 

We have studied point-like and one-dimensional-extended sources in the context of 
scalar-tensor gravity and we have computed the effects of fields self-interactions to the 
renormalization of the effective energy-momentum tensor. 

The calculations have been performed within the framework provided by the effect- 
ive field theory methods applied to gravity [6l [83], exploiting the powerful tool of a 
systematic expansion in terms of Feynman diagrams. 

The classical "dressing" of the sources by long range interactions has the effect 
of smearing the source, consistently with coordinate covariance, and implies energy- 
momentum tensor conservation. We obtained perturbative solutions valid to first post- 
Newtonian order or, in the case of extra scalar fields, up to the trilinear interaction 
terms. 

In the case of a string source we reviewed the renormalization of both its effective 
energy-momentum tensor and its tension, which has been subject of investigation with 
apparently conflicting results in the past |109l 1110] . We exposed the fully satisfactory 
explanation of the discrepancy given by Buonanno and Damour |110] and conflrmed 
that the renormalization of the energy-momentum tensor and the renormalization of 
the string tension differ by source-localized contributions. 
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3.6 Appendix 

To second order the metric (I3.10p can be rewritten as 



where 7,^ = 6ij + <jjj (exact). It is also useful to have the form of the inverse metric 






(3.49) 



To second order one has 




6ij + 2(p6ij - c^ij + 2ip {(p6, 




(3.50) 
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The relevant integrals for computing Feynman diagrams like the one represented in 
fig. 13.41 (see for instance [124] and |117j ) are 



(27r)=^P(k + q)2 
(27r)=^P(k + q)2 
The integral relevant for fig. 13.51 is 

d^k 1 



64 1-^"' + ^ 
16g 



(3,51) 



(27r)M'^^ + ^^)[(k + q)2 + M2] 47rg V2M 
and to reconstruct the metric out of the effective EMT we used 



(3.52) 
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27r2|x|2 V 



(3.53) 



The relevant integrals for computing Feynman diagrams like the one represented in 
fig. 13.61 are (see again |124j ) 
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(3.54) 



qi 
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X [x{l — x)] 2 dx , 
, ,^ — ^Ffl--) / X \x (I - x)]^~^ dx . (3.55) 



(27r)2+^F(fc + g)2 (g2)W2(4^ 

Other useful formulas to anti-Fourier transform the string effective EMT at next- 
to-leading order, are 



rf2g 



27r2 



log(g)e 



iqx 
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(3.56) 



rf2g 1 



-— log(xg,) + 



2^2 



ln2 — 7 r q^ 



l-K 167r -' 



,^(27r)'g2 27r 

where a disk of radius q^ around the origin has been cut out of the integral. Moreover 

"277 



e'^^°'^de = 27rJo(x) 



(3.58) 
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where Jo is the Bessel function of zero-th order. To derive the metric out of the string 
effective EMT the following integral 



It 



<Pq 1 



X 

2^ 



2qlr 



1 li t 2 2\ 7 — In 2 — 1 , 
^ + ^ log (g,r ) + ^ + o{xq. 



is helpful. 

The effective EMT (13.441) in coordinate space is 



4 



(3.59) 



IT 



(3.60) 
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(a2-/32) (C'5(2)(a;'^) + i/r2) 
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— TT^ab H 7r~ 



where r denotes the distance to the string in the transverse two-dimensional space. 
Here C denotes the g- independent part of the quantity defined in text in fl3.39p . 

To explicitly check conservation in the Fourier space of the string effective EMT 
fl3.44p . let us write down the conservation of the EMT in g-space, keeping only the 
components transverse to the string world-sheet: 

d^TM = J d'x [d'^Tatix)] e'^^ = J d^x [9" (T^.e^"^) - tq'^TMe'^''] , (3.61) 

which has an extra piece with respect to fl3.47p . Let us restrict for simplicity to the 
total derivative term and let us fix the index 6 = 2. To make sense of the integral 
we have to integrate over a region Q obtained by cutting out of the plane the the two 
regions r < and r > R, and we will finally (but after taking the other limits first) 
let — )• and i? — )■ 00. 

By changing coordinates from y, z to p,6 according to y = r cos 6, z = r sin 6 and using 
the Green-Gauss theorem one obtains 



TT 



(fx [T«2(a;)e*«-^''] 



'an 2r 



2R 
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(3.62) 



cos9e"'''^'°'^d9. 



e JO 



The first integral is clearly vanishing in the limit i? — > 00. Expanding the exponential 
in the second integral, taking the limit — )■ and finally plugging this result into 
fl3.6ip . one has 



qed. 
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Chapter 4 



Extracting the three- and 
four-graviton vertices from 
binary pulsars and coalescing 
binaries 

In Sec. 11.11 1 have discussed the confrontation of scalar-tensor theories with GR by 
means of tests of relativistic gravity. This confrontation followed the studies of DEF 
and is summarized in Fig. II. 8[ where future GW detections seem not to be competitive 
with binary pulsar observations already at hand. Surprisingly, even the detection of a 
neutron star-black hole system with the most sensitive space-based interferometer LISA 
is found to be less constraining than the electromagnetic observation of the same system 
performed by means of pulsar timing. In the course of presenting these conclusions by 
DEF I have pointed out how they strongly depend on the assumptions made. Notably, 
DEF parametrized a certain class of scalar-tensor theories by considering the couplings 
of matter to the scalar field only up to second order in perturbation. Within this 
choice, DEF confronted different tests of gravity in the radiative regime only to the 
lowest orders probed by binary pulsars, i.e. up to terms of Oiyjcf in the dynamics. 
On the contrary, GW detections will probe gravity up to (9(f/c)^^! Therefore, there 
is the necessity to envisage testing frameworks for gravity in the dynamical regime 
probed by GWs. In Sec. 11.41 I have described one such test where, extending the PPN 
approach to the radiative regime, one applies a phenomenological attitude towards the 
phasing coefficients of the GW formula. Another example of extended PPN tests is the 
one I present in the present chapter. This is a study that I have conducted with my 
supervisor and the rest of his group at the University of Geneva: inspired by the NRGR 
method of ref. [6], that I have described in Chap. [21 we constrained GR non-linearities 
from a field-theoretical perspective |125j . Here I report our publication, while in the 
following chapter I discuss more details in the context of a follow-up of our work. 
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4.1 Introduction 
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Extracting the three- and four-graviton vertices 
from binary pulsars and coalescing binaries 

Umberto Cannella, Stefano Foffa, Michele Maggiore, Hillary Sanctuary 

and Riccardo Sturani 



abstract 

Using a formulation of the post-Newtonian expansion in terms of Feynman graphs, 
we discuss how various tests of General Relativity (GR) can be translated into meas- 
urement of the three- and four-graviton vertices. In problems involving only the con- 
servative dynamics of a system, a deviation of the three-graviton vertex from the GR 
prediction is equivalent, to lowest order, to the introduction of the parameter /SppN in 
the parametrized post-Newtonian formalism, and its strongest bound comes from lunar 
laser ranging, which measures it at the 0.02% level. Deviation of the three-graviton 
vertex from the GR prediction, however, also affects the radiative sector of the theory. 
We show that the timing of the Hulse- Taylor binary pulsar provides a bound on the 
deviation of the three-graviton vertex from the GR prediction at the 0.1% level. For 
coalescing binaries at interferometers we find that, because of degeneracies with other 
parameters in the template such as mass and spin, the effects of modified three- and 
four-graviton vertices is just to induce an error in the determination of these paramet- 
ers and, at least in the restricted PN approximation, it is not possible to use coalescing 
binaries for constraining deviations of the vertices from the GR prediction. 

DDI: 10.1103/PhysRevD. 80.124035 PACS numbers: 04.30.-w, 04.80.Cc, 04.80.Nn 



4.1 Introduction 

Binary pulsars, such as the Hulse- Taylor [13] and the double pulsar [14], are won- 
derful laboratories for testing General Relativity (GR). They have given the first ex- 
perimental confirmation of the existence of gravitational radiation |126[ [60] , provide 
stringent tests of GR and allow for comparison with alternative theories of gravity, 
such as scalar-tensor theories [1271 ED HSU [Ml |371 |36l [IH |33l [El [M] (see [1301 HH EU] 
for reviews). Another very sensitive probe of the non-linearities of GR is given by the 
gravitational wave (GW) emission during the last stages of the coalescence of compact 
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binary systems made of black holes and/or neutron stars, which is one of the most prom- 
ising signals for GW interferometers such as LIGO [38] and Virgo [39], especially in their 
advanced stage, and for the space interferometer LISA [40j. Various investigations have 
been devoted to the possibility of using the observation of coalescing binaries at GW 
interferometers to probe non-linear aspects of GR [371 IM flM im 1761 [T32l fT33l flSl 116] . 

Compact binary systems probe both the radiative sector of the theory, through 
the emission of gravitational radiation, and the non-linearities intrinsic to GR which 
are already present in the conservative part of the Lagrangian. In a field-theoretical 
language, these non-linearities can be traced to the non-Abelian vertices of the theory, 
such as the three- and four-graviton vertices. It is therefore natural to ask whether 
from binary pulsars or from future observations of coalescing binaries at interferometers 
one can extract a measurement of these vertices, much in the same spirit in which the 
triple and quartic gauge boson couplings have been measured at LEP2 and at the 
Tevatron ttSSl [Ml IT37] . 

In this paper we tackle this question. The organization of the paper is as follows. 
In Section 14.21 we discuss how to "tag" the contribution of the three- and four-graviton 
vertices to various observables in a consistent and gauge-invariant manner, and we com- 
pare it with other approaches, such as the parametrized post-Newtonian (PPN) formal- 
ism [19]. In particular, we find that the introduction of a modified three-graviton vertex 
corresponds - in the conservative sector of the theory and at first Post-Newtonian or- 
der (IPN) - to the introduction of a value for the PPN parameter /3ppN different from 
the value /3ppN = 1 of GR. However, a modified three-graviton vertex also affects the 
radiative sector of the theory, which is not the case for the PPN parameter /3ppN- We 
also discuss subtle issues related to the possible breaking of gauge invariance which 
takes place when one modifies the vertices of the theory. In Section H73] we present our 
computations with modified vertices, and in Section 14.41 we compare these computa- 
tions with experimental results obtained from the timing of binary pulsars and with 
what can be expected from the detection of gravitational waves (GWs) at ground- 
based interferometers or with the space interferometer LISA. Section 14.51 contains our 
conclusions. 



4.2 Tagging the three- and four-graviton vertices 

Our aim is to quantify how well the non-linearities of GR can be tested by various 
existing or planned experiments/observations. Historically, there have been several ap- 
proaches to this problem and, basically, one can identify two complementary strategies. 
The first is to develop a purely phenomenological approach in which deviations from 
GR are expressed in terms of a number of parameters, without inquiring at first whether 
such a deformation of GR can emerge from a fundamental theory. An example of such 
an approach is the parametrized post-Newtonian (PPN) formalism. In its simpler ver- 
sion, it consists of writing the IPN metric generated by a source, treated as a perfect 
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fluid with density p(x) and velocity fleld f (x), in the form 



(700 = -l + 2t/-2/3ppNt/% (4.1) 
^70* = -^(47ppN + 3)V^i, (4.2) 
= (1 + 27ppNt/)5ij , (4.3) 



where 

f/(x) = / d'x' ^lin^ , (4.4) 



V.(x)= Id^x'P-P^^ (4.5) 



and the standard PPN gauge has been used [HI ETj (we use units c = 1). General 
Relativity corresponds to /3ppN = 1 and 7ppn = 1- More phenomenological parameters 
can be introduced by working at higher PN orders, see [19] . One then investigates how 
deviations of /3ppN and 7ppn from their GR values affect various experiments. Writing 
/SppN = 1-1-/9 and 7ppn = 1 + 7, the best current limits (at 68% c.l.) are 

7 = (2.1 ±2.3) X 10"^ (4.6) 

from the Doppler tracking of the Cassini spacecraft, and 

4/3-7 = (4.4± 4.5) X 10"^ (4.7) 

from lunar laser ranging. This bound comes from the Nordtvedt effect, i.e. from the 
fact that, in a theory with /3 and 7 generic, the weak equivalence principle is violated 
and the Earth and the Moon can fall toward the Sun with different accelerations, 
which depend on their gravitational self-energy. The effect is studied by monitoring 
the Earth-Moon distance with lunar laser ranging. The perihelion shift of Mercury 
gives instead the bound < 3 x 10"^ [T91IT38] . 

To get these bounds, we do not need to know the fundamental theory that gives rise 
to values of /SppN and 7ppn that differ from their GR values. However, it is of course 
interesting to see that consistent fleld theories exist that give rise to values of /3ppN and 
7ppn different from one. For instance, a Brans-Dicke theory with parameter cjbd gives 
/3ppN = 1 and 7ppn = (1 + i^bd)/(2 -I- wbd), with the GR value 7ppn = 1 recovered 
for wbd 00, while more general tensor-scalar theories can produce both 7ppn 7^ 1 
and /3ppN 7^ 1- However, in the PPN approach, one can also explore other possibilities, 
such as PPN parameters that correspond to preferred-frame effects or to violation of 
the conservation of total momentum, which are not necessarily well-motivated in terms 
of current fleld-theoretical ideas on possible extensions or UV completions of GR. It is 
also important to observe that the parameters /3ppN and 7ppn are gauge-invariant, and 
therefore observables, because they have been deflned with respect to a speciflc gauge, 
namely the standard PPN gauge in which the metric takes the form fl4.ip -f Hy3|) . 

A second, complementary, approach to the problem is to study a speciflc class of 
fleld-theoretical extensions of GR. A typical well-motivated example is provided by 
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multiscalar-tensor theories. These have been studied in detail and compared with 
experimental tests of relativistic gravity in refs. [HH |37l [Ml [HI [33l [121 IS]- This 
approach has the advantage that one is testing a specific and well-defined fundamental 
theory. On the other hand, an experimental bound on the parameters of a given scalar- 
tensor theory, such as for instance the bound cjbd > 40000 on the parameter cjbd of 
Brans-Dicke theory obtained from the tracking of the Cassini spacecraft |32] is, strictly 
speaking, only a statement about that particular extension of GR and not about GR 
itself. 

In this paper we quantify how well GR performs with respect to experiments of 
relativistic gravity by studying how much these experiments constrain the values of 
the non-Abelian vertices of the theory, in particular the three-graviton vertex and the 
four-graviton vertex. We proceed as follows. After choosing a gauge (the De Bonder 
gauge, corresponding to harmonic coordinates) we multiply the three-graviton vertex 
by a factor (1 -|- /^s) and the four-graviton vertex by a factor (1 -|- (34), with constants 
(3^ and /34. For (33 = (3^ = we recover GR. Observe that, since (3^ and (3^ are defined 
with respect to a given gauge choice, they are gauge-invariant by definition. This is 
in fact the same logic used to define in a gauge-invariant manner the PPN parameters 
/3pPN and 7ppn. 

We then use a Feynman diagram approach to compute the modifications induced by 

and on various observables in classical GR. Diagrammatic approaches and field- 
theoretical methods have been in use in classical GR for a long time, see e.g. [T2| 
I139[ 1140] . We make use of the effective field theory formulation proposed in [6], 
which provides a clean and systematic separation of the effects that depend on (model- 
dependent) short-distance physics from long-wavelength gravitational dynamics, and in 
the non-relativistic limit (after performing a multipole expansion) has manifest power 
counting in the typical velocity v of the source. 

In the next section we see how these deformed vertices give additional terms in the 
PN effective Lagrangian. In particular we find that, in the conservative sector of the 
theory, the introduction of (3^ is phenomenologically equivalent, at IPN level, to the 
introduction of a non-trivial value of /3ppN given by /3ppN = 1 + /^s- However, (3^ also 
affects the radiative sector of the theory, i.e. the Lagrangian describing the interaction 
between the matter fields and the gravitons radiated at infinity. 

Before entering into the technical aspects, however, let us clarify the meaning of the 
introduction of P^ and (34. In ordinary GR, with ^3^4 = 0, coordinate transformation 
invariance ensures that the negative norm states decouple. After gauge fixing (in the 
De Donder gauge for instance), the kinetic terms for all of the ten components of 
the metric are invertible, but four of them have the wrong sign, i.e. they give rise 
to negative norm states. In the De Donder gauge the six positive-norm states that 
diagonalize the kinetic term are 

hij = hij + ^Sijihoo - S^'^hi^) , (4.8) 

while the four "wrong-sign" components are given by the spatial vector hoi and by the 
scalar 

h^ = hoo-5'^hi^. (4.9) 
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In standard GR the existence of these negative-norm states do not create difficulties 
because they are coupled to four integrals of motion (energy, and the three components 
of angular momentum), so they cannot be produced. In contrast, the remaining six 
"healthy" components couple to the source multipole moments. After complete gauge 
fixing one finds that among the six positive-norm states, four obey Poisson-like equa- 
tions, so they do not radiate (even though they are non-radiative physical degrees of 
freedom), while the remaining two are the radiative degrees of freedom representing 
GW's [lIT] . 

Allowing /Ss 7^ has the effect that the negative norm state Tin now couples, already 
at lowest order, to a non conserved quantity, namely to a combination of the Newtonian 
kinetic and potential energy of the binary system. This means that, in general, a 
modification of GR in which we just change the strength of the three-graviton vertex 
cannot be taken as a fundamental field theory, neither at the quantum level, nor even 
at the classical level, since the negative-norm state contributes to the classical radiated 
power (a related concern is that, for ^ 0, the energy-momentum tensor is in general 
not conserved). A consistent classical and quantum field theory could in principle 
emerge from a simultaneous modification of all the vertices of the theory, such as the 
three-, four- and higher-order graviton vertices, together with a related modification of 
the graviton-matter couplings. As a trivial example, an overall rescaling of the gauge 
coupling in a Yang-Mills theory, or of Newton's constant in GR, results in a combined 
modification of all the vertices, but obviously introduces no pathology. Anyway, our 
approach to the problem is purely phenomenological. We introduce (5^ and (i^ simply 
as "tags" that allow us to track the contribution of the three- and four-graviton vertices 
throughout the computations. As long as I/Ssi <^ 1 and \[3i\ -C 1, the corrections that 
they induce to the radiated power are small compared to the standard GR result, so 
the total radiated power is given by the GR result plus a small correction, and in 
particular the total radiation emitted is positive. At this phenomenological level the 
introduction of modified vertices is therefore acceptable, and provides a simple and, 
most importantly, gauge invariant manner of quantifying how well different observations 
constrain the non-linear sector of GR, in a way which is intrinsic to GR itself, without 
reference to any other specific field theory. 

In this sense, our approach is close in spirit to the phenomenological PPN approach, 
and can be seen as an extension of it where the radiative sector of the theory is also 
modified. Another approach which is related to ours is the one proposed in ref. [T5|[T6]. 
They consider the phase of the GW emitted during the coalescence of compact binaries, 
which up to 3.5PN has the form 

7 

vl/(/) = 27r/t, - + ^[^, + In , (4.10) 

where / is the GW frequency, and tc and $c are the time and the phase at merger. The 
seven non-zero coefficients ipk with k = 0, 2, 3, ... 7 and the two non-zero coefficients 
ipki with k = 5,6 are known from the PN expansion, in terms of the two masses mi 
and 1712. In ref. [151 IB] they study how the template is affected if these coefficients are 
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allowed to vary, so that two of them, the OPN coefficient ifjo and the IPN coefficient 
"02, are used to fix the masses mi and m2 of the two stars, while varying any of the 
remaining coefficients with k > 3 provides a test of GR. In the case of coalescing 
binaries at interferometers, our introduction of and (3^ is a particular case of a more 
general analysis in which one treats the quantities ipk and ipki as free parameters, but 
it has a sharper field-theoretical meaning since and (3^ measure the deviation of 
the three- and four-graviton vertices from the GR prediction. For the same reason, 
we are also able to compare the effect of on the waveform of coalescing binaries 
with its effect on binary pulsar timing and on solar system experiments, while in the 
phenomenological approach in which the parameters ipk and ipki of the GW phase are 
taken as free parameters, a modification of the waveform of coalescing binaries cannot 
be related to a modification of the binary pulsar timing formula. 

Another issue is whether a modification of the vertices of this form (typically with 
Ps, P4,, etc. not independent, but related to each other by some consistency conditions) 
could emerge from a plausible and consistent extension of GR. Actually, a typical 
UV completion of GR at an energy scale A will rather generate corrections to the 
vertices that are suppressed by inverse powers of A, so it would give rise to an energy 
dependent P^, e.g. P^ = E'^/A'^, which furthermore, at the energy scales that we 
are considering and for any sensible choice of A, would be utterly negligible. Still, 
let us remark that this kind of behavior is not a theorem. It assumes the UV-IR 
decoupling typical of effective field theories, and one can exhibit counterexamples. For 
instance, in non-commutative Yang-Mills theories there is a UV-IR mixing, such that 
low-energy processes receive contributions from loops where very massive particles are 
running, and these contributions are independent of the mass of these particles |142] . 
Anyway, again our aim here is not to test any given consistent extension of GR, but 
rather provide a simple and phenomenologically consistent way of quantifying how well 
various experiments can test the non-linearities of GR, and quantify how the results of 
different experiments compare among themselves. 

It is also interesting to observe that, even when P^, and /34 are non-zero, the graviton 
remains massless at the classical level, since P^ and affect interaction terms, but not 
the kinetic term. The breaking of diffeomorphism invariance induced by P^ and 
could in principle generate a graviton mass at the one-loop level. However, even if we 
are using the language of quantum field theory, in the end we are only interested in the 
classical theory, since quantum loops are suppressed by powers of h/L, where L is the 
angular momentum of the system, so they are completely negligible for a macroscopic 
system^]. 



1. At the quantum level, if a mass is generated, it will be power divergent with the cutoff, and 
will have to be fine tuned order by order in perturbation theory. However, the fact that quantum 
divergences have to be subtracted order by order is a generic problem of the standard quantum 
extension of GR, independently of ^3. In any case, again, our approach is purely phenomenological, 
and the introduction of f]^ is simply a tool for tracking a specific contribution to the computation. 
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4.3 Effective Lagrangian from a modified three-gra- 
viton vertex 

To perform our computations we use tlie effective field theory formulation pro- 
posed in ref. |6]. Computations of the conservative dynamics at 2PN level have been 
performed using this effective field theory technique |143t 1144] and the results are in 
agreement with the classic 2PN results of refs. |52l [53], while the full 3PN result for 
non-spinning particles to date has only been obtained with the standard FN formalism 
using dimensional regularization |139] 1140] (see ref. |107] , or Chapter 5 of ref. [9^ , for 



a pedagogical introduction to the PN expansion and for a more complete list of refer- 
ences). Spin-spin contributions at 3PN order in the conservative two-body dynamics 
have recently been computed both with the effective field theory techniques |145l 1146] 
and with the ADM Hamiltonian formalism |147] (see also |105l 1148] for other applica- 
tions of the EFT technique related to gravitational radiation). 

In the formalism of ref. j6] , after integrating out length-scales shorter than the size 
of the compact objects, the action becomes 

5 = 5eh + 5pp, (4.11) 
where S'eh is the Einstein-Hilbert action and 

Spp = -J2'^- [ (4- 12) 

is the point particle action. Here a = 1, 2 labels the two bodies in the binary system and 
dTa = \J g^i,{xa)dxadx\^. One then observes that, in the binary problem, the gravitons 
appearing in a Feynman diagram can be divided into two classes: the forces between 
the two bodies with relative distance r and relative speed v are mediated by gravitons 
whose momentum scales as {k^ ~ f /r, |k| ~ 1/r). These are called "potential 
gravitons" and are off-shell, so they can only appear in internal lines. The gravitons 
radiated to infinity rather have (A;°~t>/r, |k|~f/r). One then writes g^y = t]^^ + h^^ 
and separates h^^^, into two parts, h^^{x) = h^u{x) + H^^{x) with h^^{x) describing 
the radiation gravitons and H^y{x) the potential gravitons. One fixes the de Bonder 
gauge and, expanding the action in powers of /i^j^ and H\f_fj_t,{xo), one can read off the 
propagators and the vertices, and write down the Feynman rules of the theory. Then, 
using standard methods from quantum field theory, one can construct an effective 
Lagrangian that, used at tree level, reproduces the amplitudes computed with the 
Feynman graphs. For a classical system, whose angular momentum L ^ h, only tree 
graphs contribute, and reproduce the classical Lagrangian that is usually derived from 
GR using the PN expansion. For instance, the conservative dynamics of the two-body 
problem, at IPN level, is given by the Einstein-Infeld-Hoffmann Lagrangian, 



1 4 1 4 GTvmims 
^EiH = -rriiv^ + -m2V2 + 



-7vrV2 - (r-Vi)(f-V2) - 



2r 



?,{vl + vl) 

(4.13) 
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(1 + /?3)a 



(a) (b) 

Figure 4.1: The diagrams that give the terms proportional to Gj^ in the IPN Lag 
rangian. Dashed hnes denote potential gravitons, solid lines the point-like sources. 



In the language of the effective field theory of ref. [6], this result is obtained from Feyn- 
man diagrams involving the exchange of potential gravitons. In particular, the terms 
linear in Gn in eq. fl4.13p are obtained from a single exchange of potential gravitons 
between two matter lines, see Fig. 4 of ref. [6], while the term proportional to Gj^ is 
obtained from the sum of the two graphs shown in Fig. 14.11 In the derivation of the 
conservative IPN Lagrangian, the three-graviton vertex only enters through the graph 
in Fig. 14.1b . Multiplying this vertex by a factor (1 + P^) and repeating the same com- 
putation as in ref. [6] we get the additional contribution to the conservative part of the 
Lagrangian 

A£.„„.^-fe ^"""";'r' + (4.14) 

Comparing this result with the Lagrangian whose equations of motion are the same 
as the equations of motion of a test particle in the PPN metric fl4.ip - fl4.3p (see eq. (6.80) 
of ref. [27]) we find that, to IPN order and as far as the conservative dynamics is 
concerned, the introduction of /Ss gives rise to a PPN theory with P = 1 + P^, i.e. 
P = P^, while 7 = 1 as in GR. Therefore the bound on P from the perihelion of 
Mercury translates into \P^\ < 3 ■ 10~^, while eq. (14. 7p translates into the bound (at 
68% c.l.) 

I/33I <2-10-^ (4.15) 

This bound reflects the fact that /Ss, just as the PPN parameter /3, violates the weak 
equivalence principle. The introduction of P^, however, also affects the radiative sector 
of the theory, something that is not modeled in the phenomenological PPN framework 
since the latter by definition is only concerned with the motion of test masses in a 
deformed metric, and therefore only modifies the conservative part of the dynamics. 
Note however that, in the framework of multiscalar-tensor theories, the extension of 
the PPN formalism introduced in Ref [llj allows for a consistent treatment of both 
the conservative dynamics (including the case of strongly self-gravitating bodies) and 
of radiative effects. 

It is clearly interesting to see what bounds on P^ can be obtained from experiments 
that probe the radiative sector of GR, such as the timing of binary pulsars or the 
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(a) 



(b) 



(c) 



Figure 4.2: The diagrams that contribute to the matter-radiation Lagrangian. Dashed 
hues denote potential gravitons, wiggly hues radiation gravitons, and sohd hues the 
point-hke sources. 



observation of the coalescence of compact binaries at interferometers. The effective 
Lagrangian describing the interaction of the binary system with radiation gravitons is 
obtained by computing the three graphs in Fig. 14.21 (corresponding to Fig. 6 of ref. [6]), 
and the introduction of affects the HHh vertex in Fig. 14.2b . 

Computing these graphs as in ref. [6], but with our modified three-graviton vertex, 
we find 

Aad = [QijRoiOj + qRom + Psi^Vhoo + Z'^hij)] , (4.16) 

ZIVlpi 

where Qij is the quadrupole moment of the source and we define 



(4.17) 



r 

Z-(r) = ^^^!^^i^, (4.19) 

where r = xi — X2. The term QijRoioj in eq. fl4.16p is the usual quadrupole interaction. 
The second term, qRoioi, is non-radiating when (3^ = 0, but we will see that for /Sa 7^ 
it contributes to the radiated power when the orbit is non-circular. The last two terms 
in eq. f l4.16p are the explicit /^s-dependent terms induced by the modification of the 
three-graviton vertex in Fig. 14.2b . 

Finally, we have omitted a term where /zqo is coupled to the conserved energy (at 
this order) and hoi is coupled to the conserved angular momentum, since these terms 
do not generate gravitational radiation. 

The back-reaction of GW emission on the source can be computed as usual from 
the energy balance equation 

Pgw = -E, (4.20) 

where Pgw is the power radiated in GWs and E the orbital energy of the system. 
To obtain the expression for the radiated power for /^s 7^ we cannot simply use the 
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quadrupole formula of GR, since the introduction of /S^ generates new contributions. 
To take them into account we proceed as in ref . [6] , and we compute the imaginary part 
of the graph shown in Fig. 14.31 The vertices of the graph can be read from eq. (14.16^ . 
When /Ss = the only relevant vertex comes from QijRoioj, and the computation of 
the imaginary part gives back the Einstein quadrupole formula |6] . In our case we have 
various possible vertices, and we must compute the imaginary part of 

^ / dtidt2/i;(ti)/L(i2)(5,';.(ti)5^,(t2)) , (4.21) 

®^^^P1 a,h=l -J 

where If^ = {Qij,q6ij, P^VSij, P^Zij) depends on the matter variables and S^j = {Roioj, 
RoiOjySijhoo,hij) on the gravitational field. When both vertices of the diagram in 
Fig. 14.31 are proportional to the quadrupole, one obtains the usual GR result 

Pqq = ^(Q^A,)^ (4-22) 

as already found in |6j. Computing the other contributions we find that the terms 
Pqq and Pgg vanish identically. In fact, the Qq and qq graphs vanish because Qij 
is traceless, while the qq graph vanishes because 6ij6ki gives zero when contracted 
^ikSji + Sii6jk — ^Sijdki, which is the tensor that comes out from the two-point function 
(RoiOjRokoi) ■ The QV, qZ and VZ graphs vanish for similar reasons, so the only 
relevant contributions come from the QZ and qV graphs, and we find 

Pqz = -2(3sGn(Q.jZ,,), (4.23) 

and 

PgV = -6/33Gr,(qV). (4.24) 

As for the VV and ZZ graphs, they give a contribution that, from the point of view 
of the multipole expansion, is of the same order as the quadrupole radiation but pro- 
portional to and can be neglected. 

We can now use these results to perform the comparison with binary pulsars and 
with interferometers. 

4.4 Comparison with experiments 

As we already saw in eq. fl4.15p . solar system experiments, and in particular lunar 
laser ranging, give the bound l/^al < 2 ■ 10""^. In this section we study the bounds on /^s 
that can be obtained from binary pulsars and from the detection of coalescing binaries 
at interferometers. 

4.4.1 Binary pulsars 

Since (3^ modifies the emitted power already at Newtonian order, the energy of the 
orbit in eq. (14.201) can now be directly computed using the Keplerian equations of mo- 
tion. We see that this test of GR is conceptually different from the tests based on solar 
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Figure 4.3: The self-energy diagram whose imaginary part gives the radiated power. 
The wiggly hne can refer either to hij or to hoo, and the vertices to any of the four 
Jf S'f- with a = 1, . . . , 4, see the text. 



system experiments. The latter only probe the conservative part of the Lagrangian, 
i.e. the /Sa-dependent term given in eq. (14.141) . while binary pulsars are sensitive to the 

dependence given in the radiation Lagrangian fl4.16p (even if the effect of on the 
conservative dynamics will also enter, through the determination of the masses of the 
stars from the periastron shift, see below). 

Using the Keplerian equations of motion for an elliptic orbit of eccentricity e we get 

32G%n'^M^ ( 73 2 37 4 



P,, . hl^fM^ il ^ + 'lA . (4.26) 



5a5(l - e2)7/2 24 96 



where M = mi + m2 is the total mass, fi = 17111712/ M is the reduced mass, and we will 
also use the notation u = 77117712/ M'^ for the symmetric mass ratio. From the energy 
balance equation we then get the evolution of the orbital period Pb, 

g = -^Gf (g)"'^' [/(e) + PMe)] , (4.28) 

where 

«w - jrz^ {I ^ f »^ - w/) ■ («o) 

The term proportional to /(e) is the standard GR result |149] . while the term propor- 
tional to g{e) is the extra contribution due to f3s. 

In order to measure from the dynamics of binary pulsars, we must also determine 
the dependence on P^ of the periastron shift u and of the Einstein time delay 7s, since 
these two observables are used to determine the masses of the two compact stars. In 
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particular the periastron shift fixes the total mass M of the system, while the Einstein 
time delay 7^; measures a different combination of masses, see e.g. eqs. (6.56) and 
(6.93) of ref. j90j. Using the conservative Lagrangian with the modification (I4.14p and 
repeating the standard textbook computation of the periastron shift cj, we find that 
the value up.^ computed in a theory with /^s 7^ is related to the GR value cjqr by 

while the Einstein time delay is unchanged because it is not affected by the post- 
Keplerian parameters. So, if /Ss 7^ 0, the true value of the total mass M of the binary 
system, that enters in eq. fl4.28p . is not the one that would be inferred from the peri- 
astron shift using the predictions of GR, but rather we get 

Mp, = (1 + y) ^GR . (4.32) 

Similarly, using eq. (6.93) of ref. [90], for the symmetric mass ratio v we get 

vp, = (l + w/33)z/GR, (4.33) 

where 

^ + 1 

3 yryy^ {i + aky/^-{i + k) ^ > 

and 

/o \ 1/3 

K = — 



(^) (G^vMgr)-^/^ . (4.35) 



Putting everything together and keeping only the linear order in /^s we finally find that 
the ratio between the value of Ph computed at /^a 7^ and the value of Ph computed in 
GR is 



p(/33) 

j_b 

pGR 



l + /33^(e), (4.36) 



where 



/(e) 6 

Observe that the term g{e)/f{e) comes from the effect of P3 on the radiative sector of 
the theory, while the term {5/6) + w comes from the effect of ^3 on the conservative 
sector, i.e. on the mass determination. Inserting the numerical values for the Hulse- 
Taylor pulsar, we get g{e) ~ 3.21. Note that g{e)/ f{e) ~ 2.38, so g{e) is dominated 
by the effect of on the radiative sector of the theory. For this binary pulsar, after 
correcting for the Doppler shift due to the relative velocity between us and the pulsar 
induced by the differential rotation of the Galaxy, the ratio between the observed value 
Pb°^" and the GR prediction P^^ is P^^^yP^^ = 1.0013(21). Interpreting this as a 
measurement of /Ss we finally get 3.21/33 = 0.0013(21), i.e. 

/33 = (4.0 ± 6.4) ■ 10-\ (4.38) 
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so the three-graviton vertex is consistent with the GR prediction at the 0.1% level. 
This bound is slightly worse, but comparable, to the one from lunar laser ranging, 
eq. (I4.15p . It should be stressed, however, that eq. ( 14.38^ is really a test involving 
the radiative sector of GR, while eq. (14.15^ only tests the conservative sector. For 
comparison, observe that in the Standard Model the triple gauge boson couplings are 
measured to an accuracy of about 3% |137j . 

For the double pulsar we find g{e) ~ 3.3. Since Pb for the double pulsar is presently 
measured at the 1.4% level |129] . we get a larger bound compared to eq. ( 14.38^ . How- 
ever, further monitoring of this system is expected to bring the error on Pb down to 
the 0.1% level. 

4.4.2 Binary coalescences at interferometers 

We now compare these results with what can be expected from the detection of a 
binary coalescence at GW interferometers. In this case one can determine the physical 
parameters of the inspiraling bodies, by performing matched filtering of theoretical 
waveform templates. In the matched filtering method any difference in the time beha- 
vior between the actual signal and the theoretical template model will eventually cause 
the two to go out of phase, with a consequent drop in the signal-to-noise ratio (SNR). 
The introduction of (3^ and /34 affects the template, in particular the accumulated phase 



where f{t) is the time- varying frequency of the source, and the subscript min (max) 
denotes the values when the signal enters (leaves) the detector band-width. Thus in 
principle a detection of a GW signal from coalescing binaries could be translated into 
a measurement of the three- and four-graviton vertices. In this section we investigate 
the accuracy of such a determination. 

With respect to the timing of binary pulsars, there are at least three important 
qualitative differences that affect the accuracy at which these systems can test the non- 
linearities of GR. First, coalescing compact binaries in the last stage of the coalescence 
reach values of f/c ~ 1/3, and are therefore much more relativistic than binary pulsars, 
which rather have f/c ~ 10~^. Second, the leading Newtonian result for (p is of order 
{v/c)~^ so it is much larger than one, and to get the phase with a precision A0 <C 1, as 
needed by interferometers, all the corrections at least up to 0{v^/c^) to the Newtonian 
result must be included, so higher-order corrections are important even if they are 
numerically small relative to the leading term. In other words, even if PN corrections 
are suppressed by powers of f/c with respect to the leading term, they can be probed 
up to high order because what matters for GW interferometers is the overall value of 
the PN corrections to the phase, and not their value relative to the large Newtonian 
term. These two considerations should suggest that interferometers are much more 
sensitive than pulsar timing to the non-linearities of GR. 

On the other hand, for binary pulsars we can measure not only the decay of the 
orbital period due to GW emission, but also several other Keplerian observables, that 




(4.39) 
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provide a determination of the geometry of the orbit, as well as post-Keplerian obser- 
vables, such as the periastron shift and the Einstein time delay, which fix the masses 
of the stars in the binary system. This is not the case for the detection of coalescences 
at interferometers. With interferometers the parameters that determine the waveform, 
such as the masses and spins of the stars, must be determined from the phase of 
the GW itself, and one must then carefully investigate the degeneracies between the 
determination of (3^ (or of (34) and the determination of the masses and spins of the 
stars. This effect clearly goes in the direction of degrading the accuracy of parameter 
reconstruction at GW interferometers, with respect to binary pulsar timing, so in the 
end it is not obvious a priori which of the two, GW interferometers or binary pulsar 
timing, is more sensitive to the non-linearities of GR. This question is answered in what 
follows. 

Repeating with ^ the standard computation of the orbital phase for a cir- 
cularized orbit, we find that /Ss modifies the orbital phase (j)(t) already at OPN (i.e. 
Newtonian) level, where to linear order in we get 

95/8 

0°™ = -^(l + &o/33), (4.40) 
with 60 = —5/2, and is defined as 

where tc is the time of coalescence. Combining the factors u and M which enter in the 
definition of O with the explicit factor 1/z/ in eq. f l4.40p we recover the well-know result 
that the Newtonian phase depends on the masses of the stars only through the chirp 
mass Mc = z/^/^M. From eq. (I4.40p we immediately understand the crucial role that 
degeneracies have for interferometers. In fact, since Mc is determined from eq. (I4.40p 
itself, using only the OPN phase fl4.40p it is impossible to detect the deviation from the 
prediction of GR induced hj (S^. A non-zero value of would simply induce an error 
in the determination of Mc- 

The same happens at IPN level. In fact, at IPN order and with Ps ^ the phase 
has the general form 

^iPN ^ ^ ^^^^^ ^ ^ hPs)e-'/'] , (4.42) 

where, as before, Bq = —5/2 is the OPN correction proportional to P^, while 

, , 3715 55 

ai(iy) = + —U 4.43 

^ ' 8064 96 ^ ^ 

is the IPN GR prediction |107[ [90]. and bi (which is possibly i/-dependent) parametrizes 
the IPN correction due to P3. (For simplicity, we only wrote explicitly the term linear 
in P3 since \Ps\ is much smaller than one, but all our considerations below can be 
trivially generalized to terms quadratic in P3, just by allowing the function &i(z^) to 
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depend also on /^s). In general we expect bi to be 0(1), and we will see below that for 
our purposes this estimate is sufficient. 

Using Mc and u as independent mass variables, in place of mi and m2, we see that 
while the effect of on the OPN phase can be reabsorbed into Mc, its effect on the IPN 
phase can be reabsorbed into a rescaling of u. Observe that, in the detection of a single 
coalescence event, GW interferometers do not measure the functional dependence of u 
of the IPN phase, but only its numerical value for the actual value of u of that binary 
system, so we cannot infer the presence of a term proportional to (3^ from the fact that 
it changes the functional form of the i/-dependence from the one obtained by eq. ( 14.43^ . 
Thus, even at IPN order, it is impossible to detect the deviations from GR induced by 
P^. A non-zero (3^ would simply induce an error on the determination of and z/, i.e. 
on the masses of the two stars. 

We then examine the situation at 1.5PN order. Let us at ffist neglect the spin of 
the two stars. Then the 1.5PN phase with ^ has the generic form 



Q5/8 

= [(1 + feo/Ss) + ai(i/)(l + 61/33)6-1/^ 

+ 02(1 + 62/33)6-3/^ 



(4.44) 



where 02 = —Sir /A is the 1.5PN GR prediction and 62 is the (possibly //-dependent) 
1.5PN correction due to P^. Again, we will not need its exact value, and we will simply 
make the natural assumption that it is 0(1). 

However, for the purpose of determining /^s, neglecting spin is not correct. Indeed, 
at 1.5PN order the spin of the bodies enters through the spin-orbit coupling, and the 
evolution of the GW frequency / with time is given by |150] 



24 

1-— ai(z/)x + (47r-/3Ls)a;3/^ 
5 



(4.45) 

where x = (vrM/)^/^, while /3ls describes the spin-orbit coupling and is given by 



1 ^ ' 



a=l 



2 

Tfl 

113-^ + 75. 



L-X. , (4,46) 



where L is the orbital angular momentum, Xa — ^a/^'i and is the spin of the a-th 
body. In principle /3ls evolves with time because of the precession of L, Si and S2. 
However, it turns out that in practice it is almost conserved, and can be treated as 
a constant [151j . Integration of eq. (14.45^ then shows that, in the 1.5PN phase, /3ls 
is exactly degenerate with [3^ in eq. fl4.44p . Furthermore, observe that /Slsj depending 
on the spin configuration, can reach a maximum value of about 8.5 |151j (and its 
maximum value remains large even in the limit z/ — )■ 0), while [i^ is already bound by 
laser ranging at the level of 2 x 10"^ and by pulsar timing at the level of 10^'^ (which 
tests the radiative sector, as do GW interferometers). Thus, the effect of (3^ at 1.5PN 
is simply reabsorbed into a (very small) shift of /3ls- 
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(1 + P4j>: 
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(a) (b) 

Figure 4.4: The diagrams that contribute to the conservative dynamics, which are 
affected by a modification of the four-graviton vertex. 



At 2PN order fij, is degenerate with the parameter a that describes the spin-spin 
interaction 

a = ^[721(L-Xi)- 247Xi-X2(L-X2)] 

(4.47) 

" ~ a=l 

The first term is the one which is usually quoted in the literature, first computed in 
|150j (see also |152[ 1153] ). The term in the second line, computed recently in [154] . is 
however of the same order, and must be included. 

The first term is proportional to and reaches a maximum value o^axiy) — lOi^. 
In a coalescence with very small value of this term is therefore suppressed; e.g. in an 
extreme mass-ratio inspiral (EMRI) event at LISA [ID] where a BH of mass mi = IOMq 
falls into a supermassive BH with = IO^Mq, one has u = 10^^ and the term in the 
first line has a maximum value ~ 10~^. If this standard term gave the full answer, a 
value of in excess of this value could therefore give an effect that cannot be ascribed 
to a. However, the presence of the new term recently computed in [154j spoils this 
reasoning, since it is not proportional to z/. The conclusion is that, just as with /3sl at 
1.5PN order, the effect of (3^ at 2PN order is just reabsorbed into a small redefiniton of 
0", and therefore simply induces an error in the reconstruction of the spin configuration 
(observe also that fixing /3ls does not allow us to fix the spin combinations that appear 
in ex.) 

One could in principle investigate the effect of on higher-order coefficients, such 
as the 2.5PN term ip^ and the 3PN term i/jq in eq. (I4.10p . which at LISA can be 
measured with a precision of order 10~^ [16]. Unfortunately, it is difficult to translate 
them into bounds on because at 2.5PN order one finds that P^ is degenerate with 
a different combination of spin and orbital variables, which is not fixed by the 1.5PN 
spin-orbit term (see Table II of ref. |155] ). while at 3PN order the spin contribution is 
not yet known. The conclusion is therefore that interferometers cannot measure the 
three- and higher-order graviton vertex, since the effect of a modified vertex is simply 
reabsorbed into the determination of the masses and spin of the binary system. The 
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Figure 4.5: The diagram that contributes to the radiative dynamics, which is affected 
by a modification of the four-graviton vertex. 



conclusion that interferometers are not competitive with pulsar timing for measuring 
deviations from GR was also reached in ref. [36], although in a different context. In fact, 
ref. [36] was concerned with multiscalar-tensor theories, whose leading-order effect is 
the introduction of a term corresponding to dipole radiation (a "minus one"-PN term). 

We now examine what can be said about the four-graviton vertex, parametrized by 
/34. In the conservative part of the Lagrangian contributes through the diagrams 
of Fig. 14. 4[ However, these contributions only affect the equations of motion at 2PN 
order, so there is no hope to see them in solar system experiments, where the velocities 
at play are very small. For the same reason, no significant bound can be obtained from 
binary pulsars; from a simple order of magnitude estimate we find that the Hulse- Taylor 
pulsar can only give a limit /34 < 0(10), which is not significant. 

At GW interferometers, enters into the phase for the first time at IPN order, 
through the diagram in Fig. 14.51 However, it suffers exactly of the same degeneracy 
issues as /^s, so it cannot be measured to any interesting accuracy at present or future 
interferometers, at least with the technique discussed here. Note however that in 
this paper we have worked in the restricted PN approximation, in which only the 
harmonic at twice the source frequency is retained. Higher-order harmonics however 
break degeneracies between various parameters in the template |156j . and it would be 
interesting to investigate whether their inclusion in the analysis allows one to put a 
bound on /^a and fi^ from binary coalescences. 

4.5 Conclusions 

We have proposed to quantify the accuracy by which various experiments probe 
the non-linearities of GR, by translating their results into measurements of the non- 
Abelian vertices of the theory, such as the three-graviton vertex and the four-graviton 
vertex. This is similar in spirit to tests of the Standard Model of particle physics, where 
the non-Abelian vertices involving three and four gauge bosons have been measured at 
LEP and at the Tevatron. 

We have shown that, at a phenomenological level, this can be done in a consistent 
and gauge-invariant manner, by introducing parameters /Ss and /34 that quantify the 
deviations from the GR prediction of the three- and four-graviton vertices, respectively. 
We have found that, in the conservative sector of the theory, i.e. as long as one neglects 
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the emission of gravitational radiation at infinity, the introduction of /3s at IPN order is 
phenomenologicaUy equivalent to the introduction of a parameter /SppN = 1 + /Ss in the 
parametrized PN formalism. Strong bounds on /^s therefore come from solar system 
experiments, and most notably from lunar laser ranging, that provides a measurement 
at the 0.02% level. 

The modification of the three-graviton vertex however also affects the radiative 
sector of the theory, and we have found that the timing of the Hulse-Taylor pulsar 
gives a bound on fS^ at the 0.1% level, not far from the one obtained from lunar laser 
ranging. Conceptually, however, the two bounds have different meanings, since lunar 
laser ranging only probes the conservative sector of the theory, while pulsar timing is 
also sensitive to the radiative sector. 

We have then studied the results that could be obtained from the detection of 
coalescences at interferometers, and we have found that, even if already modifies 
the GW phase at the Newtonian level and at IPN order, their effect can always be 
reabsorbed into other parameters in the template, such as the mass and spin of the two 
bodies so, rather than detecting a deviation from the GR prediction, one would simply 
make a small error in the estimation of these parameters. 
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Chapter 5 



The effects of modified graviton 
vertices in the templates 
for binary coalescences 

5.1 The possibility of measuring the non-Abehan 
vertices of General Relativity with interfero- 
meters 

In the previous chapter I have described a parametrized framework to constrain 
gravity non-hnearities from an effective field theory point of view [I25j. In this work, 
performed with the rest of my group, we tagged the three- and four-graviton vertices 
of GR with some parameters /Sa and /94, respectively, and we investigated the bounds 
that could be put on them by means of experiments of relativistic gravity. One of the 
conclusion of our study is that the timing of binary pulsars can provide better tests 
of GR than what is possible with GW observations because the latter are plagued by 
degeneracies among parameters. A similar conclusion in favor of binary pulsar tests has 
been derived by DEF in their study of scalar-tensor theories, which I have discussed 
in Sec. 11.11 As pointed out in that section, the validity of this type of conclusions is 
limited to the PN orders probed by pulsar timing, i.e. up to terms of (9(f/c)^: these 
comprise only the lowest orders of the radiative regime, while GW observations aim 
at probing the PN structure at least until the 3.5 next-to-leading order order, i.e. up 
to terms of 0{v/cY'^. Moreover, in confronting the probing power of different testing 
grounds, one should also consider how the extraction of physical parameters is done in 
the various experiments. Through the timing of binary pulsars one can measure the 
decay of the orbital period, which is a consequence of GW emission; on top of the decay, 
one has access to other post-Keplerian parameters, like the periastron shift and Einstein 
time delay, which are used to fix the unknown masses of the stars in a binary system. 
The situation is radically different for what concerns the detection of coalescences 
at interferometers: the entire physical information on the emitting system has to be 
extracted from the phase of the GW signal that has spent many cycles in the detector. 
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5.2 Constraining alternative theories with spin effects in the template 



Therefore, every type of parameters, be them the constituent masses or the coupling 
constants of alternative theories of gravity, must be determined from the a unique 
measured quantity. These differences between binary pulsar timing and coalescences 
at interferometers are at the heart of the conclusions of the work with my group [125J. 
In fact, on one hand we were able to constrain the three- and four-graviton vertices 
to interesting accuracies using the timing of Hulse- Taylor pulsar; on the other, by 
looking at the GW phasing formula we could not derive any bound: possible deviations 
of the non-Abelian vertices from their GR values were found to be degenerate with 
numerically larger phasing coefficients of the GW template, notably those due to the 
spins of the objects. Because we worked in the restricted PN framework, we concluded 
that these degeneracies could be broken by the inclusion of amplitude corrections in the 
template, as higher harmonics enrich the functional dependence of the parameters |156j . 
Improvements in parameter estimation brought in by the use of full waveforms have 
been discussed in Sec. II. 4[ based on the study of ref. [17]. However, a richer functional 
dependence could offer the possibility of bounding small parameters like /Ss and /34 
even in the restricted PN framework. An example along this direction is the work of 
ref. |132] . where Berti, Buonanno and Will (BBW) investigated the bounds on some 
alternative theories that could be put with LISA [40] observations if templates include 
spin-orbit and spin-spin effects. I will now describe the analysis of BBW by focusing 
on the aspects which are more salient to the problem of estimating (3^ and (3 4^. In doing 
so, I will partially follow the presentation given in ref. |132] . 

5.2 Constraining alternative theories at interfero- 
meters with spin effects in the template - 
The case of the space-borne experiment LISA 

The focus of BBW was on Brans-Dicke (BD) scalar-tensor theory and on the phe- 
nomenological effect of a massive graviton. In the course of this thesis we have already 
encountered the extra contributions due to these alternatives to GR: for a scalar-tensor 
theory the relevant term is the dipolar flux (11.15^ . for a massive graviton it is the shift 
in the IPN phasing coefficient f ll.62p . Within the restricted post-Newtonian approxim- 
ation, BBW used the following expression for the waveform in the frequency domain^ 

M/) = ^^r'/'e^^(/) (5.1) 
1 W'^ 

A -- 



/307r2/3 

where is the luminosity distance given by 



Dl = ^^ Tn- 5.2 

Jo [fiM(l + ^')' + ^A] ^ 



1. Other assumptions behind eq. (|5.ip will be discussed in what follows 
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in which one has assumed a zero-spatial-curvature Universe with Qm the energy density 
of matter and the one of a cosmo logical constant. 

Concerning the phase ip in eq. (15. ip . BBW adopted an expression valid up to order 
2PN: this is indeed the accuracy with which the spin terms of the phasing were known. 
With a slight change of notation with respect to the 2PN phasing formula 

reads 



3 \l{l + z) ■ \^ 756 ' 9 

-levr (7rM/)+ 4/3^5 (vrM/) 



/15293365 27145 3085 A ^ ..r^An , ,..x4/0 /.ox 

+ + V + z/2 (nMf)^/^ -10 a iuMfY'^ \ 5.3 

V 508032 504 72 y ^ ^ \ J ) ^ \ ) 



where: 



The first two terms are related to the time tc and phase $c of coalescence, so they 
basically establish where the waveform begins or ends; M is the total mass of the 
binary, v its mass-ratio, / the frequency variable. Because LISA will be able to 
probe sources at cosmological distances, M includes a redshift factor, i.e. M is 
the observed mass and it is related to the mass Ms measured in the rest frame of 
the source by M = M^(l + z). 

The combination (irMf) sets the PN order of the expansion as f = (ttM/)^/'^. 

The first term inside the braces, the "1", together with the prefactor of the expression, 
comes from the GR quadrupole. 

The second term is the contribution of dipole gravitational radiation in BD theory, 
where S is proportional to the difference of the scalar charges of the bodies. 
From the discussion of Sec. 11.11 we remember that a dipole term is larger than 
the quadrupole by a factor t>~^; however, S can be at most 1, while ubd is at least 
4 X 10^ (as we know from Solar-System tests p2]): therefore, the dipole phasing 
coefficient is numerically small. 

The third term is the IPN phasing coefficient that phenomenologically parametrizes 
a massive graviton by means of the graviton Compton wavelength A^. I have 
discussed it in Sec. II. 4[ where the distance D was given for small redshift only. 
The full expression of D reads 

^ = —E— 5.4 

Ho Jo {l + z'y[QM{l + z'y + QAf^ 



2. Here I use (ttM/) as an expansion parameter for consistency with Sec. 11.41 BBW adopted the 
combination (irMf) where the chirp mass is used instead of the total mass M: this choice brings a 
dependence of the phasing coefficients on the mass ratio v which is not present in Sec. 11.41 
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so it is not a conventional cosmological distance measure as the luminosity dis- 
tance (15. 2p . It emerges when studying how a massive graviton propagates in a 
cosmological background [7^ . For the Hubble constant one can assume Hq = 72 
km s~^ Mpc~^, according to the present observational estimates; however, as we 
will see, the analysis of BBW reveals that the choice of a specific cosmological 
model affects the bound on Xg only weakly. 

From the discussion of Sec. 11.41 we remember that the effect of a massive graviton 
is to alter the arrival time of waves of a given frequency. Moreover, let us men- 
tion also the graviton Compton wavelength is already bound from Solar-System 
tests [157] : Ag> 10^^km~ 3 x 10"^ pc, which corresponds to a bound on the 
graviton mass mg<2 x lO^^^eV (cfr. the discussion of eq. ( 15.17P later on). 

All the remaining terms are those predicted by GR. Notably, the quantities /3ls and a 
represent the spin-orbit and spin-spin contributions to the phasing, respectivelyEl: 



12 

i=l 



2 

L-Si, (5.5) 



2 



113^ + 75z. 



V 

a = — 
48 



XiX2 (-247Si • S2 + 721L • SiL • S2) , (5.6) 



where Sj and L are unit vectors in the direction of the spins and of the orbital 
angular momentum, respectively, and Sj = Xi^i Sj. For black holes, the dimen- 
sionless spin parameters Xi must be smaller than unity, while for neutron stars 
they are generally much smaller than unity. It follows that |/3ls| ^ 9.4 and 
\a\ < 2.5, which makes the phasing coefficients due to spin effects numerically 
large. 

As one can see, the phasing formula ( 15. 3p contains two different types of parameters: 
those which are intrinsic to the source, like the masses, the distance D and the spins 
of the objects, and those which are extrinsic, like the parameters wbd and which 
are characteristic of the theory of gravity. Denoting the full set of parameters by {Qa\i 
one can describe a waveform in the time domain as hi = hi[t; {6a}] where the index i 
labels the waveform and t is the time variable. To detect a GW signal and estimate 
the parameters of the emitting source, it is convenient to adopt the matched filtering 
technique, i.e. to maximize the convolution of a signal with a putative waveform; this 
convolution is expressed by the inner product 

(/^i|M = 2 -±J^-^dj, (5.7) 

where /?-(/) = /i(t) exp{2TTift)dt is the Fourier transform of the GW waveform, 
h*{f) its complex conjugate and S'„(/) is the noise spectral density of the detector. 



3. The expression for the spin-spin contribution differs from the one reported in eq. (|4.47p of 
Chap, m because BBW was pubhshed in 2005 while the extra term in eq. (|4.47p has been calculated 
in 2009 [154] . 
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a quantity that describes the sensitivity of an interferometer as a function of the fre- 
quency. By means of the inner product (15.70 . one can define the signal-to- noise ratio 
(SNR) for a given hi by 

p[h,] = {K\h,y/^ (5.8) 

and the so-called Fisher matrix Tab 



fdh. 


dhi\ 


[dea 





(5.9) 



which can be used to estimate the root mean square error A6a in measuring the para- 
meter 6a- In fact, in the limit of large SNR, these errors can be obtained by taking the 
square root of the diagonal elements of the inverse of the Fisher matrix, 

Aea = ^/^a, S = r-i; (5.10) 

analogously, the coefficients quantifying the correlation between two parameters da and 
9b are given by 

Cab = ^ab/ V^aa^bb ■ (5.11) 

By making use of the Fisher matrix, BBW derive the errors and correlations coef- 
ficients on the various parameters to assess the role that spin terms have in estimating 
parameters with detections at LISA. The candidate events considered by BBW are 
from binary systems in circular orbits with spins aligned or anti-aligned with the or- 
bital angular momentum, i.e. the spins of the objects are non-precessing. The choice 
of neglecting eccentricity is to simplify the treatment. Other simplifications stem from 
modeling LISA as a single interferometer and not taking into account the orbital mo- 
tion of the spacecrafts. The latter assumptions are justified by the fact that the angles 
describing the position of the source and the orientation of LISA are parameters that 
modify the amplitude of the signal and not its phase: for this reason, there is no cor- 
relation among the angles and the intrinsic or extrinsic parameters of the phase (15.31) . 
In the part of the work where they ignore LISA's orbital motion, BBW find that the 
SNR takes the form 

where angular braces indicate the average over LISA's angular response and where the 
limits of integration in the frequency are set as follows. The upper limit is chosen to be 
/fin = niin(/isco5 /end)- Here /isco is twice the conventional (Schwarzschild) frequency 
of the innermost stable circular orbit for a point mass, namely /isco = (6^/^7rM)~^. 
On the other hand, /end = 1 Hz is a conventional upper cutoff on the LISA noise curve. 
The initial frequency /in is determined by assuming that we observe the inspiral over 
a time Tobs before the ISCO and by selecting a cutoff frequency below which the LISA 
noise curve is not well characterized. BBW's default cutoff is /low = 10~^ Hz. The 
choice of these limits of integration also applies to the frequency integrals that are 
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needed to compute the Fisher matrix (see eqs. (15. 9p and (15. 7p ). The expression of the 
SNR (I5.12P will be useful to understand the bounds on the parameters. 

When placing bounds on alternative theories of gravity, BBW treated only spin- 
orbit terms. The reason is due to the dimensionality of the Fisher matrix, which gets 
too large in presence of two types of spin contributions at the same time as gravity 
modifications. Concerning the role of spin terms, the main conclusions of the work of 
BBW are the following ones. 



How the presence of spin terms in the template generally affects parameter 
estimation 

There is a degradation in parameter estimation caused by spin contributions even 
in the absence of any modification to GR. This has been known since the work of 
ref. |152j and is due to the fact that if the spins are non-precessing, they are not 
modulated by the orbital motion: as a consequence, the spins are correlated with other 
intrinsic parameters, like the masses and distances. Indeed, ref. |158] has shown that 
spin modulation allows a better estimation of masses and distances with respect to 
the case when spins are aligned or anti-aligned with the orbital angular momentum. 
Therefore, as a general rule, the presence of a higher number of parameters dilutes the 
available information. 



Constraining Brans-Dicke theory 

To constrain BD theory the best sources are asymmetric systems, which maximize 
the difference of the scalar charges of the bodies in the parameter S of eq. (15. 3p . In fact, 
black holes have vanishing coupling to the scalar field because of the no-hair theorem, 
so black hole binaries are not useful to bound wbd- Neutron star binaries are not 
optimal either because their scalar charge is only weakly dependent on the equation 
of state (cfr. the discussion of Fig. 11.81 in Sec. II. ip . Finally, binaries comprising white 
dwarves can be excluded from consideration because these objects are not so compact, 
so they introduce tidal effects that complicate the analysis. BBW then focus their 
treatment on the inspiral of a neutron star into a black hole with the following masses: 
Mns = 1-4M0 and Mbh ■ (4x 10^- 1 x 10^^)71^0. They set an SNR=10, which, according 
to eq. (I5.12p . implies that a higher total mass M corresponds to a larger luminosity 
distance Dl probed by LISA: notably, one has ~ 3 x 10^ Mpc for M ~ 4 x IO^Mq . 
However, a high total mass worsens the bound on cjbd ioT the following reason. The 
derivative of the waveform (15. ip with respect to wbd is 

^'-u-\rrMf)-'/'h (5.13) 



du-BB 3584 tj|p 

^BD 



2 ^^^NS ^^^BH 



^BD 
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where in the second passage we have included the dependence of the amplitude of 
the waveform ( 15. ip on the mass parameters and in the last passage we have expanded 
M and v with the assumption M^s Mbh- With this behavior one can roughly 
estimate the error A9'^ = a/E^ on the BD parameter. In fact the diagonal element 
of the Fisher matrix T^^^ is quadratic in the derivative f l5.13p . If we invert only this 
element we obtain a bound on cjbd that BBW call un- correlated because this is the 
bound that could be put in the ideal case in which all the parameters were known and 
not correlated among themselves. The bound from T^^^ is explicitly quoted by BBW 
in their tables: it is seen to be always two orders of magnitude higher than the one 
comprising correlations. Here I will only use T^^^ to show the functional dependence of 
the error on the parameters. A part from the frequency dependence in the first line of 
eq. f l5.13p we have 



where we have ignored the dependence on M^s because this mass is kept constant in 
the analysis of BBW. The lower bound on a parameter like wbd can be chosen to be the 
error committed in estimating the parameter itself; denoting by ujbd the lower value 
corresponding to the error (15.14^ one has 



which explains what found by BBW concerning the behavior of the bound cjbd with 
Mbh, i-e. why the most constraining system is a 1.4 Mq neutron star inspiraling into 
a 400 Mq black hole. 

Concerning the dilution of information in the presence of spin, BBW find that the 
spin-orbit term reduces the bound on cjbd by factors of order 10 - 20: for example, in 
the case of the most constraining system, the bound is degraded from wbd ^ 8 x 10^ to 
i^BD ~ 3.9 X 10"^. This last constraint is almost identical to the Solar-System bound from 
Cassini measurements of the Shapiro time delay [32]. However, a bound coming from 
binary coalescences detected at LISA is of greatest importance even if numerically 
equivalent to the constraint derived in a different testing ground. This situation is 
similar to the one concerning the bounds on scalar-tensor theories derived by DEF: 
as I have discussed in Sec. 11.11 the dynamical regimes probed by interferometers are 
very different from those characterizing the Solar-System and binary pulsars. First 
and foremost, the LISA bound obtained by BBW comes from a 2PN template, i.e. 
from considering GR non-linearities that are two orders higher than those probed by 
binary pulsars. In a second place, if one of the binary constituents of a LISA source 
is a neutron star, the companion object will be a black hole of mass Mbh^^^^Mq: 
therefore, the LISA bound would come from the strongest gravitational fields possible. 




(5.14) 
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Constraining the presence of an effective mass term in graviton propagation 

To place bounds on the graviton Compton wavelength, BBW do not fix a thre- 
shold on the SNR. In fact, to investigate this non GR behavior one is not limited to 
consider binary systems which comprise neutron stars, as it was needed to constrain BD 
theory. Therefore, one can take advantage of the high sensitivity of LISA to sources at 
cosmological distances; probing the Universe at such a deeper level, one awaits to detect 
black hole-binaries with constituent masses from W^Mq to WMq. This is indeed the 
mass range investigated by BBW, who considered these sources at a luminosity distance 
Dl = 3 Gpc. 

The derivative of the waveform (15. ip with respect to is 

D M'/' 



oc 



il + z)DL Xl 



V 



where in the second passage we have included the dependence of the amplitude of the 
waveform (15. ip on the mass parameters. The behavior of eq. (15.151) is different with 
respect to eq. (15.131) . notably the derivative is now directly proportional to the total 
mass. For this reason, BBW obtained the best bound on from the system with 
highest total mass, notably an equal mass binary with M = 2 x 10'^ Mq , whose SNR 
is 2063. 

Now I repeat the analysis based on the diagonal element of the Fisher matrix in 
order to estimate the error and understand its dependence on the other parameters. 
From eq. (15.151) we have 



V 1/2 p -1/2 
^AA ~ i AA 

D 



dh 

-1 



{1 + z)Dl 



^ AJM-^/^ (5.16) 



in which we have omitted the dependence on the mass ratio because we are considering 
equal mass binaries. The corresponding lower bound on A^ reproduces the result of an 
earlier investigation by Will [76] : 



where M encodes various contributions from the experimental noise and where the 
dependence on the mass accounts for that of the SNR (cfr. eq. (I5.12p ). As noted in 
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ref. [76], the bound on depends only weakly on distance, via the redshift dependence 
of the factor [-0/(1 + Z)Dl]^^'^: in fact, this ratio varies from unity at 2; = to 0.45 
at 2; = 1.5. The mild dependence of the bound (15.17^ on redshift is due to the fact 
that while the signal strength and the accuracy decrease with distance, the size of the 
arrival-time effect increases with distance. This is connected to the distortion effects 
caused by the presence of an effective mass term in the dispersion relation. 

For all the systems investigated the degradation factor due to the spin-orbit term 
never exceeds 5. The best constraint on the graviton Compton wavelength in the 
presence of spin-orbit is > 2.2 x lO^^km, which is four orders of magnitude larger 
than the Solar-System one [157] and corresponds to A^ > 7 x 10^ pc. It is interesting 
to convert this bound in a constrain on the graviton mass as we did in Sec. II. 4[ In 
units h = c = 1 one has 1 = /ic ~ 200 MeVfm"^, where fm=fermi=10~^^m: with a 
conversion factor of lm=(2 x lO^^ev)"^, one can finally obtain rUg < 10~^^eV. 

5.2.1 The study of gravitational- wave constraints on (3^ and (5^ 
and its comparison with other analysis 

The conclusions of BBW imply that the phasing coefficients corresponding to ubd 
and Xg are both numerically small. In fact, from the bounds discussed just above one 
has 

5.4X10- (5,18) 



84wbd ~ 84(4 X 104) 
while the phasing coefficient due to a massive graviton in the case of LISA reads 

DM {3 X 10^ pc) ■ {3 X W km) 



(AJ2 (2xl0i6A;m)2 

^ (3 X 10^2 A;m) • (3 x W km) 
~ 4 X 1032 km^ 



2 X 10"^ (5.19) 



This characteristic is especially relevant to assess if one can possibly estimate small 
parameters like (S^ and /34 with GW templates comprising spin terms. In this respect, 
another result of BBW is of crucial importance. Both cjbd and Xg contribute a signi- 
ficant number of cycles to the template: from 10 to 1000. This is what really matters 
for GW astronomy: it is the individual value of the PN corrections to the phase, not 
their value relative to a numerically larger term. Therefore, the conclusions of BBW 
suggest that it would be possible to constrain /^a and (3^ even in the presence of spin 
terms. 

It is interesting to compare the work of BBW with the studies of Arun, Iyer, Sathya- 
prakash and collaborators (AISC) |T5l|T6l[T3 that I have discussed in Sec. 11.41 AISC did 
not include spin in their treatment so that a direct comparison between their analysis 
and the one from BBW can not be drawn. However, the results of the methodologies can 
be compared for what concerns the issues of degeneracies and the estimation accuracy 
achievable on individual parameters. In the most powerful version of the test, AISC 
consider a reduced set of PN phasing coefficients and use two of them as a basis to 
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5.2 GW constraints on /^s and /34 in comparison with other analysis 




Figure 5.1: Plot showing the relative errors in the test parameters i/jt = 
4'3,ip4,, ip5i, ipe, ipeh "^7 as a function of the total mass M. The left panel displays the case 
of the third generation ground-based interferometer ET [79] for a stellar mass compact 
binary at luminosity distance Dl = 200 Mpc. The right panel shows the situation 
for a supermassive black hole binary at a redshift of z = 1 as observed for a year by 
LISA [4(I|. The test parameters are analyzed one at a time. Figure taken from ref. [16] . 



express all the others: in so doing, AISC "take advantage" of the correlations among the 
coefficients. One of the conclusions of AISC is that it will be possible to estimate the 
various PN phasing coefficients at future interferometers with the accuracies reported 
in Fig. 15.11 As we can see, the best accuracies are awaited with LISA [40]: even in this 
case the estimation precision is worse than 10~^, apart from one parameter. One could 
wonder how such accuracies might help in distinguishing an error on a GR coefficient 
from the actual presence of, say, that is already known to be smaller than 10~^. 
In answering this question, one should note that the focus of AISC is on estimating 
the PN phasing coefficients ip/s; on the other hand, in BBW one tries to estimate the 
physical parameters of the template, both the intrinsic and the extrinsic ones. The 
difference lies in the fact that a single physical parameter like the mass ratio enters 
almost all the PN coefficients (apart from those that are constant, like tps): therefore, 
the very aims of the two studies are different. As a consequence, for the purpose of 
estimating and P^, the analysis implemented in BBW is the optimal one. On the 
other hand, for the sake of testing gravity in general, the best strategy is to implement 
more than one approach: on top of those I have confronted here, one could employ the 
ones suggested by AISC (see footnote[7]in Sec. II. 4p and the parametrized post-Einstein 
framework proposed by Yunes and Pretorius [80] . 

So far, I have not discussed how the constraints could be affected by a joint ana- 
lysis of many systems. Detection rates for ground-based interferometers will definitely 
improve in passing from the initial configuration to the advanced one, whose imple- 
mentation is already underway. The most up-to-date estimates for detection rates 
have been reported in ref. [H] , from which Table 15.11 is adapted. The work of ref. [H] 
included results appeared until September 2009 concerning various types of compact 
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Interferometer 


Source 


Nlow 






-^max 










yr'i 


yr-1 




NS-NS 


2 X 10^^ 


2 X 10-2 


0.2 


0.6 


Initial 


NS-BH 


7 X 10-^ 


4 X 10-3 


0.1 






BH-BH 


2 X 10-"^ 


7 X 10-3 


0.5 






NS-NS 


0.4 


40 


400 


1000 


Advanced 


NS-BH 


0.2 


10 


300 






BH-BH 


0.4 


20 


1000 





Table 5.1: Detection rates for compact binaries: quantities are quoted in number of 
events per year; the physical significance of the various estimates in the columns is 
explained in Table 15.21 Table adapted from ref . |18] . 



Symbol 


Rate statement 


Physical significance 




Plausible pessimistic estimate 


Rates could be as low as... 


Nre 


Realistic estimate 


Rates are likely to be... 


-^high 


Plausible optimistic estimate 


Rates could be as high as... 


N 

^ ' max 


Upper limit 


Rates should be no higher than... 



Table 5.2: Rate statement terminology. Table adapted from ref. |48] . 



binary coalescences, among which I focus on the following ones: 

neutron star-binaries (NS-NS); 

neutron star-black hole binaries (NS-BH); 

black hole-binaries (BH-BH). 

For a given binary type in a LIGO- Virgo search, the table contains different estimates 
of the detection rate N = R x Nq, where R is the coalescence rate of that type of 
binary per galaxy and Nq is the number of galaxies accessible with a search for the 
relevant binary type. The physical significance of the various estimates is explained in 
Table [El 

For the purpose of calculating the rates in Table 15. H ref. [H] assumed that all 
neutron stars have a mass of 1.4M0 and all black holes have a mass of lOM©. Although 
both types of objects will rather cover a range of masses (see, e.g., |159i 1160] ). the 
present knowledge of the mass distribution is not sufficient to provide detailed models; 
moreover, the uncertainties in the coalescence rates dominate errors from the simplified 
assumptions about component masses. Within these assumptions, ref. |18] estimated 
the event rates detectable by the LIGO- Virgo network by scaling to an average volume 
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5.2 GW constraints on /^s and /S^ in comparison with other analysis 



Interferometer 


Source 


DhorizoniMpc) 




NS-NS 


33 


Initial LIGO 


NS-BH 


70 




BH-BH 


161 




NS-NS 


445 


Advanced LIGO 


NS-BH 


927 




BH-BH 


2187 



Table 5.3: According to the source type, a given interferometer can probe different 
volumes of space, characterized by the so-called horizon distance. With a fiducial 
signal-to-noise threshold of 8, Initial and Advanced LIGO are sensitive to the indicated 
values of the horizon distance [H] . 

within which a single detector is sensitive to coalescences above a fiducial SNR threshold 
of 8. For a single interferometer, this volume probed is characterized by the so-called 
horizon distance: according to the type of source and the interferometer involved in 
the search, the horizon distance assumes the values reported in Table 15.31 

As it is evident from Table 15.11 with the advanced versions of LIGO [78] and 
Virgo |161j one can realistically await as many as 70 detections per year, i.e. more 
than one per week! This large statistics will constitute a major improvement in es- 
timating parameters like and P^. in fact, the values of extrinsic parameters depend 
only on the theory of gravity and not on the system detected; on the other hand, each 
source will be unique for what concerns the intrinsic parameters like the spin orienta- 
tion. Therefore, if N events are detected, each one can be interpreted as providing an 
independent measurement of an extrinsic parameter, for example P^: the total error 
on /Sa would then be a fraction 1/y/N of the error coming from a single observation. 
Detection of many systems has a further advantage in measuring extrinsic parameters: 
as for PPK tests with binary pulsars, there might be systems that are better suited 
for measuring one parameter instead of another. For example, the scalar charge of a 
neutron star is better measurable if one can detect a dipolar flux from an asymmetric 
system like a pulsar-white dwarf binary (see eq. fll.l5p in Sec. II. ip . Moreover, the 
bound derived on a given parameter from one system might be used as a prior when 
testing a different parameter with another source: in this way, one should be able to 
disentangle the effects of various extrinsic parameters. A final note concerns the pos- 
sibility to further increase the accuracy of GW experiments by cross- correlating the 
outputs of many interferometers: this enables one to dig signals lower than the sens- 
itivity of an individual instrument and has proven crucial to obtain the bound on the 
GW background of stochastic origin with first-generation instruments [2]. 

The considerations discussed in the previous sections motivate a dedicated analysis 
to extend the study of Ps and (34 to PN orders in the GW phasing. I have started 



The effects of modified graviton vertices in GW templates 



135 



this investigation thanks to a collaboration with Emanuele Berti, one of the authors of 
BBW. In what follows I describe part of the ongoing work in this direction. We expect 
that the analysis presented in this chapter will be finalized soon |162j . 



5.3 Feynman diagrams for the phase shift induced 
by the modified three- and four-graviton ver- 
tices up to first post-Newtonian order 

As we saw in Sec. 11.41 the starting point to compute the evolution of the orbital 
phase is the energy balance equation ( 11.49^ that we re-write here as 



P 



-E. 



(5.20) 



where P is the power emitted in GWs and E is the mechanical energy of the system. 
The PN expansions of P and E that we met in Sec. 11.41 eq. fll.50p can be conveniently 
dealt with by means of a gauge invariant adimensional quantity 



.2/3 



(5.21) 



where Ug is the frequency of the source bulk movement, for example that of its orbital 
motion; up to IPN, the expressions of P and E in terms of x read 



E = Eox{l + eix) , P = Pox^no + TTix) 



(5.22) 



where ei and tti indicate the IPN corrections. The energy balance equation (I5.20p then 
translates as the time evolution of the variable x, which is given by 



X 



Pq-^o 5 

-X 



E^ 



Til - 27roei 
1 H X 



(5.23) 



and can be easily solved perturbatively {Eq < 0,Po > 0). By redefining the time 
parameter as 



e 



64^0 

the solution of eq. (15.231) reads: 

0-1/4 



(te - t) with (6)"^/^ ~ C»(? 



X 



127rn 



(5.24) 



(5.25) 



Finally, defining the source phase as dcp/dt = Us, one gets the time evolution of the 
orbital phase in terms of the parameter G 



V 



1 + 5 2^1^0 -7ri Q_i/4 
247ro 

+ 02 e-^/^] + const 



const 



(5.26) 
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5.3 Diagrams for the modified non-Abelian vertices of GR up to IPN 





Figure 5.2: The Feynman diagrams that, in NRGR, represent the next-to-leading order 
contributions to the conservative dynamics, i.e. corrections to the Newtonian potential 
which are 0{mGN /r), upper row, and O(f^), lower row, with respect to it. Figure 
taken from ref. [6]. 



where I have introduced the compact notation for the phasing coefficients: 0o for the 
leading order term and 02 for the 0{v'^) IPN correction. As I have discussed in Sec. 11.41 
every theory of gravity has its own prediction for these coefficients: in the case of GR, 
this prediction reads 

^°-"2' '^-"4"12' 

32z/2 1247 35 

Po = —77- , VTo = 1 , TTi = V 5.27 

" 5G ' ' 336 12 ^ ' 

where jU is the reduced mass of the system /i = mim2/M and we remind that Eq and 
Po contribute to the phase through the time parameter 9 (see eq. (15.241) ). Plugging 
this prediction in eq. (I5.26P gives the GR result for the phase as a function of time 



95/8 



V 8064 96 



+ const ; (5.2J 



therefore, in GR one has 



/^3715 55 , 

(8064 + 96" I f^-^^' 



which coincide, respectively, with the coefficients and 02 in eq. (I1.57p . 

Looking at eq. (I5.26p . one can see that the phasing coefficients originate from both 
the conservative and the radiative dynamics of the system; notably, the IPN correction 
to the phase 02 is made of ttq and tti, coming from the expansion of P, and ei coming 
from the expansion of E (see eq. (I5.22p ). In the EFT terminology, these expansions can 
be obtained by means of Feynman diagrams. For what concerns the energy function, 
the IPN correction ei is calculable from the diagrams of Fig. I5.2l that represent the EIH 
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Figure 5.3: NRGR diagram whose imaginary part gives rise to the quadrupolar gravit- 
ational radiation power spectrum (see text and Sec. l2.3.4l for discussion). Figure taken 
from ref . [6] . 



Lagrangian, like in Sec. 12.3.11 On the other hand, the EFT prescription to calculate 
the power emitted in GWs is to take the imaginary part of the self-interaction diagram 
of Fig. 15.31 We remind from Sec. 12.3.41 that this Feynman diagram describes a process 
in the effective theory which is valid at scales i larger than the orbital radius r : for this 
reason, the binary system is un-resolved in the diagram and indicated with a double 
line. In the figure, the individual matter-gravity vertex corresponds to many vertices 
of the theory valid at scales I < r. These vertices are described by the diagrams of 
Fig. 15.41 they are 0{v'^) with respect to the coupling of a radiation graviton with 
the mass monopole of the system (cfr. discussion in Sec. 12.3.41) . In GR both this 
coupling and the coupling to the dipole vanish: therefore, the leading order in the 
radiation dynamics of GR is represented by the 0{v'^) terms of Fig. 15. 4[ Looking at 
these diagrams, one can see that the three-graviton vertex enters the power at leading 
order. A possible deviation of this vertex from its GR value is reflected in a shift in 
the power coefficient ttq : tagging by [i^ the three-graviton vertex, we can keep track of 
this deviation and of its contribution to the phase. The shift produced by /Ss at leading 
order in the power is what enabled my group and me to put the bound (5^ < 10"'^ from 
the data on Hulse- Taylor pulsar. This calculation in the strong-field/radiative regime 
did not need to take into account the shift produced by /^a in the conservative dynamics. 
In fact, as one can see from Fig. 15.21 this last shift affects the EIH Lagrangian, which 
is related to the IPN correction to the energy ei. In other words, the three-graviton 
vertex enters the energy at next-to-leading order: this produces a shift in ei which 
affects the phase (]5.26p at IPN order. In ref. |125j . we rather used the IPN change in 
the conservative part to constrain /^s by means of Solar-System tests. 

For the purpose of estimating (3^ and through GW observations, it is necessary 
to see how the phase is affected by the modified three- and four-graviton vertex at least 
to order IPN, i.e. one has to calculate the modification of the coefficient 02 (15.261) . For 
this one needs the shifts in tti and ei: because I already have the latter, I only have to 
calculate the former. The shift in tti can be obtained from the Feynman diagrams that 
contain three- and four-graviton vertices in the radiative dynamics at IPN. I display 
these diagrams in the following sections. 
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5.3 IPN phase shift induced by the modified three-graviton vertex 



nnnnr 



Figure 5.4: Feynman diagrams describing the leading order (OPN) radiation Lagrangian 
in GR (see text for discussion). 



5.3.1 The phase shift induced by the modified three-graviton 
vertex at first post-Newtonian order 

We start the investigation of the radiative sector from the contributions due to the 
three-graviton vertex. We have to consider which are the possible modifications of 
O(v^) or of 0{Gn) to the three-graviton vertex diagram of Fig. 15.41 

A first type of modifications comes from the PN corrections of the matter-graviton 
vertices with factors of the velocity: these come from the expansion of the point particle 
action 



5, 



m 



pp 



Ml 



PI 



dt 



~-Hoo — HoiVi — -i^oov^ 



(5.30) 



and produce the diagrams of Fig. 15.5 



1 + ps) nrmnr 



(a) 



(b) 



Figure 5.5: IPN modifications of the three-graviton vertex coming from the velocity 
expansion in the point-particle action eq. fl5.30p . 



In the spirit of the PN expansion, one treats time derivatives as perturbations. For 
what concerns the EFT approach we met an example of this attitude in Sec. 12.3.31 
when we were dealing with the propagator of a potential graviton (see the derivation 
of eqs. f l2.46p and fl2.55p ). We report excerpts from that discussion here for a more 
complete treatment. By considering potential modes in Fourier space one can expand 
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ID '2D 



02D 



(a) (b) 



Figure 5.6: IPN modifications of the three-graviton vertex coming from the inclusion 
of time derivatives. 



their propagator as: 

1 1 



k2 ^ 



^ [l + 0{v')] , (5.31) 



k2 



whose leading order term k ^ gives the 1/r dependence of Newtonian potential once 
Fourier-transformed according to 



d'k 1 



re 



5{x^) . (5.32) 



(27r)4 k2 47r|x 
Therefore, to leading order, the propagator for ifk/xi/ reads 

(Ti7kM.(a;°)^qa/3(0)) = - j^(27r)35=^(k + q)<5(a;°)P^,;„^ , (5.33) 

where T stands for time- ordering and the tensor structure is given by eq. ( 12.17^ . On 
the other hand the 0{v'^) terms give a propagator which reads 

(i/kM.(^°)^qa/3(0))® = -(27r)353(k + q)^5(xO)P^,^„^, (5.34) 

where the subscript is to distinguish this propagator from the one in eq. fl5.33p and 
where we stress that 



(5.35) 



k4 k4 k4 k2 ' 

so that eq. f l5.34p is indeed suppressed by with respect to eq. fl5.33p . 

In the present context, time derivatives are to be considered also in graviton self- 
interactions like the hHH vertex of eq. f l2.82p in Sec. 12.3.41 The resulting diagrams are 
those of Fig. (15. 6p . where the blobs and the acronyms 1-D, 2D indicate corrections due 
to the inclusion of one or two time derivatives, respectively. 

There is another way of getting PN modifications from derivatives that is peculiar 
to radiative dynamics: it is the multipole-expansion of the radiative field at the level of 
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5.3 IPN phase shift induced by the modified three-graviton vertex 



X X 



Figure 5.7: IPN modification of the three-graviton vertex coming from the multipole- 
expansion of the radiative field, a necessary ingredient of the EFT treatment. 



/ \(i + A) 



Figure 5.8: IPN contribution to the radiative sector coming from 0{Gn) modifications 
of the diagrams in Fig. (15.41) . 

the action, a necessary ingredient of the EFT treatment [96] . Condensing the discussion 
of Sec. I2.3.4| one can take the origin of the reference frame to be the center of mass of 
the system and write the radiative field as 

V(x, t) ~ V(0' t) + x'd,h^u{0, t) + ^x'x^did,h^,{0, t) , (5.36) 

where x*'s are coordinates with respect to the center of mass and where the (spatial) 
derivative scales as di ~ v/r so that x^di ~ v. To have a correction of 0{v'^) we need 
to consider the last term of eq. ( 15.36p . as it is the case for the diagram in Fig. ( 15. 7p . 

Finally, we have the 0{Gn) modifications due to graviton self-interactions. A first 
example is the one represented in Fig. (15. 8p . 

After Fig. (15. Sp it is easy to build the diagrams of Fig. (15. 9 p and Fig. (I5.10p . 
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1 



2 




(a) (&) 

Figure 5.9: IPN contributions to the radiative sector coming from 0{G]\[) modifications 
of the diagrams in Fig. (15. 4p . 



(1 + /93) 



.\(l + /?3) 



(a) (b) 

Figure 5.10: IPN contributions to the radiative sector coming from 0{G]\f) modifica- 
tions of the diagrams in Fig. (15. 4p . 
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5.3 IPN phase shift induced by the modified four-graviton vertex 



5.3.2 The phase shift induced by the modified four-graviton 
vertex at first post-Newtonian order 

Because of the scahngs of graviton self-interactions, a diagram containing a four- 
graviton vertex is of order IPN with respect to diagrams where a three-graviton vertex 
is present. Therefore, the four-graviton vertex contributes for the first time to the 
phasing directly at IPN order through the diagram of Fig. (15. lip . 



xSLSLSiJLy 



Figure 5.11: The only diagram containing a four-graviton vertex that contributes to 
the radiative sector at order IPN. 



5.3.3 The phase shift induced by the modified three-graviton 
vertex at half post-Newtonian order 



1 



2 



•ID 



(a) (6) 

Figure 5.12: Feynman diagrams representing possible modifications of a three-graviton 
vertex at half PN order in the phase. (See text for complete discussion.) 

With scaling arguments in mind, one can also build the diagrams of Fig. 15. 121 which 
are 0{v) with respect to the quadrupole. In GR both these diagrams do no contribute. 
Let us see in more detail why this is so. 

We start from the diagram in panel (a). The velocity factor attached to the worldline 
is contracted with a field H^i, which in turn requires a couple of indices Oj in the 
three-graviton vertex, otherwise the propagator vanishes. There are two possible ways 
to further contract these indices. The first possibility is to have kokj i.e. one time 
derivative in the vertex; however, this would bring another factor of order v. The 
diagram in Fig j5. 12( a) would then be incomplete as it stands and would not belong to 
a 0.5PN order: if one draws the ID symbol for the derivative, the result is one of the 
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IPN diagrams we met in the preceding section. A second possible way of contracting 
the indices Oj in the vertex is offered by the polarization of the radiative graviton: this 
mode is non-radiating in GR to obey a tenet of the theory, so it makes sense to consider 
it in presence of non- vanishing and P4. 

Now we consider Fig. 15.12( b). A time derivative in the vertex amounts to having 
a term kokj which we can contract with a radiative h^j mode: this mode is forbidden 
in GR but not necessarily in alternative theories. Incidentally, we remind here that a 
similar case of terms which are forbidden in GR allowed us to put a constrain on 
(see the discussion at the end of Sec. 14. 3[ starting from eq. fl4.16p ). 

At the same 0.5PN order there is also a term from the multipole expansion (cfr. 
Fig. 15. 7p . The bottom line is that the presence of a shift in the three-graviton vertex 
contributes terms of order 0.5 PN in the phase on top of those of order IPN described 
in the previous sections. Therefore, the full set of diagrams that should be investigated 
at order 0.5 PN in the phase are those of Fig. I5.13[ 



v'H, 



Oi 



hoj jTTinnr h kj /lo, ^^nnnr da nnnrrr 



(a) {b) (c) 

Figure 5.13: Terms of half-PN order in the phase excited by the three-graviton vertex 
modification P^. The 0.5PN order is absent in GR. 



5.4 Final remarks 

In this chapter I have presented an extension of the work I performed with the 
rest of the gravity group at the University of Geneva [125j : there we proposed a field- 
theoretical extension of the PPN, where the parameters are used to tag the strength of 
graviton vertices, notably the deviation of the three- and four-graviton vertex from the 
GR predictions. With available data from experiments in the Solar System and from the 
timing of binary pulsars, we could constrain these deviations to interesting accuracies. 
Concerning the bounds that could be put with future detections at interferometers, we 
qualitatively concluded that the numerical smallness of Ps and P4 with respect to other 
parameters in the waveform at the same PN order could have prevented us from getting 
significant constraints. For this reason we decided not to go further in calculating the 
PN modifications that Ps and introduce in the phase of GWs and to publish our 
results. Because we worked in the restricted PN framework, we suggested that these 
degeneracies could be broken by the inclusion of amplitude corrections in the template. 
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5.4 Final remarks 



as higher harmonics enrich the functional dependence of the parameters [156] . However, 
a richer functional dependence could offer the possibility of binding small parameters 
like and Pi even in the restricted PN framework. An example along this direction is 
the work of ref . |132j , where Berti, Buonanno and Will investigated the bounds on some 
alternative theories that could be put with LISA observations if templates include spin- 
orbit and spin-spin effects. I have described this analysis in Sec. 15. 2[ focusing on the 
aspects which are more salient to the problem of estimating and P^; in this context, 
the most notable conclusion is the following. Deriving the Fisher matrix, Berti et 
al. found that the inclusion of spin-orbit contributions in the waveform degrades the 
accuracy in estimating the parameters of alternative theories; however, the possibility 
of deriving meaningful constraints on these parameters in the realistic case of spinning 
objects does not get spoilt. On these grounds, it seems likely that one can constrain 
P^ and Pi even in the presence of spin terms. In collaboration with Emanuele Berti I 
am conducting this follow-up of the work of ref. |125j . The details of this investigation 
have been discussed in Sec. 15. 3[ where I have presented the Feynman diagrams needed 
to calculate how the modified three- and four-graviton vertices affect the phase up to 
IPN order. Once the phase modification is obtained, an analysis of the type of ref. [132] 
should enable Berti and me to derive the constraints that future GW observations will 
put on the parameters P^ and P^ |162j . 

The field-theoretical extension of the PPN that characterizes refs. |125t 1162] sets it- 
self in the active research domain of testing gravity with GWs. In sections 11.41 and 
15.2.11 I have discussed two other phenomenological approaches that have appeared 
in the literature to parametrize the gravitational waveform: the series of studies by 
Arun et al. [151 [IS 113 and the parametrized post-Einstein framework put forward by 
Yunes and Pretorius in ref. [80]. These studies and the work discussed in the present 
chapter constitute relevant phenomenological frameworks to test gravity in the strong- 
field/radiative regime represented by GWs. 



Conclusions 



This thesis is based on the effective field theory approach that has been recently 
proposed to tackle the two-body problem in General Relativity. Such a framework 
has been borrowed from particle physics studies of non-relativistic bound states: for 
this reason it has been dubbed Non-Relativistic General Relativity (NRGR) by their 
authors, Goldberger and Rothstein. Through NRGR one can reformulate the post- 
Newtonian expansion of binary dynamics in a streamlined fashion. In fact, an effective 
field theory approach is best suited to treat a problem like the inspiral of compact 
binaries, where the disparity of physical scales allows one to separate the relevant 
degrees of freedom in subsets with decoupled dynamics. At each scale, by means 
of power counting rules, one can assess the number of terms to be included in the 
computation of a given observable: in a field-theoretical context, this perturbative 
calculation is conveniently performed by means of Feynman diagrams. 

The use of these diagrammatic techniques is quite natural when discussing alternat- 
ive theories of gravity with respect to General Relativity, as first shown by Damour and 
Esposito-Farese. A major example in this sense is represented by scalar-tensor theories 
where the dynamics of the helicity-2 gravitons characterizing General Relativity is sup- 
plemented by a sector of hehcity-0 modes. In collaboration with Riccardo Sturani of 
the gravity group at the University of Geneva, I apphed NRGR to both point-like and 
string-like sources in a scalar-tensor theory of gravity. Notably, wc treated the issue 
of the renormalization of the energy-momentum tensor and reproduced two types of 
results which had been obtained previously: with the insight of NRGR we were able to 
re-derive and confirm the resolution of an apparent discrepancy between those results. 

Beside representing a convenient computational tool, diagrammatic techniques en- 
able one to shed a different light on the information extracted from experiments of 
relativistic gravity. This field-theoretical perspective is the standpoint taken in the 
present thesis and it has been exposed in Chapters 4 and 5. There I have reported 
investigations which I developed in two steps: in a first stage I worked with the en- 
tire gravity group at the University of Geneva composed of Stefano Foffa, Michele 
Maggiore, Hillary Sanctuary, Riccardo Sturani and me; afterwards, I started a collab- 
oration with Emanuele Bcrti at the University of Mississippi. In these works we adopt 
NRGR description of non-linearities in General Relativity: these are represented by 
a perturbative series of Feynman diagrams in which classical gravitons interact with 
matter sources and among themselves. Focusing on the purely gravitational sector, we 
tag self-interaction vertices with parameters: in such a way, my group and I were able, 
for example, to translate the measure of the period decay of Hulse- Taylor pulsar in 
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a constraint on the three- graviton vertex at the 0.1% level. What emerges from this 
attitude is an extension of the parametrized post-Newtonian framework of Will and 
collaborators to the strong-field/radiative regime of gravity. Such a formulation is very 
useful in view of gravitational wave detection, which will constitute a unique source 
of knowledge in the domains of gravity, in particular, and physics, in general. In this 
context, our parametrized field-theoretical framework constitutes a complementary al- 
ternative to other approaches in the field: the reference example is the series of studies 
put forward by Arun, Iyer, Sathyaprakash and collaborators to test General Relativity 
by measuring the coefficients of the gravitational wave phase. 



Abbreviations List 



BD = in Chap. [3l Buonanno & Damour, to refer to their work on the effective action 
of cosmic strings [110] . in Chap. |5l Brans & Dicke, to refer to their scalar-tensor 
theory [26] 

BH = Black Hole 

BBH = Binary Black Hole 

DEF = Damour and Esposito-Farese, as in their treatment of scalar-tensor theories, 
described in Chap. [1] 

DH = Dabholkar and Harvey, to indicate their study |109j discussed in Chap. [3] 

EEP = Einstein Equivalence Principle, i.e. universality of free-fall together with 
local Lorentz invariance and local position invariance (see Will's review [19] for 
a comprehensive discussion) 

EFT = Effective Field Theory 

EH = Einstein-Hilbert (referred to the action) 

EIH = Einstein-Infeld-Hoffmann (referred to the action) 

EMT = Energy-Momentum Tensor 

ET = Einstein Telescope ^79] 

GR = General Relativity 

GW(s) = Gravitational Wave(s) 

IR = Infra Red 

LIGO = Laser Interferometer Gravitational-wave Observatory [38] 
LISA = Laser Interferometer Space Antenna |40J 

LLR = Lunar Laser Ranging, referred to the experiment which uses retroflecting 
mirrors left by the first men on the Moon to monitor the distance of the Earth 
from its natural satellite (see for example ref. [29] ) 

NRGR = Non Relativistic General Relativity, which is how the EFT of ref. [6] has 
been coined by the authors in analogy with the treatment of two-body problems 
in QED and QCD 

NS = Neutron Star 

PK = Post-Keplerian 

PM = Post-Minkowskian 
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PN = Post-Newtonian 

PPK = Parametrized Post-Keplerian, referring to the phenomenological approach 
proposed in ref. [50] to analyze pulsar data in a theory-indipendent way 

PPN = Parametrized Post-Newtonian, referring to the original phenomenological 
framework used to interpret Solar-System experiments on gravity (see ref. [27j 
for a comprehensive treatment) 

PSR = Pulsar 

SEP = Strong Equivalence Principle, i.e. the version of the principle which is valid 
for self-gravitating objects 

UV = Ultra Violet 
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